Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 10.1% 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.1% 418.1 / 4.0 219intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (47/117);  new base C (19/6);  total (78/141)
Fisher's exact test for biased strand distribution p-value = 8.32e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326
                                                                                                                                                      |                                                                                                                                                     
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCa                                                                                                                                                       >  2:2218244/1‑150 (MQ=255)
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCa                                                                                                                                                       >  2:6954190/1‑150 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAAGCACTGGCGCAAGGAGCGC‑GGAggg                                                                                                                                                    >  1:5715377/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTACAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCGGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑gggggggg                                                                                                                                                  >  1:4461039/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAggggg                                                                                                                                                  >  2:49172/1‑150 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCgggggg                                                                                                                                                 >  2:6385265/1‑149 (MQ=255)
       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGggcg                                                                                                                                                >  2:1189113/1‑150 (MQ=255)
       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcg                                                                                                                                                <  2:9191416/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGggcgg                                                                                                                                               >  2:2964690/1‑150 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcgg                                                                                                                                               >  2:7565674/1‑150 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCAGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGggcgg                                                                                                                                               <  2:6973025/150‑1 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCa                                                                                                                                            >  2:9142298/1‑150 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGccca                                                                                                                                           >  2:7799265/1‑148 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATc                                                                                                                                         <  1:6236348/150‑1 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATc                                                                                                                                         <  2:2617115/150‑1 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCgggggggggccaaac                                                                                                                                        >  2:3504774/1‑143 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGGGCGGccccacc                                                                                                                                        >  1:423007/1‑145 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcg                                                                                                                                        <  2:819216/150‑1 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcg                                                                                                                                        <  2:8537832/150‑1 (MQ=255)
                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgcc                                                                                                                                      <  2:1681286/150‑1 (MQ=255)
                   aTTTAATATTATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATACTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATcgccgc                                                                                                                                    <  1:4787420/150‑1 (MQ=255)
                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCGATTGgccgccccc                                                                                                                                <  2:1049044/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:1189932/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:8983279/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGcctcct                                                                                                                            <  2:5609155/150‑1 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGcccccac                                                                                                                           >  2:6651451/1‑148 (MQ=255)
                               atTTTAAAATTAATTAATGTCATCTGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTGGTATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAc                                                                                                                        <  2:5316240/150‑1 (MQ=255)
                               atTTTAAAATTAATTAATGTCATCAGGTCCTAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAc                                                                                                                        <  1:2849249/150‑1 (MQ=255)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCCTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCCCCGCCCCCCCCAGCAccc                                                                                                                      >  2:5030146/1‑150 (MQ=255)
                                  ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcccc                                                                                                                     <  2:6717044/150‑1 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCACGGGGCGGCCAAACGGCGCCCCCCCCAGCAccccc                                                                                                                    >  2:5836028/1‑150 (MQ=38)
                                    aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCg                                                                                                                   <  2:4738933/150‑1 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGACAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                 <  2:706052/150‑1 (MQ=255)
                                        ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGCCCACCCCCCCCCCCCCCCCCAcccccccccc                                                                                                               >  1:2004584/1‑145 (MQ=255)
                                           attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGGGCGGCCCAAACCCCCCCCCCCCAGCACCCCCCGGCGCTg                                                                                                            >  2:9026398/1‑150 (MQ=255)
                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCG‑CCCCTGCACCCCCGGGCTCTGGCg                                                                                                         >  1:6598038/1‑150 (MQ=255)
                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGGCGGGCCAACGCCCCCCCCCCCAGCCCCCCCCGCCGCAGGCGAACa                                                                                                     >  1:5380755/1‑150 (MQ=37)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACaaaa                                                                                                  <  2:2740055/150‑1 (MQ=255)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCGTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTGGCGAACaaaa                                                                                                  >  1:8233853/1‑150 (MQ=255)
                                                         gTCCCAAAATAACTATAATATTTCATTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGATTCATATAATCACTGGCGCATGGTGCGC‑GCCGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                              <  1:5556223/150‑1 (MQ=21)
                                                           ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGCCCCAGCCCGCCCCCCCCACCACCCCCCGGCCCGGGCGCCCaaaacacccg                                                                                            >  1:4564146/1‑144 (MQ=255)
                                                                 aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcgc                                                                                      <  2:2811167/150‑1 (MQ=255)
                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCGCACAAAAACGCCGCGGCGCCGTGc                                                                                 >  2:1438658/1‑150 (MQ=21)
                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCGCCCCCTGCAACCCCCGGCTCTGGCCAACAAAATCGCCCCTGCGCGCTGc                                                                                 >  2:6058626/1‑150 (MQ=21)
                                                                          atatTTCAGTCTCTCATCCTGTGGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCACGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:4107633/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:444683/150‑1 (MQ=255)
                                                                              ttCTGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCTTTGCTTAATATAAACCCTGGCGCACGGAGGGC‑GCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCt                                                                         <  1:387127/150‑1 (MQ=17)
                                                                                  gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGTGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  2:889277/150‑1 (MQ=255)
                                                                                  gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  2:8379792/150‑1 (MQ=255)
                                                                                   tctctcATCCTGTTGCGCTCCTGTCTGGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  1:147867/150‑1 (MQ=17)
                                                                                   tctctaTTCGTGTTGCGCTCCTGTCTTGTGCATTGCTTAATATAATCAGTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  1:7282421/150‑1 (MQ=17)
                                                                                      ctcATCCTGTTGCGCTCCTGTCTTGTGTATTGCTTCATATAATCACGGGCGCACGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                 <  1:1096264/150‑1 (MQ=18)
                                                                                       tcttcttGTTGGGCTCGTGTATTGTGATTTGCTTCATATAATCCTGGGCGCATGGTGGGC‑GCTGGGGGCGGCAATTCGCCGCGGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTTTCCCTTa                                                                <  1:8423685/144‑1 (MQ=14)
                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:3102246/150‑1 (MQ=255)
                                                                                        cATCGTGTTGCGCTCGTGTCTTGTGATTTGCTTCATATAATCACGGGCGCATGGAGCGC‑GCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                               <  1:6385265/150‑1 (MQ=17)
                                                                                         aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCGGCGCGGTGCCCGCGGCGTATCCCGTACg                                                              >  2:6204183/1‑150 (MQ=18)
                                                                                             tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCCCCGCCCCCAGCACCCCCCGGCGATGGCGAACAAAATAGCCGCGGCGCGGTGCCCACGGCATATACCATaacgcga                                                          >  2:2251449/1‑144 (MQ=17)
                                                                                              gTTGCGCTCCTGTCTTGTGCTTTGCTACATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAAGCGCGGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  1:4159410/150‑1 (MQ=14)
                                                                                               ttGCGCTCCTGTCATGTGCTTTGATTCATATAATCACGGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  1:8927110/150‑1 (MQ=255)
                                                                                                 gcgcTCCTGTCTTGTGCATTGCTTCATATAAACACGGGCGCACGGAGCGC‑GCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  1:2577450/150‑1 (MQ=17)
                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  2:6613477/150‑1 (MQ=255)
                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  2:6874286/150‑1 (MQ=255)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGGCAAAAGCCCCCGCCCCCAGCACCCCCCGGCTCTGGCGAACAAAAAAGCCGCGGCGCGGTGCCCGCGGCTTATACCATACGGCGACCGGGTCggg                                                >  2:8347851/1‑150 (MQ=17)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                <  2:7669303/150‑1 (MQ=255)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCCCCCCCCCCCGCACCCCCCCGGGCGGGCGAACAAAAAAGCCCCGGCGCGGTGGCCCCCGCGTATACCATAACGCAACCCGGTCggg                                                >  2:3785211/1‑150 (MQ=12)
                                                                                                            aTGTGCATTGCTTCATATAATCCCTGGCGCATGGAGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:8330404/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCATGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:7074512/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:7446804/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:1884955/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGCGCGC‑GCCGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:34561/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:3182691/150‑1 (MQ=255)
                                                                                                               tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCGCCCCCTGCACCCCCCGGCTCATGCGAACAAAATCCCCGCAGCGCGGTGCCCTCGGCGTATCCCATACGGCTACCCGGTCGGGCGCGAGGt                                        >  2:5227795/1‑150 (MQ=14)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCAAAACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:3142513/150‑1 (MQ=255)
                                                                                                                 cATTGCTTAATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  1:5952317/150‑1 (MQ=255)
                                                                                                                 cATTGATTCATATAATCACTGGCGAACGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCTAACAAAATCGCCTCTGCGCGGTGCCCTCTGCTTTTCCCTTACGGCTACCGGGTCGGGCGCGATGTaa                                      <  1:5033963/150‑1 (MQ=14)
                                                                                                                  aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAAc                                     <  2:9113473/150‑1 (MQ=255)
                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACa                                    <  2:5757286/150‑1 (MQ=255)
                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACa                                    <  1:9081233/150‑1 (MQ=255)
                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACAGCTACCGGGTCGGGCGCGAGGTAACa                                    <  2:4303639/150‑1 (MQ=255)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGGGGGGCAAAAGCCGCCGCCCCCAGCACCCCCCGGCGCTGGGGTACAAAAAAGCCGCTGGGCGGTGGCCGAGGGGTTTACCATAAGGCTACCGGGTGGGGGGCGGGGttaaatga                                 >  2:5146256/1‑143 (MQ=9)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCCAAACCCCCCCCCCCCAGCCCCCCCCCGCGGGGGCCAACAAAAAAGCCCCGGCCCGGGGGCCACCGCCGATACCATAAAGCGACCCGGGCGGGCGCGGGGGAAccacc                                 >  2:6494511/1‑146 (MQ=255)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:5898032/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:4820841/150‑1 (MQ=39)
                                                                                                                         cATATAATCACTGGCGCAAGGAGCGC‑GCACGGGGCGGCCAAAAGCCCCCCCCCCCAGCACCCCCCGGCGCTTGCGAACAAAAAACCCCCGGCGCGGTGCCCGCCGCGTAAACCATACGGCGACCCGGTCGGGCGCCGGGTTACAAcccag                              >  2:1363176/1‑148 (MQ=2)
                                                                                                                              aaTCACTGGCGCAAGGAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:6948043/150‑1 (MQ=21)
                                                                                                                              aaTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCTTGCACCCCCGGGCTCTGGCGAACAAAAGCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGCGGGAACCTCCCTGTacaa                         >  2:4191062/1‑150 (MQ=11)
                                                                                                                               aTCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:7351287/150‑1 (MQ=18)
                                                                                                                               aTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  2:8319422/150‑1 (MQ=33)
                                                                                                                                tCACTGGCGCAAGGAGCGC‑GCACGGGGCCGCCCAACGCCGCCCCCCCCAGCCCCCCCCGGCGCGTGCGAACAAAATCGCCCCGGCGCCGGGCCCACCGCATAACCCATACCGCGACCCGGGAGGGCGCCGGGTAACAAACCTGAaaaaag                       >  2:6018580/1‑148 (MQ=255)
                                                                                                                                 cACGGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:423007/150‑1 (MQ=17)
                                                                                                                                   cTGGCGCAAGGTGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:2401432/150‑1 (MQ=17)
                                                                                                                                   cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:7191712/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:1782847/150‑1 (MQ=32)
                                                                                                                                     ggCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  2:8722022/150‑1 (MQ=32)
                                                                                                                                         cAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:2064745/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  1:9098820/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  1:4489609/150‑1 (MQ=32)
                                                                                                                                           aGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:5313402/150‑1 (MQ=32)
                                                                                                                                             gAGCGC‑GCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  1:49172/150‑1 (MQ=17)
                                                                                                                                               gcgc‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca        <  2:6770645/150‑1 (MQ=32)
                                                                                                                                                      |                                                                                                                                                     
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326

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