Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 12.8% 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,1770AC12.8% 366.7 / 12.3 174intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (40/95);  new base C (21/4);  total (64/110)
Fisher's exact test for biased strand distribution p-value = 4.38e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326
                                                                                                                                                      |                                                                                                                                                      
 aCTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACTAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAg                                                                                                                                                       <  1:1900238/150‑1 (MQ=255)
 aCTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAg                                                                                                                                                       <  2:276006/150‑1 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggg                                                                                                                                                     <  1:1453170/150‑1 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggg                                                                                                                                                     <  2:1136905/150‑1 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggg                                                                                                                                                     <  2:2884726/150‑1 (MQ=255)
    ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCgggg                                                                                                                                                    >  2:4753655/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCggggg                                                                                                                                                   >  2:1027142/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggggg                                                                                                                                                   >  1:4484999/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAggggg                                                                                                                                                   <  1:3190664/150‑1 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGggc                                                                                                                                                  >  1:4055098/1‑150 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggc                                                                                                                                                  <  1:1488527/150‑1 (MQ=255)
      tggtcACTGGCGTATTTAAGATGATATTATAGAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGGCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCggggggg                                                                                                                                                  >  1:4547573/1‑149 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGggcgg                                                                                                                                                >  1:5968742/1‑150 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcgg                                                                                                                                                <  1:2054920/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcgg                                                                                                                                                <  1:479888/150‑1 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcggc                                                                                                                                               <  2:3868726/150‑1 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcggc                                                                                                                                               <  2:5735918/150‑1 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcggc                                                                                                                                               <  2:2486677/150‑1 (MQ=255)
          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGcc                                                                                                                                              <  2:3727125/150‑1 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGC‑GCCCGGGGGGGccc                                                                                                                                             >  2:4634430/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcgggccaa                                                                                                                                           >  2:1281644/1‑145 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAt                                                                                                                                           <  2:5319933/150‑1 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAcg                                                                                                                                         >  2:3018941/1‑150 (MQ=255)
                   aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgccgc                                                                                                                                     <  1:902142/150‑1 (MQ=255)
                    tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgccgcc                                                                                                                                    <  2:3542479/150‑1 (MQ=255)
                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgccgccgc                                                                                                                                  <  2:2653268/150‑1 (MQ=255)
                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                                <  2:2069877/150‑1 (MQ=255)
                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                                <  1:287762/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCCCCGcccc                                                                                                                               >  1:1002125/1‑150 (MQ=255)
                          atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGccccc                                                                                                                              <  1:4282985/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCCCTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                             <  2:953692/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGTGGCCCATCGCCGCCGCCCCCt                                                                                                                             >  1:1845075/1‑150 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGTGGCCAAACGCCGCCGccccca                                                                                                                             >  1:1844982/1‑149 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTg                                                                                                                            <  2:4066897/150‑1 (MQ=255)
                             atatTTTAAAATTAATTAATGTCATCTGGTCCGAAAATAACGAGAATATTTCTGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCATGGAGGGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGc                                                                                                                           <  1:5891913/150‑1 (MQ=38)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGCGGGCCAAAGCCCCCCCCCCCCGCc‑ccc                                                                                                                      >  1:5255778/1‑150 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGGCAAAAGCCCCCCCCCCCCGCc‑ccccc                                                                                                                    >  2:5200161/1‑150 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCGCCGCCGCCCCCGGCA‑ccccc                                                                                                                    >  2:3928162/1‑150 (MQ=255)
                                       attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCGCCGCCGCCCCCAGCA‑CCCCCCGGc                                                                                                                >  2:4415368/1‑150 (MQ=255)
                                        ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCCCCCCCCCCGGCA‑CCCCCCGGca                                                                                                               >  2:5903814/1‑149 (MQ=38)
                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCACTGGc                                                                                                          >  2:4007537/1‑150 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGaa                                                                                                       <  2:3612396/150‑1 (MQ=255)
                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGGCCAACGCCCCCCCCCCCCCCC‑CCCCCCCGCCCGGGCGGACa                                                                                                     >  2:83440/1‑150 (MQ=255)
                                                     tcaGGTCCGACACTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACaaaa                                                                                                  <  2:5582464/150‑1 (MQ=255)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACaaaa                                                                                                  <  2:262907/150‑1 (MQ=255)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCCCCCCCAGCA‑CCCCCCCGCGCTGGCGAACaaaa                                                                                                  >  2:1521155/1‑150 (MQ=25)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCCGGGCTCTGGCGAACAAAATc                                                                                                <  2:673145/150‑1 (MQ=255)
                                                        ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCCGGGCTCTGGCGAACAAAATc                                                                                                <  2:673005/150‑1 (MQ=255)
                                                          tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGCCCAACGCCCCCCCCCCCAGCA‑CCCCCCGCCGCAGGGGAaaaaaaaagcc                                                                                             >  1:3425282/1‑145 (MQ=255)
                                                           ccGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGTGGGCCGCCCAAAGCCGCCCCCCCCGGCA‑CCCCCCGGCGCTTGGGAACaaaaaaaccc                                                                                            >  2:4210586/1‑144 (MQ=255)
                                                            cgcacaTAACGCGAATATTTCAGACCCTCACCCTGTTGCGCCCCTGCCATGTGCATGGCTTCATACAAGCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCGCCCCCTGCC‑CCCCCGGGCTCTGGCTAACAAAATcgccgc                                                                                           <  2:232622/145‑1 (MQ=18)
                                                            cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAA‑AGCCCCCCCCCCCCCCC‑CCCCCCCGCCCGGGCGCACaaaaaaacccc                                                                                           >  2:5612696/1‑143 (MQ=255)
                                                            cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGGGGGGCACAAGGCCCCCCCCCCAGCC‑CCCCCCCGCACGGGGGAACaaaaaaaaccc                                                                                           >  1:1803877/1‑143 (MQ=255)
                                                                aaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCGATTGgccgccccc                                                                                                                                 <  2:1868271/109‑1 (MQ=37)
                                                                 aTAACGATAATATTTCTGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGATTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcgc                                                                                      <  2:1045701/150‑1 (MQ=255)
                                                                    aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                   <  2:3691708/150‑1 (MQ=255)
                                                                       agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACCCCCCCGCCCCCCGCA‑CCCCCCGGCACAGGCCAACAAAAAACCCCCAGCGCGGTGcc                                                                                >  2:3324287/1‑150 (MQ=21)
                                                                        gaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCGCCCCCGCCCCCGGCA‑CCCCCCGGCTCTGGCGAACAAAAAAGCCGCGGCGCGGTGccc                                                                               >  2:1008466/1‑150 (MQ=21)
                                                                        gaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCA‑CCCCCCGGCTCAGGCGAACAAAATAGCCGCGGCGCGGTGccc                                                                               >  2:1008656/1‑150 (MQ=21)
                                                                        gaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCACGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCCGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGccc                                                                               >  2:3969956/1‑150 (MQ=255)
                                                                          atatTTCAGCCTCCCATCCTGTGGCGCTCCTGTCATGTGCAGGGCTTCATATACGCCCGGGCGCACGGCGCGC‑GCCGGGGGCGGCCAAGCGCCGCCGCCCCCGGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  1:3317508/150‑1 (MQ=17)
                                                                           tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                            <  2:1537803/150‑1 (MQ=255)
                                                                           tatTTCAGTCTCTCATCCTGTTGCGCTCCGGTCATGTGCTTTGCTTCATATAATCCCTGGCGCACGGCGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACACAATCGCCGCTGCGCGGTGCCCTCg                                                                            <  2:5640809/150‑1 (MQ=18)
                                                                                 aGTCTCTCTTCCTGTTGCGCTCCTTTCTTGTGCTTTGCTAAATATAACCACTGGCGAACGGCGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt                                                                      <  1:753419/150‑1 (MQ=14)
                                                                                     tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCACGGGGCGGCCCAACGCCCCCCCCCCCAGCA‑CCCCCCGGCTCAGGCGAACAAAAAACCCGCGGCGCCGGGGCCCCGGCGTATacca                                                                  >  2:3189470/1‑149 (MQ=18)
                                                                                        cATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                               <  2:3359822/150‑1 (MQ=255)
                                                                                         agccTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:5266520/148‑1 (MQ=255)
                                                                                              gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  2:4021149/150‑1 (MQ=255)
                                                                                              gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  2:4250418/150‑1 (MQ=255)
                                                                                                 gcgcgcCGGTCTTGTGCTTTGCTTCTTCTCATCACGGGGGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGGTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCGTCGGCTTATCCCTTACGGGTACCgg                                                      <  2:4547573/150‑1 (MQ=14)
                                                                                                  cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCACGGGGCGGCCCAACGCCGCCGCCCCCCGCA‑CCCCCCCGCGCAGGCGAACAAAATACCCCCAGCCCGGGGGCCTCGGCGTATCCCAAACGGcgaacccg                                                     >  2:3810922/1‑143 (MQ=12)
                                                                                                    cTCCTGTCTTGTGCTTTGCTTCATATAATCACTGGCGCAAGGTGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  1:1947278/150‑1 (MQ=255)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                <  1:1497582/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:1672819/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGCGGCCAAAACCCCCCCCCCCCCGCA‑CCCCCCCGCGCAGGCCAACAAAAAACCCCCAGCGCGGGGCCCGCCGCAGATACCGTAACGCTAACCGGTCCGGcgcgc                                           >  2:4886712/1‑149 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCCCTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:2238259/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:5417794/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGTAGCGC‑GCAGGGGGCGGCCGATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  1:4925044/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:1471098/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCACGGGGCCGCCAAAAGCCGCCCCCCCCAGCC‑ACCCCCGGCTCGGGCCCACAAAAACGCCCCAGCGCGGGGGCCGCCGCGGATACCATACCGCCGCCCGGTCGGGGGCGGGGGaa                                      >  2:5375593/1‑150 (MQ=9)
                                                                                                                 cATTGCTTCATAGAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCACTCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACCAACTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:4863553/150‑1 (MQ=255)
                                                                                                                      cTTCATATAATCACTGGCGCATGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:3516324/150‑1 (MQ=255)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:4415368/150‑1 (MQ=40)
                                                                                                                       ttCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCCAAAGCCCCCGCCCCCAGCA‑CCCCCCGGCGCAGGCGGACAAAAAACCCGCGGCCCCGTGCCCACCGCATATCCCATAACGCAAACCGGTCGGGCGCGCGGGAACaaacc                                >  2:4917409/1‑146 (MQ=255)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:472204/150‑1 (MQ=39)
                                                                                                                          atatAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                             <  1:3788394/150‑1 (MQ=38)
                                                                                                                          atatAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                             <  1:3788027/150‑1 (MQ=38)
                                                                                                                           tataATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTa                            <  2:4720722/150‑1 (MQ=37)
                                                                                                                                tCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                       <  2:601340/150‑1 (MQ=32)
                                                                                                                                 cACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  1:5225787/150‑1 (MQ=17)
                                                                                                                                      gcgcAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                 <  1:5518847/150‑1 (MQ=32)
                                                                                                                                         cAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  2:5145860/150‑1 (MQ=32)
                                                                                                                                           tggAGCGC‑GCTGGGGGCGGCAAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:5242977/149‑1 (MQ=17)
                                                                                                                                                 gc‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc      >  2:5048764/1‑150 (MQ=32)
                                                                                                                                                 gc‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc      >  2:776827/1‑150 (MQ=32)
                                                                                                                                                 gc‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc      >  2:4988362/1‑150 (MQ=32)
                                                                                                                                                     cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  2:6424/150‑1 (MQ=32)
                                                                                                                                                     cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  2:4814034/150‑1 (MQ=32)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGATACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:5409755/1‑150 (MQ=17)
                                                                                                                                                      aGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCGTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTACCATCCCTGTAAAACGCGCCCTCAGCCCCCCTCCAtg  >  2:2661607/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCA‑CCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  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: ‑