Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,918,022 A→G 6.2% intergenic (+209/+23) barA → / ← gudD hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,918,0220AG 6.2% 715.9 / 20.5 219intergenic (+209/+23)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base A (101/104);  new base G (5/9);  total (106/113)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

Only 100 of 224 total aligned reads displayed.

TGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAA  >  NC_000913/2917752‑2918276
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                              
tgctgGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATagag                                                                                                                                                                                                                                                             <  1:67032/274‑1 (MQ=255)
                          aaTGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                                                                                                                                  <  2:603473/275‑1 (MQ=255)
                                                             aTGGGTTTCTTTTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGGCAAAGTGCCTGTGGTTATGCCGGATGCGGCGTTAACGCCTTAGCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCGGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                                                                                                                                                <  4:118556/274‑1 (MQ=255)
                                                                        gTCGAGTGGCAATCGGATTTTTATGACGCTGGTTTCAGATTGGTGACAAGGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                                                                                                                                              <  2:305523/217‑1 (MQ=255)
                                                                           gggTAGAATTCGGCTTTTTATGAGCCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTg                                                                                                                                                                                 <  4:570207/273‑1 (MQ=255)
                                                                                                                 ggTGACAAAGTGCC‑‑TGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                         >  2:326468/1‑275 (MQ=255)
                                                                                                                            gCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATaga                                                                                                                                                                                                                                                              >  5:736300/1‑149 (MQ=255)
                                                                                                                                gtgtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCg                                                                                                                                                                                                                                                          >  7:1050490/1‑149 (MQ=255)
                                                                                                                                   tttATGCCGGATGCGGCGTTAACGCCTTATCCGGGCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCTATCTGGAGCCggtt                                                                                                                                                                                                                                                       <  6:470124/149‑1 (MQ=39)
                                                                                                                                   tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGGTGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCTATCTGGAGCCggtt                                                                                                                                                                                                                                                       >  8:331108/1‑149 (MQ=38)
                                                                                                                                   tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGGTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggtt                                                                                                                                                                                                                                                       >  7:269644/1‑149 (MQ=255)
                                                                                                                                   tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCTATCTGGAGCCggtt                                                                                                                                                                                                                                                       <  6:96038/149‑1 (MQ=255)
                                                                                                                                   tttATGCCGGAGGCGGCGTTAGCGCCTTATCCAGCCTACAAAACCAATTAGATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggtt                                                                                                                                                                                                                                                       >  8:926038/1‑149 (MQ=255)
                                                                                                                                    ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCTATCTGGAGCCggttt                                                                                                                                                                                                                                                      <  7:331108/149‑1 (MQ=39)
                                                                                                                                     tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggtttt                                                                                                                                                                                                                                                     >  7:678281/1‑149 (MQ=255)
                                                                                                                                      aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCTATCTGGAGCCggtttga                                                                                                                                                                                                                                                    <  5:70181/149‑3 (MQ=39)
                                                                                                                                      aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggttttg                                                                                                                                                                                                                                                    <  8:439363/149‑1 (MQ=255)
                                                                                                                                      aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggtttcg                                                                                                                                                                                                                                                    <  6:770994/149‑3 (MQ=255)
                                                                                                                                       tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTGAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑TGACACGCTACAGATAGAGCCggttttgg                                                                                                                                                                                                                                                   >  8:169377/1‑149 (MQ=38)
                                                                                                                                       tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggttttgg                                                                                                                                                                                                                                                   <  5:885392/149‑1 (MQ=255)
                                                                                                                                       ggCCGGATAAGGCGTTAACGCCGCATCCGACCTACAAAACCAATTAAATGCAATGACTCAAAAAATGATGTATACCTGATACATAGCACAAAAGGCTATTTTGCGATTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggctttgg                                                                                                                                                                                                                                                   >  8:1060468/2‑149 (MQ=255)
                                                                                                                                        gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggttttggt                                                                                                                                                                                                                                                  <  5:237631/149‑1 (MQ=255)
                                                                                                                                          cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                                                                                                                                                >  7:194663/1‑149 (MQ=255)
                                                                                                                                             ataaggcgTTAACGCCGCATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                             >  7:853303/5‑149 (MQ=255)
                                                                                                                                              tGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                             >  8:848781/1‑148 (MQ=255)
                                                                                                                                                 ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt                                                                                                                                                                                                                                         >  8:129397/1‑149 (MQ=255)
                                                                                                                                                  gcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                                                                                                                                                >  4:716789/1‑141 (MQ=255)
                                                                                                                                                  gcgTTAAAGCCTTATCCGGCCTACAAAACCAATTAAATTAAATGAATTAAAAAATGATTTAGACCGAGTTCATAGAACAACAGGCAATTTTGGGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGTGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACGGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGTTACAGCTCATGGGCTTTCATCACTTGAt                                                                                                          <  4:447932/275‑1 (MQ=255)
                                                                                                                                                    gTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTaa                                                                                                                                                                                                                                      <  5:955943/149‑1 (MQ=255)
                                                                                                                                                    gTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTaa                                                                                                                                                                                                                                      <  5:584734/149‑1 (MQ=255)
                                                                                                                                                    gTTAACGCCCTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGCTAGAGCCGGTTTTGGTTTTCTGTCTTat                                                                                                                                                                                                                                      >  8:694439/1‑148 (MQ=255)
                                                                                                                                                      tAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                                                                                                                                            <  2:716789/141‑1 (MQ=255)
                                                                                                                                                          gCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAcc                                                                                                                                                                                                                                >  5:535253/1‑149 (MQ=255)
                                                                                                                                                             ttATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCtt                                                                                                                                                                                                                                        >  1:45266/1‑138 (MQ=255)
                                                                                                                                                                  cGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAAACCGGCTCTATCTGGAGCCggttt                                                                                                                                                                                                                                                      >  5:567037/1‑120 (MQ=255)
                                                                                                                                                                   ggCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACggg                                                                                                                                                                                                                       >  5:259805/1‑149 (MQ=255)
                                                                                                                                                                      cTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                                                                                                                                    <  5:782763/149‑1 (MQ=255)
                                                                                                                                                                          aaaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGt                                                                                                                                                                                                                >  8:793497/1‑149 (MQ=255)
                                                                                                                                                                                aaTTAAATTCAATTAATTAAAAAAGGATGTAGACCTGACACATAGCACAACAGGCTTTTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGTTAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGtt                                                                            <  4:326468/275‑1 (MQ=255)
                                                                                                                                                                                 aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                                                                                                                                  >  2:446508/1‑204 (MQ=255)
                                                                                                                                                                                  ttAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        <  6:645477/149‑1 (MQ=255)
                                                                                                                                                                                  ttAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        >  7:575957/1‑149 (MQ=255)
                                                                                                                                                                                       ttCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTcca                                                                                                                                                                                                   >  6:1048171/1‑149 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGt                                                                                                                                                                                                                                           >  4:197570/1‑106 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAt                                                                                                                                                                                                             >  5:109203/1‑136 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        >  7:454906/1‑141 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        <  8:454906/141‑1 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                                                                                                                                                >  6:895651/1‑149 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                                                                                                                                                >  5:409561/1‑149 (MQ=255)
                                                                                                                                                                                           gtttatcagagactgattTAGACCTGTNANATATCACAACAGGCTTTTTTGCGTTTCGCCTCAGTCTC‑AAACCGGNCCCAGATAGCGCCTGTTTTGGTTTTCTGTCGTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCCNCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCAc                                                                  <  3:60845/261‑1 (MQ=255)
                                                                                                                                                                                           aTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcca                                                                                                                                                                                               <  7:299992/149‑1 (MQ=255)
                                                                                                                                                                                            tGAATTAAAAAATGATGTAGACCTGGTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                                                                                                                              <  8:269644/149‑1 (MQ=255)
                                                                                                                                                                                            tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                                                                                                                              <  6:16599/149‑1 (MQ=255)
                                                                                                                                                                                              aaTTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTa                                                                                                                                                                                                                                       <  2:197570/106‑1 (MQ=255)
                                                                                                                                                                                                ttAAAAAATGATGTAGACCCGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAAt                                                                                                                                                                                          <  8:400650/149‑1 (MQ=255)
                                                                                                                                                                                                 tAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATc                                                                                                                                                                                         <  8:655802/149‑1 (MQ=255)
                                                                                                                                                                                                 tAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATc                                                                                                                                                                                         <  6:507735/149‑1 (MQ=255)
                                                                                                                                                                                                 tAAAAAACGATGTAAACCTGATACATAGCACAACAGGCTATTTTGCGTTTATCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTAAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGaaac                                                                                                                                                                                         <  6:574793/149‑3 (MQ=255)
                                                                                                                                                                                                      aaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATa                                                                                                                                                                                    <  6:722742/149‑1 (MQ=255)
                                                                                                                                                                                                             tAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCAtt                                                                                                                                                                             <  5:919801/149‑1 (MQ=255)
                                                                                                                                                                                                               gACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        >  7:827606/1‑120 (MQ=255)
                                                                                                                                                                                                                  cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcat                                                                                                                                                                        >  5:552384/1‑149 (MQ=255)
                                                                                                                                                                                                                  cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcat                                                                                                                                                                        >  6:702283/1‑149 (MQ=255)
                                                                                                                                                                                                                   tGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatc                                                                                                                                                                       <  6:409561/149‑1 (MQ=255)
                                                                                                                                                                                                                     aTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgc                                                                                                                                                                     <  6:818036/149‑1 (MQ=255)
                                                                                                                                                                                                                     aTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCCCTTGATc                                                                                                         >  1:192009/1‑209 (MQ=255)
                                                                                                                                                                                                                       aCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGc                                                                                                                                                                                                                            >  6:586070/1‑92 (MQ=255)
                                                                                                                                                                                                                        cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcatc                                                                                                                                                                  <  5:838945/149‑1 (MQ=255)
                                                                                                                                                                                                                        cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcatc                                                                                                                                                                  <  5:357544/149‑1 (MQ=255)
                                                                                                                                                                                                                        cATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcgcatc                                                                                                                                                                  <  7:140095/149‑1 (MQ=255)
                                                                                                                                                                                                                         aTAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                                                                                                                                                <  2:687255/118‑1 (MQ=255)
                                                                                                                                                                                                                           aGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTc                                                                                                                                                               <  8:218089/149‑1 (MQ=255)
                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATANTGCATTCCCATCGCATCGTCACCNCCGCNAAGCCCCTGTTTCTGATACAGCTCATTGGCTTTCCTCACTTGATCCataaccc                                                                                                 >  4:653163/1‑202 (MQ=255)
                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgcg                                                                                                                                                          >  6:942736/1‑149 (MQ=255)
                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgcg                                                                                                                                                          >  6:634720/1‑149 (MQ=255)
                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgcg                                                                                                                                                          >  5:1016099/1‑149 (MQ=255)
                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgcg                                                                                                                                                          >  7:408000/1‑149 (MQ=255)
                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                                                                                                                                        >  1:152238/1‑167 (MQ=255)
                                                                                                                                                                                                                                 acaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCt                                                                       <  2:100627/231‑1 (MQ=255)
                                                                                                                                                                                                                                  caGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTaa                                                                                                                                                                                                                                      >  1:108848/1‑71 (MQ=255)
                                                                                                                                                                                                                                     gCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCa                                                                                                                                                     >  5:473763/1‑149 (MQ=255)
                                                                                                                                                                                                                                      cTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                                                                                                                                  <  3:108848/71‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTTTCCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGGGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGgcgc                    <  2:530188/275‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTc                                                                                                                                         <  5:406888/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGATTCTGa                                                                                                                                      <  6:781582/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATa                                                                                                                                    >  7:126447/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTCATACAg                                                                                                                                 <  6:535253/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           tctc‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCt                                                                                                                               >  6:143979/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             tc‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTccc                                                                                                                                                                          >  5:714839/1‑104 (MQ=255)
                                                                                                                                                                                                                                                              caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCa                                                                                                                             >  8:729819/3‑149 (MQ=255)
                                                                                                                                                                                                                                                               aaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                                                                                                                                      >  7:445986/2‑139 (MQ=255)
                                                                                                                                                                                                                                                                 aaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATgg                                                                                                                          <  7:926038/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATgg                                                                                                                          >  8:881149/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      gCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACt                                                                                                                                                                                  >  3:174442/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                      gCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCt                                                                                                                       >  4:206460/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                       cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtc                                                                                                                    >  6:382354/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtc                                                                                                                    <  5:677546/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                   >  8:619481/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        tCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                   >  7:105667/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                   <  2:206460/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                              
TGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAA  >  NC_000913/2917752‑2918276

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑