Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 19.7% intergenic (+112/+261) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,1770AC19.7% 154.1 / 7.8 114intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (26/50);  new base C (18/2);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 7.24e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326
                                                                                                                                                     |                                                                                                                                                      
aCTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTACTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAg                                                                                                                                                       <  1:3965047/150‑1 (MQ=255)
  tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggg                                                                                                                                                     >  2:2241644/1‑150 (MQ=255)
    gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggggg                                                                                                                                                   <  2:1900584/150‑1 (MQ=255)
     tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGACTATTTCAGTCTCTCATACAGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGggc                                                                                                                                                  >  1:1691547/1‑150 (MQ=255)
     tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggc                                                                                                                                                  <  2:2669978/150‑1 (MQ=255)
     tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggc                                                                                                                                                  <  2:136935/150‑1 (MQ=255)
      ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcg                                                                                                                                                 >  1:2989547/1‑150 (MQ=255)
       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGC‑GCGGGGggcgg                                                                                                                                                >  1:270951/1‑150 (MQ=255)
       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGC‑GCCGGGggcgg                                                                                                                                                >  1:1405535/1‑150 (MQ=255)
        tgACTGTCGTATTTTAGATGATATTTTAAAATTAATTAATTTAATGAGGTCCTTATCTTACGAGAATATTTCTGTCTGTCATCCTGTTGCGCTCGTGTCTTGTGCATTGCTTAATATAATCACTGGCTGAATGATCGC‑GCAggggtattc                                                                                                                                               <  2:4184098/150‑6 (MQ=255)
        tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcggc                                                                                                                                               <  1:1414255/150‑1 (MQ=255)
        tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcggc                                                                                                                                               <  2:4191081/150‑1 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCa                                                                                                                                             <  1:1046504/150‑1 (MQ=255)
           cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCggggggggggca                                                                                                                                            >  1:3280854/1‑146 (MQ=255)
           cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCaa                                                                                                                                            <  1:3027595/150‑1 (MQ=255)
           cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGccc                                                                                                                                             >  1:3280794/1‑149 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGCGGGGGGGcccca                                                                                                                                           >  2:368618/1‑147 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCgggggggggcgca                                                                                                                                           >  1:2747745/1‑145 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGGGGCCaaa                                                                                                                                           >  1:3893777/1‑149 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGcccag                                                                                                                                           >  2:2777426/1‑147 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCCGCCaaaac                                                                                                                                         >  2:2073024/1‑147 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGGGCGGcccaacc                                                                                                                                         >  2:4021262/1‑145 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAtag                                                                                                                                         >  2:310242/1‑148 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCACGGGgccgcccaaac                                                                                                                                         >  2:4021093/1‑145 (MQ=255)
                gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgcc                                                                                                                                       <  2:445492/150‑1 (MQ=255)
                gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgcc                                                                                                                                       <  1:3310360/150‑1 (MQ=255)
                  aTTTAAGATGATATTTTAAAATTAATTAATGTCATCGGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgccgc                                                                                                                                     <  1:122586/150‑1 (MQ=255)
                   tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCGATTggccgcc                                                                                                                                    <  1:2002086/150‑1 (MQ=255)
                      aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGCGGGCCAACgcccccccc                                                                                                                                 >  1:3786887/1‑150 (MQ=255)
                      aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGGGGGGCAAAAgcccccgcc                                                                                                                                 >  2:570778/1‑150 (MQ=255)
                      aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGGCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGGCCAGCgcccccccc                                                                                                                                 >  1:1107705/1‑150 (MQ=255)
                       aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                                <  1:3693046/150‑1 (MQ=255)
                        tatgatATTTTAAAATTAATTAATTTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCGTGTCATGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  1:3270625/149‑1 (MQ=255)
                        gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  2:2701651/150‑1 (MQ=255)
                        gatgatATTTTAAAATTAATTAAGGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCCGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  1:4165814/150‑1 (MQ=255)
                            atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCAAAACCCGCCGCCCCCAGc                                                                                                                           >  1:321795/1‑150 (MQ=38)
                            atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCGCCGCCGCCCCCGGc                                                                                                                           >  1:321892/1‑150 (MQ=255)
                              tttttAAAATTAATTATGTTCATCAGGTCCTAAAATAACGAGAATATTTCAGTCTCTCTTCCTGTTGCGCTCCTTTCTTGTGCTTTGCTTCATATAATCACTGGCGCAATGAGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCAc                                                                                                                         <  1:1172848/149‑1 (MQ=25)
                               ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCCAccccccccccccccccccc                                                                                                                         >  1:1939388/1‑145 (MQ=255)
                                tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCCAACGCCCCCCCCCCCCCCAccc                                                                                                                       >  2:1563970/1‑150 (MQ=38)
                                  tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGGGGGCCAAAGCCGCCGCCCCCCGCAccccc                                                                                                                     >  2:2549304/1‑150 (MQ=38)
                                  tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAGACCCCCCCCCCCCCGCAccccc                                                                                                                     >  2:2549331/1‑150 (MQ=38)
                                  tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCGTCATATACTCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCAAAAGCCCCCCCCCCCCCCAccccc                                                                                                                     >  2:2549329/1‑150 (MQ=255)
                                    aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGCCCAACGCCGCCGCCCCCCGCCCCCCCgg                                                                                                                   >  2:3794012/1‑150 (MQ=38)
                                    aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGCCCAAAGCCCCCCCCCCCAGCAccccccc                                                                                                                   >  2:2648347/1‑148 (MQ=38)
                                     aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                  <  2:4127298/150‑1 (MQ=255)
                                      attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGCCCAACCCCCCCCCCCCCAccccccccccgc                                                                                                                 >  1:1962185/1‑146 (MQ=255)
                                      attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGGCGGCCCAACGCCGCCCCCCCCCGCACCCCCCGGc                                                                                                                 >  1:2472256/1‑150 (MQ=255)
                                      attaataaATGTCATCAGGTCCGAAAATAACGAGAATAATACAGACTATCATCCAGTTGCGCTCATGACATGTGCATTGCTTCATATAATCACTGTCGCAAGGAGCGC‑GGAGGGCGCGGGCCAACACACCCCCCCCCCCCACCCCCcggg                                                                                                                 >  1:1258343/1‑149 (MQ=255)
                                           ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGGCCATACCCCCCGCCCCCCGCACCCCCCGGcgcggg                                                                                                            >  1:163922/1‑145 (MQ=255)
                                           ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGGCCAACGCCCCCCCCCCCCGCACCCCCCGGCTCggg                                                                                                            >  1:3838555/1‑150 (MQ=255)
                                            taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGCCCAACCCCCCCCCCCCCCCCCCCCCCCCGCACTGGc                                                                                                           >  2:2065111/1‑150 (MQ=37)
                                            taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCCAACGCCCCCCCCCCCCGCCCCCCCCGGCGCCGGc                                                                                                           >  2:2065043/1‑150 (MQ=37)
                                            taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGGGCGGCCACACCCCCCCCCCCCCCGCCCCCCCCGGCGCCGGc                                                                                                           >  2:2065028/1‑150 (MQ=37)
                                             aatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAACCCCCCCCCCCCAGCAACCCCCGGCTCTGGCg                                                                                                          >  2:561410/1‑150 (MQ=37)
                                              atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGCGGGGGGGCCCAACGCCGCCCCCCCCACCCCCCCCCCGCCCTGGCG‑a                                                                                                        >  1:3473816/1‑150 (MQ=255)
                                               tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATATTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑aa                                                                                                       <  2:1940216/150‑1 (MQ=255)
                                               tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑aa                                                                                                       <  2:1164053/150‑1 (MQ=255)
                                                  catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAACGCCCCCCCCCCCACCACCCCCCGGCGCTGGCG‑ACCaa                                                                                                    >  2:244501/1‑150 (MQ=37)
                                                      tgttCCTAAAATATCTATAATTTTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGATTTGCTTAATATAATCCGGGGCGCAAGGAGGGC‑GCTGGGGGCGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATc                                                                                                <  1:2211668/147‑1 (MQ=18)
                                                       ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATcg                                                                                               <  2:2288497/150‑1 (MQ=255)
                                                         tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCAGGGGCGGCCCAAACCCCCCCCCCCCCGCACCCCCCCGCCCAGGCGCAACaaaaaccc                                                                                              >  2:2194902/1‑146 (MQ=25)
                                                         tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGGGGCCCAAAGCCGCCCCCCCCAGAACCCCCCGGCACGGGCG‑AACaaaaaagcc                                                                                             >  2:571923/1‑145 (MQ=25)
                                                                    cGATAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCATGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTg                                                                                  <  2:915749/150‑1 (MQ=18)
                                                                    cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTg                                                                                  <  2:3036537/150‑1 (MQ=255)
                                                                     gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCACGGGGCGGCCCAAAGCCCCCCCCCCCAGCCCCCCCCGGCAAAGGCC‑AAAAAAAACCCCCCAGCGCCGTGc                                                                                 >  2:2698271/1‑150 (MQ=21)
                                                                          tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                            <  2:3538902/150‑1 (MQ=255)
                                                                          tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                            <  2:2194365/150‑1 (MQ=255)
                                                                           atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:2989547/150‑1 (MQ=255)
                                                                            tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑ACAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                          <  2:1844113/150‑1 (MQ=255)
                                                                                  tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCACGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:3930935/150‑1 (MQ=255)
                                                                                  tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGAGTCATATAATCACTGGCGCCAGGAGCGC‑GCAGGGGGCGGCCATTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:1084219/150‑1 (MQ=255)
                                                                                      tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGGGGCCCAAACCCCCCGCCCCCCGCACCCCCCGGCGCGGGCG‑GACAAAAACCCCGCGGCGCGGTGGCCGCGGCGTATCCCCTa                                                                >  2:1335426/1‑150 (MQ=18)
                                                                                      tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGGCAAAAGCCGCCCCCCCCAGCACCCCCGGGCTCTGGCG‑AACAAAATAGCCCCGGCCCGGTGGCCGCGGCGTATcccaaa                                                                >  2:1347769/1‑147 (MQ=18)
                                                                                                     ccTGTCATGTGAATTGCTTCATATAAGCACTGGCGCATGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  1:368618/150‑1 (MQ=255)
                                                                                                           aTGTGCTTTGATTCATATAATCACTGGCGCACGGAGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:36376/150‑1 (MQ=255)
                                                                                                              tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCGGGCACTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCAAGGCTTATCCCATACCGCTAACCGGTCGGGCGCGAGGt                                        >  2:3192074/1‑150 (MQ=14)
                                                                                                                      ttCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATccc                                <  1:244501/150‑1 (MQ=255)
                                                                                                                       tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCGCCGCCGCCCCCAGCACCCCCCGGCTCTGGCG‑AACAAAAACCCCGCTGCGCGGGGCCCGCGGCTTATCCCATACGGCACCCCGGTCGGGCGCGCGGTAACAAcccg                               >  2:2469697/1‑149 (MQ=14)
                                                                                                                       tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:786136/150‑1 (MQ=39)
                                                                                                                         atatAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                             <  1:4081494/150‑1 (MQ=38)
                                                                                                                         atatAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCCCCGCCCCCAGCACCCCCCGGCTCAGGCG‑AACAAAATCGCCGCGGCGCGGTGCCCACCGCGTATACCATACGGCTAACCGGTCCGGCGCGAGGTTACAACCCAGt                             >  2:4214936/1‑150 (MQ=11)
                                                                                                                            taATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaa                          <  1:255652/150‑1 (MQ=35)
                                                                                                                             aaTCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:2538518/150‑1 (MQ=21)
                                                                                                                             aaTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:3193163/150‑1 (MQ=34)
                                                                                                                             aaTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCCCCCCCCCCGGCACCCCCCGGCGCTGGCG‑AACAAAAACGCCCCGGCGCCGGGGCCACGGCGTATACCATACCGCTACCCGGTCGGGCGCCAGGTAAAAACCCAGTaaaa                         >  2:2311631/1‑150 (MQ=2)
                                                                                                                                 aCTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  2:3514735/150‑1 (MQ=32)
                                                                                                                                  cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:3845465/150‑1 (MQ=32)
                                                                                                                                    ggCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  2:3085726/150‑1 (MQ=32)
                                                                                                                                        cAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:193942/150‑1 (MQ=32)
                                                                                                                                        cAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:292467/150‑1 (MQ=32)
                                                                                                                                        cAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  2:3955827/150‑1 (MQ=32)
                                                                                                                                           ggAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  2:3074802/150‑1 (MQ=32)
                                                                                                                                                 c‑gcAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  2:413428/150‑1 (MQ=17)
                                                                                                                                                 c‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  1:2819192/150‑1 (MQ=32)
                                                                                                                                                 c‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  1:2073081/150‑1 (MQ=32)
                                                                                                                                                 c‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  2:882577/150‑1 (MQ=32)
                                                                                                                                                    cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  1:3869089/150‑1 (MQ=32)
                                                                                                                                                     aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:66527/1‑150 (MQ=32)
                                                                                                                                                     aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:290355/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                      
ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326

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