Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 10.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,1770AC10.7% 417.2 / 3.3 201intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (35/119);  new base C (20/4);  total (61/140)
Fisher's exact test for biased strand distribution p-value = 1.49e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326
                                                                                                                                                      |                                                                                                                                                      
 aCTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAg                                                                                                                                                       <  2:1079778/150‑1 (MQ=255)
    ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgggg                                                                                                                                                    <  1:6544722/150‑1 (MQ=255)
    ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgggg                                                                                                                                                    <  1:3350466/150‑1 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                   <  1:2832279/150‑1 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                  <  2:1121334/150‑1 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                  <  2:3628855/150‑1 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                  <  2:2652837/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  1:5112003/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               >  1:3276453/1‑150 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcc                                                                                                                                             >  1:376594/1‑150 (MQ=255)
          gACTGGCGTATTTAAGATGATAATTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATAGTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCgggggggggc                                                                                                                                             >  1:2801428/1‑148 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGgggggggggggcc                                                                                                                                            >  1:1239874/1‑147 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCgggggggggcc                                                                                                                                            >  1:3538636/1‑147 (MQ=255)
                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgcc                                                                                                                                      <  1:6343220/150‑1 (MQ=255)
                   aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgc                                                                                                                                    <  2:6141444/150‑1 (MQ=255)
                   aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgc                                                                                                                                    <  1:1245785/150‑1 (MQ=255)
                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCCGCCAATcgccgccgc                                                                                                                                 >  2:4428442/1‑150 (MQ=255)
                      tAAGATGATAGTGTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgccgc                                                                                                                                 <  2:6409214/150‑1 (MQ=255)
                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                               <  1:6005536/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGTGGCCAATCGCCGCCGcccc                                                                                                                              <  2:2908807/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTCATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                            <  2:309599/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGGCCAAACcccccccccccc                                                                                                                             >  2:788533/1‑150 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                            <  1:6187294/150‑1 (MQ=255)
                           tgatATTGGAAAATTACGTAATGTCATCAGGTCCGCCCAGACCGCGCATCTTTCAGTCTCTCATCCTGGTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                            <  2:5348747/150‑1 (MQ=25)
                               atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGATAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAc                                                                                                                        <  1:177151/150‑1 (MQ=255)
                                ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCAAACccccccccccccccccc                                                                                                                        >  2:1332805/1‑145 (MQ=255)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAACCCCCCcgccccccgccccc                                                                                                                      >  2:4167006/1‑150 (MQ=38)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAAAGCCCCCCCCCCCAGCAccc                                                                                                                      >  1:3895366/1‑150 (MQ=38)
                                  ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCCCCCCCCCCCGCAcccc                                                                                                                     >  2:2626207/1‑150 (MQ=38)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCCATCGGCGCCGCCCCCGGCAccccc                                                                                                                    >  2:2130970/1‑150 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGTCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                    <  1:3623246/150‑1 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                    <  2:3661773/150‑1 (MQ=255)
                                    aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCg                                                                                                                   <  2:1423512/150‑1 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGACTATTTCAGGCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATACTCACTGGCGCAAGGAGCGCGCCCGGGGGGGCCCACccccccccccccccccccccccccc                                                                                                                 >  2:4758/1‑147 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAAACCCCCCGCCCCCCGCccccccccg                                                                                                                 >  2:2578487/1‑147 (MQ=38)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCATGCACCCCCggg                                                                                                                 >  1:4802541/1‑150 (MQ=255)
                                       attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAACGCCCCCCCCCCCACCCCCCCCCGGc                                                                                                                >  2:4450916/1‑150 (MQ=255)
                                        ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGCCCAGACCCGCCCCCCCCAGCACCCCCCGGcg                                                                                                               >  1:5256452/1‑149 (MQ=255)
                                        ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCAACCCCCCCCCCCCCACCACCCCCCGGcg                                                                                                               >  2:4582465/1‑149 (MQ=255)
                                        ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAAGGGGCGGCCCAAACCCCCCCCCCCCACCACCCCCCGGcc                                                                                                               >  2:4582326/1‑149 (MQ=255)
                                           attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                            <  2:5250278/150‑1 (MQ=255)
                                            ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCCACCCCCCCCCCCCCCGCCCCCCCCGGcgccgg                                                                                                           >  2:5281505/1‑145 (MQ=37)
                                            ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCCATCGCCGCCCCCCCCAGCACCCCCCGGCTCTgg                                                                                                           >  2:4272427/1‑150 (MQ=255)
                                            ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGAGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCACTgt                                                                                                           >  2:1097932/1‑149 (MQ=255)
                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                          <  2:2228584/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCGGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCCCTGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:719466/150‑1 (MQ=255)
                                                    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGGCCAAGCCCCCGCCCCCCGCCCCCCCCCGCGCTGGCggccaca                                                                                                   >  1:1926952/1‑144 (MQ=255)
                                                    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGGCACTCGCCGCCCCCCCCACCCCCCCCCGGCGCGGGCGCACaaa                                                                                                   >  2:5869576/1‑150 (MQ=255)
                                                    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAAACGCCCCCCCCCCCACCACCCCCCGGCGCAGGCGCACaaa                                                                                                   >  2:4031141/1‑150 (MQ=25)
                                                    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACaaa                                                                                                   <  2:2631662/150‑1 (MQ=255)
                                                      caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCCAGAGCCCCCCCCCCCCCCACCCCCCGGCGCAGGCGGACaaaaa                                                                                                 >  1:3817180/1‑149 (MQ=25)
                                                       aGGTCCCAACCTCACGAGAAGATTTCACTCTCTCCTCCGGTTGCGCTCCAGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                                <  2:4280418/150‑1 (MQ=18)
                                                          tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                             <  2:5333306/150‑1 (MQ=255)
                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCATGGGGCGGCCAATCGCCGCCCCCCCCAGCACCCCCCGGCGAATGCCAACAAAAACGCCCCAGCGCCGTGc                                                                                 >  2:5984594/1‑150 (MQ=21)
                                                                       ttaaTATTTCAGTCTCTCATCCTGTTGCGCTCCTTTCTTGTGATTTGCTTAATATAATCACTGGCGCATGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGcc                                                                                <  1:4340600/148‑1 (MQ=18)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCGGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:5372767/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:3860091/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCGGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:3509144/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGGTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:4827460/150‑1 (MQ=255)
                                                                             tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                          <  2:5095200/150‑1 (MQ=255)
                                                                                  gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  2:6246440/150‑1 (MQ=255)
                                                                                       tcATCCTTTTGCGCTCCTGTCTTGTGCTTTGATTCATATAATCACTGGCGCATGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:1620200/150‑1 (MQ=17)
                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:3869388/150‑1 (MQ=255)
                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:1239960/150‑1 (MQ=255)
                                                                                         acccggTTCCGCACCGGTCAGCTGCATTGCTTACTATACTCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:2930276/145‑1 (MQ=14)
                                                                                         aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:1470142/150‑1 (MQ=255)
                                                                                           ccTGTTGCGCTCCTGTCATTTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGc                                                            <  1:392858/150‑1 (MQ=255)
                                                                                             tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCACTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                          <  2:2121403/150‑1 (MQ=255)
                                                                                              gTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  2:6401934/150‑1 (MQ=255)
                                                                                              gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  2:2683821/150‑1 (MQ=255)
                                                                                               ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  2:18518/150‑1 (MQ=255)
                                                                                               ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  2:4904450/150‑1 (MQ=255)
                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  2:4877398/150‑1 (MQ=255)
                                                                                                   gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCAGCACCCCCCGGCGCTGGCGGACAAAATACCCGCAGCCCCGTGCCCACGGCGTATCCCAAACGGCAACCCGGt                                                    >  2:4327591/1‑150 (MQ=17)
                                                                                                     tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAACCCCCCCCCCCCAGCACCCCCCGGCTAAGGCCAACAAAAACCCCCCTGCCCGGGGGCCACCGCGTATACCATAAGGCTACCCGGTCg                                                  >  2:2753497/1‑150 (MQ=12)
                                                                                                      ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  1:5858331/150‑1 (MQ=255)
                                                                                                         gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                              <  2:376594/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:5435644/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCATGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:4255985/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:5682246/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:209194/150‑1 (MQ=255)
                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACa                                    <  2:5445514/150‑1 (MQ=255)
                                                                                                                      cTTCATATAATCACTGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:6434466/150‑1 (MQ=255)
                                                                                                                             taATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaa                          <  2:6570718/150‑1 (MQ=35)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGGAGGGGGCGGGCCAACGCCCCCCCCCCCAGCACCCCCCGGCGCAGGCGCACAAAAAAGCCCCCGCGCGGGTGCCGCCGCGTATACCATAATGCGACCCCGTCGGGCGCGGGGGAACAAACCAGTaaaaa                        >  2:3566828/1‑149 (MQ=255)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  2:319939/150‑1 (MQ=33)
                                                                                                                                 cACTGGCGCATGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  1:3137962/150‑1 (MQ=17)
                                                                                                                                 cACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  1:2976961/150‑1 (MQ=32)
                                                                                                                                 cACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  1:5357690/150‑1 (MQ=32)
                                                                                                                                  ccTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  1:1135672/149‑1 (MQ=17)
                                                                                                                                   cgggCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:2416497/148‑1 (MQ=17)
                                                                                                                                   cTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:648628/150‑1 (MQ=17)
                                                                                                                                   cTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:1377656/150‑1 (MQ=32)
                                                                                                                                     ggCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  2:6302651/150‑1 (MQ=32)
                                                                                                                                         cAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:2004185/150‑1 (MQ=32)
                                                                                                                                           aGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:4899969/150‑1 (MQ=32)
                                                                                                                                            ggAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  2:5623136/150‑1 (MQ=32)
                                                                                                                                             gAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  2:2484688/150‑1 (MQ=32)
                                                                                                                                                cgcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt       <  2:3560959/150‑1 (MQ=32)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:5145916/1‑150 (MQ=32)
                                                                                                                                                      |                                                                                                                                                      
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑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: