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.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326
                                                                                                                                                      |                                                                                                                                                     
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                       >  2:2218244/1‑150 (MQ=255)
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                       >  2:6954190/1‑150 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAAGCACTGGCGCAAGGAGCGCGGAggg                                                                                                                                                    >  1:5715377/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTACAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCGGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCgggggggg                                                                                                                                                  >  1:4461039/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAggggg                                                                                                                                                  >  2:49172/1‑150 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCgggggg                                                                                                                                                 >  2:6385265/1‑149 (MQ=255)
       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGggcg                                                                                                                                                >  2:1189113/1‑150 (MQ=255)
       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcg                                                                                                                                                <  2:9191416/150‑1 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGggcgg                                                                                                                                               >  2:2964690/1‑150 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               >  2:7565674/1‑150 (MQ=255)
        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCAGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcgg                                                                                                                                               <  2:6973025/150‑1 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCa                                                                                                                                            >  2:9142298/1‑150 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGccca                                                                                                                                           >  2:7799265/1‑148 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                         <  1:6236348/150‑1 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATc                                                                                                                                         <  2:2617115/150‑1 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCgggggggggccaaac                                                                                                                                        >  2:3504774/1‑143 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGGGCGGccccacc                                                                                                                                        >  1:423007/1‑145 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcg                                                                                                                                        <  2:819216/150‑1 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcg                                                                                                                                        <  2:8537832/150‑1 (MQ=255)
                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgcc                                                                                                                                      <  2:1681286/150‑1 (MQ=255)
                   aTTTAATATTATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATACTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgc                                                                                                                                    <  1:4787420/150‑1 (MQ=255)
                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCGATTGgccgccccc                                                                                                                                <  2:1049044/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:1189932/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                              <  1:8983279/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcctcct                                                                                                                            <  2:5609155/150‑1 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccccac                                                                                                                           >  2:6651451/1‑148 (MQ=255)
                               atTTTAAAATTAATTAATGTCATCTGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTGGTATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAc                                                                                                                        <  2:5316240/150‑1 (MQ=255)
                               atTTTAAAATTAATTAATGTCATCAGGTCCTAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAc                                                                                                                        <  1:2849249/150‑1 (MQ=255)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCCTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCCCCGCCCCCCCCAGCAccc                                                                                                                      >  2:5030146/1‑150 (MQ=255)
                                  ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcccc                                                                                                                     <  2:6717044/150‑1 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCAAACGGCGCCCCCCCCAGCAccccc                                                                                                                    >  2:5836028/1‑150 (MQ=38)
                                    aaaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCg                                                                                                                   <  2:4738933/150‑1 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGACAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                 <  2:706052/150‑1 (MQ=255)
                                        ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCCACCCCCCCCCCCCCCCCCAcccccccccc                                                                                                               >  1:2004584/1‑145 (MQ=255)
                                           attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGGGCGGCCCAAACCCCCCCCCCCCAGCACCCCCCGGCGCTg                                                                                                            >  2:9026398/1‑150 (MQ=255)
                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCTGCACCCCCGGGCTCTGGCg                                                                                                         >  1:6598038/1‑150 (MQ=255)
                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGCGGGCCAACGCCCCCCCCCCCAGCCCCCCCCGCCGCAGGCGAACa                                                                                                     >  1:5380755/1‑150 (MQ=37)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACaaaa                                                                                                  <  2:2740055/150‑1 (MQ=255)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCGTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTGGCGAACaaaa                                                                                                  >  1:8233853/1‑150 (MQ=255)
                                                         gTCCCAAAATAACTATAATATTTCATTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGATTCATATAATCACTGGCGCATGGTGCGCGCCGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                              <  1:5556223/150‑1 (MQ=21)
                                                           ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGCCCCAGCCCGCCCCCCCCACCACCCCCCGGCCCGGGCGCCCaaaacacccg                                                                                            >  1:4564146/1‑144 (MQ=255)
                                                                 aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcgc                                                                                      <  2:2811167/150‑1 (MQ=255)
                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCGCACAAAAACGCCGCGGCGCCGTGc                                                                                 >  2:1438658/1‑150 (MQ=21)
                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCGCCGCCCCCTGCAACCCCCGGCTCTGGCCAACAAAATCGCCCCTGCGCGCTGc                                                                                 >  2:6058626/1‑150 (MQ=21)
                                                                          atatTTCAGTCTCTCATCCTGTGGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:4107633/150‑1 (MQ=255)
                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:444683/150‑1 (MQ=255)
                                                                              ttCTGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCTTTGCTTAATATAAACCCTGGCGCACGGAGGGCGCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCt                                                                         <  1:387127/150‑1 (MQ=17)
                                                                                  gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGTGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  2:889277/150‑1 (MQ=255)
                                                                                  gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  2:8379792/150‑1 (MQ=255)
                                                                                   tctctcATCCTGTTGCGCTCCTGTCTGGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  1:147867/150‑1 (MQ=17)
                                                                                   tctctaTTCGTGTTGCGCTCCTGTCTTGTGCATTGCTTAATATAATCAGTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  1:7282421/150‑1 (MQ=17)
                                                                                      ctcATCCTGTTGCGCTCCTGTCTTGTGTATTGCTTCATATAATCACGGGCGCACGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                 <  1:1096264/150‑1 (MQ=18)
                                                                                       tcttcttGTTGGGCTCGTGTATTGTGATTTGCTTCATATAATCCTGGGCGCATGGTGGGCGCTGGGGGCGGCAATTCGCCGCGGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTTTCCCTTa                                                                <  1:8423685/144‑1 (MQ=14)
                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:3102246/150‑1 (MQ=255)
                                                                                        cATCGTGTTGCGCTCGTGTCTTGTGATTTGCTTCATATAATCACGGGCGCATGGAGCGCGCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTAc                                                               <  1:6385265/150‑1 (MQ=17)
                                                                                         aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCGGCGCGGTGCCCGCGGCGTATCCCGTACg                                                              >  2:6204183/1‑150 (MQ=18)
                                                                                             tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCCCCGCCCCCAGCACCCCCCGGCGATGGCGAACAAAATAGCCGCGGCGCGGTGCCCACGGCATATACCATaacgcga                                                          >  2:2251449/1‑144 (MQ=17)
                                                                                              gTTGCGCTCCTGTCTTGTGCTTTGCTACATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAAGCGCGGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  1:4159410/150‑1 (MQ=14)
                                                                                               ttGCGCTCCTGTCATGTGCTTTGATTCATATAATCACGGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  1:8927110/150‑1 (MQ=255)
                                                                                                 gcgcTCCTGTCTTGTGCATTGCTTCATATAAACACGGGCGCACGGAGCGCGCTGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  1:2577450/150‑1 (MQ=17)
                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  2:6613477/150‑1 (MQ=255)
                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  2:6874286/150‑1 (MQ=255)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAAAAGCCCCCGCCCCCAGCACCCCCCGGCTCTGGCGAACAAAAAAGCCGCGGCGCGGTGCCCGCGGCTTATACCATACGGCGACCGGGTCggg                                                >  2:8347851/1‑150 (MQ=17)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                <  2:7669303/150‑1 (MQ=255)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCCCCCCCCCCCGCACCCCCCCGGGCGGGCGAACAAAAAAGCCCCGGCGCGGTGGCCCCCGCGTATACCATAACGCAACCCGGTCggg                                                >  2:3785211/1‑150 (MQ=12)
                                                                                                            aTGTGCATTGCTTCATATAATCCCTGGCGCATGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:8330404/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCATGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:7074512/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:7446804/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:1884955/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGCGCGCGCCGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:34561/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:3182691/150‑1 (MQ=255)
                                                                                                               tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCGCCGCCCCCTGCACCCCCCGGCTCATGCGAACAAAATCCCCGCAGCGCGGTGCCCTCGGCGTATCCCATACGGCTACCCGGTCGGGCGCGAGGt                                        >  2:5227795/1‑150 (MQ=14)
                                                                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCAAAACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:3142513/150‑1 (MQ=255)
                                                                                                                 cATTGCTTAATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  1:5952317/150‑1 (MQ=255)
                                                                                                                 cATTGATTCATATAATCACTGGCGAACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCTAACAAAATCGCCTCTGCGCGGTGCCCTCTGCTTTTCCCTTACGGCTACCGGGTCGGGCGCGATGTaa                                      <  1:5033963/150‑1 (MQ=14)
                                                                                                                  aTTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAAc                                     <  2:9113473/150‑1 (MQ=255)
                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACa                                    <  2:5757286/150‑1 (MQ=255)
                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACa                                    <  1:9081233/150‑1 (MQ=255)
                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACAGCTACCGGGTCGGGCGCGAGGTAACa                                    <  2:4303639/150‑1 (MQ=255)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGGGGGCAAAAGCCGCCGCCCCCAGCACCCCCCGGCGCTGGGGTACAAAAAAGCCGCTGGGCGGTGGCCGAGGGGTTTACCATAAGGCTACCGGGTGGGGGGCGGGGttaaatga                                 >  2:5146256/1‑143 (MQ=9)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCAAACCCCCCCCCCCCAGCCCCCCCCCGCGGGGGCCAACAAAAAAGCCCCGGCCCGGGGGCCACCGCCGATACCATAAAGCGACCCGGGCGGGCGCGGGGGAAccacc                                 >  2:6494511/1‑146 (MQ=255)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:5898032/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:4820841/150‑1 (MQ=39)
                                                                                                                         cATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCAAAAGCCCCCCCCCCCAGCACCCCCCGGCGCTTGCGAACAAAAAACCCCCGGCGCGGTGCCCGCCGCGTAAACCATACGGCGACCCGGTCGGGCGCCGGGTTACAAcccag                              >  2:1363176/1‑148 (MQ=2)
                                                                                                                              aaTCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:6948043/150‑1 (MQ=21)
                                                                                                                              aaTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCTTGCACCCCCGGGCTCTGGCGAACAAAAGCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGCGGGAACCTCCCTGTacaa                         >  2:4191062/1‑150 (MQ=11)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:7351287/150‑1 (MQ=18)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  2:8319422/150‑1 (MQ=33)
                                                                                                                                tCACTGGCGCAAGGAGCGCGCACGGGGCCGCCCAACGCCGCCCCCCCCAGCCCCCCCCGGCGCGTGCGAACAAAATCGCCCCGGCGCCGGGCCCACCGCATAACCCATACCGCGACCCGGGAGGGCGCCGGGTAACAAACCTGAaaaaag                       >  2:6018580/1‑148 (MQ=255)
                                                                                                                                 cACGGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:423007/150‑1 (MQ=17)
                                                                                                                                   cTGGCGCAAGGTGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:2401432/150‑1 (MQ=17)
                                                                                                                                   cTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:7191712/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:1782847/150‑1 (MQ=32)
                                                                                                                                     ggCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  2:8722022/150‑1 (MQ=32)
                                                                                                                                         cAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:2064745/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  1:9098820/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  1:4489609/150‑1 (MQ=32)
                                                                                                                                           aGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:5313402/150‑1 (MQ=32)
                                                                                                                                             gAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  1:49172/150‑1 (MQ=17)
                                                                                                                                               gcgcgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca        <  2:6770645/150‑1 (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 < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: