Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 23.2% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC23.2% 345.5 / 49.3 253intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (79/111);  new base C (8/55);  total (87/166)
Fisher's exact test for biased strand distribution p-value = 1.48e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164433‑164727
                                                                                                                                                     |                                                                                                                                                  
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                   >  2:8881506/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                             >  1:7645810/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                             >  1:10607746/1‑150 (MQ=255)
            aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                       >  2:9828039/1‑150 (MQ=255)
              agagATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAATTACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                                     <  2:1200912/150‑1 (MQ=39)
                 ttttaaggggCTTTACCCCAAACTGGGCGGGCCGGACGCCCCGCAAAATAGTG‑CACCGACGCGACGGACTAAAATTTTTTGTTATGTTTTTTTTTCCCCCTCCCCGGGGGGGTGTGGGCGGGGGTTTGGGCATCAGCGCGCACGTtcacc                                                                                                                                  <  2:9390445/143‑1 (MQ=255)
                   tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggccca                                                                                                                                >  2:5771373/1‑143 (MQ=255)
                    gAATGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGAGGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                               <  2:6747533/150‑1 (MQ=255)
                     aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                              <  2:10063433/150‑1 (MQ=255)
                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATCCTG‑CCCCGACGCGACGGACGAAAAAGTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                            <  2:2561919/150‑1 (MQ=255)
                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                            <  2:6887150/150‑1 (MQ=255)
                          gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccacccccga                                                                                                                         >  2:1676694/1‑148 (MQ=255)
                            tttATCCCAAACTGGGGGGGCGGGACGACCGGGAAAATACTG‑CACCGACGCGAGGGACTAAAAAGTATTGTTATGGTTGGTTTTCTCCCTCCCCCTGTGGTTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                       <  2:9790605/149‑1 (MQ=255)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGTGACcccaaccccaac                                                                                                                       >  2:373265/1‑150 (MQ=38)
                             ttATCCCAACCAGGTCTGGCCGGCCGACCGGGCAAATACTG‑CACCGCCGCGCCGGACTAAAAATTTTTCGTACGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:7899557/150‑1 (MQ=25)
                             ttATCCCAAACTTGTCGGGCCGGACGACCCGGCAAATCCTG‑CACCGACGCGACGGACGAAAAAGTATTCTTACGGTTGGTTTTCGCCCACCCCCGGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:6976897/150‑1 (MQ=255)
                             ttAGCCCAAACTGGTCTGGCCGGACGACCGGGCAAATACTG‑CACCGCCGCGAGGGACGAAAATGTTTTCGTACGGTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:5787398/150‑1 (MQ=37)
                                  ccAACCATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                 <  2:11009585/150‑1 (MQ=255)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCGaagacacta                                                                                                                 >  2:5469871/1‑143 (MQ=255)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCACACCCGAAcccacac                                                                                                                 >  2:3186467/1‑146 (MQ=38)
                                    aaaCATGTCTGGCAGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACAACCAAAaccacacac                                                                                                               >  1:6623919/1‑144 (MQ=37)
                                       ctggTCGGGCCGGACCCCCCGGAAATCCCGG‑CCCCGACGCGCGGGCCGAAAAAGTTTTGGTTGGGTGGGTTTTCTCCCCCCCCCGGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                            <  2:2028699/147‑1 (MQ=255)
                                        aTGTCGGGCCGGACGACCCGGAAAATACTC‑CACCGACGCCACGCACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                           <  2:7989721/150‑1 (MQ=38)
                                         tggcgggCCGGACGCCCCGGAAAATACTG‑CACCGACGCGACGGACTAAAATGTATTCGTATGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                          <  2:8711731/145‑1 (MQ=255)
                                           tCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTCGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAACCCCCACCCACACggg                                                                                                        >  2:4793137/1‑150 (MQ=37)
                                              ggCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACTAAAATTTTTTGGTATGGTTGTTTTTCTCCCTCGCGCTGTGGTCGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  1:9157925/150‑1 (MQ=255)
                                              ggCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCACCCTCTCCCGGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  2:821401/150‑1 (MQ=255)
                                               gcgggccgacCCGGCAAATCCTG‑CACCGCCGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgagat                                                                                                    <  2:5423386/144‑2 (MQ=255)
                                               gCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCCAACCCACACCCACaaccgggag                                                                                                    >  2:10358322/1‑143 (MQ=37)
                                                cgggccgacCC‑GCAAAAACCGCCACCGAGGCGCCGGACGAAAATGTTTTCGTATGTTGGGTTTTCGCCCCCCCCGGGGGGGTTTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                   <  1:4417494/145‑1 (MQ=255)
                                                ccGGACGCCCCGCAAAATCCTC‑CACCGCCGCGCCGCACTAAAATTTTTGCTTTTGGTTGGTTTCCCCCCTCTCCGTGTGGGAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                   <  1:6398284/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAACCACACCCCCCAgggggggg                                                                                                   >  1:1202921/1‑145 (MQ=25)
                                                  ggtcgccccGCCAAATACGG‑CCCCGCCGCGCCGCACTAAAATGTTTTCTTATGGTGGGTTTACTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggt                                                                                                 <  1:10416188/144‑2 (MQ=255)
                                                     cgacCCGGCAAATACGG‑CCCCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCGCCCCCTCCCGGTGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                              <  2:10717633/150‑1 (MQ=255)
                                                       tcccGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCACCCCCAAACCCACGCCCACCCGGGCGGGGGGGCCgc                                                                                            >  1:1802127/2‑149 (MQ=25)
                                                       acCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCCCACCCGCACCCACACGGGGGCGGGGGcccg                                                                                            >  1:8600067/1‑148 (MQ=25)
                                                         ccGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCACAAGGGGGCGGGGGCCGCTc                                                                                          >  1:10193711/1‑150 (MQ=37)
                                                         ccGGCAAAACCAC‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  2:9692227/150‑1 (MQ=38)
                                                         ccGGCAAAAACTG‑CACCGACGCGAGGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  2:6451638/150‑1 (MQ=255)
                                                          cGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACTCCCTCAAGGGAGAGGGGGCCGATCg                                                                                         >  2:7363665/1‑150 (MQ=38)
                                                          cGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAACCACAACCCCAACGGCGAGGGGGcccaaag                                                                                         >  2:4322513/1‑145 (MQ=25)
                                                            gaaaaTCCTG‑CCCCGACGCGCCGGACTAAAAATTATTCGTAAGGTTGGTTTACTCCCTCGCCCTGTGGTAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  1:6783149/148‑1 (MQ=255)
                                                            gCAAATACTG‑CACCGACGCGACGGACGAAAAAGTGTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCACAACCCCACCCCCACCGGGGAGGGGGCCgggaggg                                                                                       >  1:2340931/1‑144 (MQ=21)
                                                                  aCTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCCAAACCGCACCCACAAGGGGGAGGGGGCCGATAGAGCCCaaa                                                                                 >  1:5416018/1‑150 (MQ=25)
                                                                  aCTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAACCCACACCCCCACGGGGGAGGGGCCCGCTCGTGCACaaa                                                                                 >  1:7103696/1‑150 (MQ=25)
                                                                       cACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCCCCCGCACCCGAACCCACACCCCCCCGGGGGGGGGGCCCCGTCGGGCCCCAAtttt                                                                             >  2:7298613/1‑150 (MQ=21)
                                                                         ccgccgCGCCGTACGAAAAAGTATTCTTCCGGTGGGTTGTCCCCCGCGCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:4262030/150‑1 (MQ=255)
                                                                         ccgccgCGAGGTACTAAAAAGTATTGGTAGGGTGGTTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:4298478/150‑1 (MQ=37)
                                                                         ccgccgCGACGGCCGAAAAATTATTGGTAGGGTTGGTTTTCTCCCTCCCCGGGTGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:5231854/150‑1 (MQ=255)
                                                                         ccgccgCGACGGACTAAAATGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:10866627/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACTAAAAAGTTTTCGTCCGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:9889806/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACTAAAAAGTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:3525055/150‑1 (MQ=255)
                                                                         ccGACGCGACGTACTAAAATTTATTCTTATGGTGGTTTTCCGCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:6938040/150‑1 (MQ=37)
                                                                         ccGACGCGACGGACTAAAACTTTTTCGTAAGTTTGGTTTTCTCCCTCCCCCTGGGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:7379325/150‑1 (MQ=38)
                                                                         ccGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:11093008/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCAAACCCAAAACCACACACACAAGGGGGAGGGGGCCGAAAGAGCCCAAAatttta                                                                           >  1:9098699/1‑148 (MQ=21)
                                                                              gcgacgGACTAAAAAGTATTCGTAAGGTTGGATTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  1:9973545/150‑1 (MQ=255)
                                                                              gcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  2:10056570/150‑1 (MQ=255)
                                                                               cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAACAGCGCGCACGTTCACCCGCACCCGAACCCACTACCGCAAGGGCGGGGGGCCCGATCGAGCACAAAATTTgagagaaa                                                                     >  2:9164962/1‑148 (MQ=37)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCCCACCCAAAACCGCAACCACCAGGGAGGGGGGGCCCGGCGAGCACAAAATTTTAGTGGTAtatt                                                                 >  1:8162743/1‑150 (MQ=16)
                                                                                    gaccAAAATGTTTTCTTACGGTGGGTTCGCCCCCTCGCCGGGGGGGTGCGGGCGGGGGGTCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                <  1:3016227/150‑1 (MQ=255)
                                                                                       gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAACCCACACCCACACCGGAGAGGGGACCCAGCGACCACCAAAATTTAGAGATAtcttattc                                                             >  2:4516856/1‑150 (MQ=21)
                                                                                          aatttATTCGTTTGGTTGTTTTTCGCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                          <  1:1637537/146‑1 (MQ=255)
                                                                                               tttCTTACGGTGGGTTTACTCCCTCTCCCGGTGGGTGTGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  2:9405387/149‑1 (MQ=255)
                                                                                                  cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                  <  2:7111709/150‑1 (MQ=255)
                                                                                                   ttttggttgtttTTCTCCCTCGCCGGGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  2:8384/146‑1 (MQ=255)
                                                                                                   gtttggtgggttTGCTCCCTCGCCGGGGGGGTGGGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  1:9501670/146‑1 (MQ=255)
                                                                                                      aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  2:5941244/150‑1 (MQ=255)
                                                                                                         tgggttTTCTCCCTCTCCCTGTGGGCGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:3423140/148‑1 (MQ=255)
                                                                                                          tggttTTCTCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:6135264/150‑1 (MQ=255)
                                                                                                          gggttGTCTCCCTCCCCCTGTGGGTGCGGGCCGGGGTGAGCGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:8967225/149‑1 (MQ=255)
                                                                                                                tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                    <  2:1355471/150‑1 (MQ=255)
                                                                                                                  gccctctcggtGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                  <  2:7645810/149‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCGAACCCGATCCCGCCAGGGAGGGGGGGCCGGTCGGGCCCCAATTTTGTGAGGTATCTTCTTCAGTCTTGTTACAGGAGCACAGAAAATCg                                 >  2:2696207/1‑150 (MQ=21)
                                                                                                                    cctcgccctGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:10038001/150‑1 (MQ=255)
                                                                                                                         ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAACCACTCCCGCAAGGGGGAGGGGGCCGGTCGGGCACAAAATTTGTGGGATATATTATTATGTCTTTTAACAGAAGAAAAGAAAATCGGGCgtt                           >  2:7262620/1‑150 (MQ=19)
                                                                                                                          cctGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                          <  1:10720002/150‑1 (MQ=255)
                                                                                                                              tgggtgggGGCGGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  2:6835720/150‑1 (MQ=255)
                                                                                                                                ggcgggggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:4932802/145‑1 (MQ=255)
                                                                                                                                 gagtGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                   <  2:284490/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAACCCACTCCCACAAGGGAGAGGGGGCCGAACGAGCACAAATTTTGTGAGATATCATCTTCTGTCGTGTAACAGAAGAACCGAAAAAAGGGCTTTTGCGCCGGa                  >  2:10987088/1‑150 (MQ=255)
                                                                                                                                   gagGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTTCGCCTGaa                 <  2:2641772/150‑1 (MQ=255)
                                                                                                                                     gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata               <  1:11131818/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACACCCACAAGGGGGAGGGGCCCCACAGAGCACAAAAGTTTGGAGAAATAATATTCTTGGGTGTAAAAGAAGTACAGAAAATCGGGCGTTTTCGCCAGtataa              >  2:2578007/1‑149 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCCCACCCGAACCCGCTCCCGCCAGGGCGAGCGGCCCCCTCCCGCCCCAATGGTGAGAGATAAATTCTTCTGGCGGGGAACAGCAGAACCGCACAGCGGGCGGTTGCGCCGGacgag              >  2:10122746/1‑146 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATGTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  2:8997364/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:7653905/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:7535004/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:10113578/150‑1 (MQ=255)
                                                                                                                                          ggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGGAGAACTGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:5110124/150‑1 (MQ=255)
                                                                                                                                          ggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:11005171/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACGCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:7557684/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:5923829/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:8331718/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:3329714/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:4676592/150‑1 (MQ=255)
                                                                                                                                            ggggcgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:1070607/145‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:8106545/150‑1 (MQ=255)
                                                                                                                                              tGAGGGCCTCAGCGCGCACGTTCACCCACACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga      <  1:5645625/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCACACCCAAAACCACAACCACAAGGGAGAGGGGACCCATCGAGCACAAAATTTGAGAGAAATAATATTAATTATTGTAACAGAAGAAAACAAAAACGGGCTTTTTCGCCAGAATATTTCGGAGaaaaa  >  2:8279411/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTG‑CACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164433‑164727

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑