Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,501,771 T→G 6.6% intergenic (+47/+25) mdeA → / ← PP_1309 methionine gamma‑lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,501,7710TG 6.6% 1165.7 / 48.8 353intergenic (+47/+25)mdeA/PP_1309methionine gamma‑lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (188/141);  new base G (11/13);  total (199/154)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 374 total aligned reads displayed.

GCATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTG  >  NC_002947/1501622‑1501920
                                                                                                                                                     |                                                                                                                                                     
 cATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCtt                                                                                                                                                      <  1:5137004/150‑1 (MQ=255)
  aTCACGGGATATCGGAGGGGTTGGTGCGTCTGTCTGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTg                                                                                                                                                     >  2:6766260/1‑150 (MQ=255)
   tCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGc                                                                                                                                                    <  2:2678547/150‑1 (MQ=255)
   tCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGc                                                                                                                                                    >  2:3048463/1‑150 (MQ=255)
    cACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCg                                                                                                                                                   >  2:3739188/1‑150 (MQ=255)
       gggATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAgg                                                                                                                                                >  2:2380849/1‑150 (MQ=255)
        ggATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTCAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGc                                                                                                                                               <  1:8605598/150‑1 (MQ=255)
            atCGGTGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCGCAAGAGGCCttt                                                                                                                                           <  1:1610658/150‑1 (MQ=255)
            atCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCttt                                                                                                                                           <  1:4722649/150‑1 (MQ=255)
               ggAGGGGTTGGTGTGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCGCAAGAGGCCtttttt                                                                                                                                        <  1:6855918/150‑4 (MQ=255)
                gAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTca                                                                                                                                       >  1:6637352/1‑150 (MQ=255)
                 aGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCGCAAGAGGCCttttttgt                                                                                                                                      <  2:7429101/150‑6 (MQ=255)
                     gTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTc                                                                                                                                  >  2:4624286/1‑150 (MQ=255)
                        ggTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCAGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGa                                                                                                                               >  1:4470876/1‑150 (MQ=255)
                         gTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGAt                                                                                                                              <  1:4629560/150‑1 (MQ=255)
                          tGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATc                                                                                                                             >  2:3638875/1‑150 (MQ=255)
                            cGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCaa                                                                                                                           >  1:2429062/1‑150 (MQ=255)
                                     gTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGc                                                                                                                  >  1:2142878/1‑150 (MQ=255)
                                     gTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATGGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGc                                                                                                                  <  2:8010065/150‑1 (MQ=255)
                                        gggCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGGTGCCCGCGAAGAGGCCCTTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCg                                                                                                               <  2:7941998/150‑1 (MQ=255)
                                            tGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGc                                                                                                           >  2:2595148/1‑150 (MQ=255)
                                                   tgtgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcgtc                                                                                                    <  1:1162775/150‑1 (MQ=255)
                                                    gtgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcgtcg                                                                                                   >  2:3318624/1‑150 (MQ=255)
                                                    gtgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTAGAGCCCGTAGCTTtcgtcg                                                                                                   >  2:7760068/1‑150 (MQ=255)
                                                     tgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                                  <  1:3054683/150‑1 (MQ=255)
                                                     tgGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGa                                                                                                  <  1:3854650/150‑1 (MQ=255)
                                                      ggaggaCCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAg                                                                                                 >  2:5450356/1‑150 (MQ=255)
                                                             ctgctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGccc                                                                                          >  2:6264181/1‑150 (MQ=255)
                                                             ctgctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGccc                                                                                          >  1:1307808/1‑150 (MQ=255)
                                                                ctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGc                                                                                       >  2:2068973/1‑150 (MQ=255)
                                                                ctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGc                                                                                       <  1:3235511/150‑1 (MQ=255)
                                                                ctgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGGTGCCCGCGAAGAGGCCAGTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACAGTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGc                                                                                       <  2:3521866/150‑1 (MQ=255)
                                                                 tgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCg                                                                                      <  1:8217133/150‑1 (MQ=255)
                                                                 tgGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCg                                                                                      >  1:7711918/1‑150 (MQ=255)
                                                                  ggCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGa                                                                                     <  2:5174165/150‑1 (MQ=255)
                                                                     aGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCCACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGtt                                                                                  <  1:680648/150‑1 (MQ=255)
                                                                       aTGTCGAACAGGCATTGCAGTCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAg                                                                                <  2:424106/150‑1 (MQ=255)
                                                                        tGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAgg                                                                               <  2:1496192/150‑1 (MQ=255)
                                                                        tGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAgg                                                                               >  1:5906161/1‑150 (MQ=255)
                                                                         gTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgc                                                                              >  2:1077218/1‑150 (MQ=255)
                                                                         gTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGgc                                                                              >  1:940385/1‑150 (MQ=255)
                                                                            gAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCt                                                                           >  2:8189764/1‑150 (MQ=255)
                                                                             aaCAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCtt                                                                          <  2:8362574/150‑1 (MQ=255)
                                                                             aaCAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCtt                                                                          >  1:6314571/1‑150 (MQ=255)
                                                                                aGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCg                                                                       >  2:1243262/1‑150 (MQ=255)
                                                                                 ggCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCgg                                                                      >  2:2036554/1‑150 (MQ=255)
                                                                                  gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGt                                                                     >  1:2633550/1‑150 (MQ=255)
                                                                                  gCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGAAGCTTTCGTCGAGCATGCCAGGCGAGTTAGGCGCTTTCGGt                                                                     >  1:108397/1‑150 (MQ=255)
                                                                                  gCATTGCAGGCCTGCTGTGAGGTTGCCCGCGAAGAGGCCATTGACTCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGCTAGGCGCTTTCGGt                                                                     <  2:7475285/150‑1 (MQ=255)
                                                                                   cATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTg                                                                    <  1:419655/150‑1 (MQ=255)
                                                                                    aTTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGc                                                                   >  2:1707295/1‑150 (MQ=255)
                                                                                     ttGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCa                                                                  >  2:6626115/1‑150 (MQ=255)
                                                                                     ttGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCa                                                                  >  2:4787770/1‑150 (MQ=255)
                                                                                         aGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGt                                                              >  2:5758122/1‑150 (MQ=255)
                                                                                           gCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcg                                                            >  1:4910981/1‑150 (MQ=255)
                                                                                           gCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcg                                                            >  2:8604352/1‑150 (MQ=255)
                                                                                             ctgctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcg                                                          <  2:7251361/150‑1 (MQ=255)
                                                                                               gctgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTcgcgcg                                                        >  1:3540414/1‑150 (MQ=255)
                                                                                                 tgTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGc                                                      >  1:8403656/1‑150 (MQ=255)
                                                                                                    tAGGTTGCCCGCGAAGAGGCCACTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCGTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACt                                                   >  1:7852337/1‑150 (MQ=255)
                                                                                                         tGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTcggc                                              <  2:7747138/150‑1 (MQ=255)
                                                                                                            ccGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgt                                           >  2:1822912/1‑150 (MQ=255)
                                                                                                              gcgAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgtgt                                         <  2:761647/150‑1 (MQ=255)
                                                                                                              gcgAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGgtgt                                         <  1:6006543/150‑1 (MQ=255)
                                                                                                                 aaGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCtt                                      >  2:2492783/1‑150 (MQ=255)
                                                                                                                 aaGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCtt                                      <  1:5225599/150‑1 (MQ=255)
                                                                                                                   gagTCCAGTGAATCCAGCCCACAAAACAGGCCTCGTGCGAGGCCTTTGTCACATTCAGATCAACTCTTGAGCCCGTAGCTTTCGTCGAGCAGGCCCGGCGAGTTAGGCGCTTTCGGTGCAGAGTCGCGCGGCACTTCGGCGGTGTCTTTa                                    <  2:6539636/150‑1 (MQ=255)
                                                                                                                    agGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAg                                   <  1:2291390/150‑1 (MQ=255)
                                                                                                                      gCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCtt                                      >  1:6414162/1‑145 (MQ=255)
                                                                                                                       ccATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTg                                <  1:3739188/150‑1 (MQ=255)
                                                                                                                        cATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                               >  2:6874196/1‑150 (MQ=255)
                                                                                                                        cATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                               >  2:4063141/1‑150 (MQ=255)
                                                                                                                        cATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgg                               >  2:6015238/1‑150 (MQ=255)
                                                                                                                         aTTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCGCATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTggg                              >  1:5686782/1‑150 (MQ=255)
                                                                                                                         aTTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTggg                              >  1:8183911/1‑150 (MQ=255)
                                                                                                                         aTTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTggg                              >  1:1660299/1‑150 (MQ=255)
                                                                                                                          ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcg                                                                                                      >  2:3483460/1‑77 (MQ=255)
                                                                                                                          ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTtcg                                                                                                      <  1:3483460/77‑1 (MQ=255)
                                                                                                                          ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgggg                             >  2:5570415/1‑150 (MQ=255)
                                                                                                                          ttGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTgggg                             >  1:5435579/1‑150 (MQ=255)
                                                                                                                               tCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAgg                        >  2:2157039/1‑150 (MQ=255)
                                                                                                                                 cAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCg                      >  2:5009512/1‑150 (MQ=255)
                                                                                                                                  aGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCgg                     <  2:8469677/150‑1 (MQ=255)
                                                                                                                                    cccACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTc                   >  2:1140841/1‑150 (MQ=255)
                                                                                                                                     ccACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCa                  >  2:2906806/1‑150 (MQ=255)
                                                                                                                                       acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcg                >  1:5827609/1‑150 (MQ=255)
                                                                                                                                       acaAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTGTAGGTGGGGTCAGGCGGTCAcg                >  2:6370315/1‑150 (MQ=255)
                                                                                                                                        caAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgc               >  1:1982612/1‑150 (MQ=255)
                                                                                                                                        caAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgc               >  1:54303/1‑150 (MQ=255)
                                                                                                                                         aaaaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCAcgcg              >  1:183662/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCgg             <  1:6708339/150‑1 (MQ=255)
                                                                                                                                           caaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCggg            >  1:7237091/2‑150 (MQ=255)
                                                                                                                                           aaaaGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCggg            <  2:3958379/150‑1 (MQ=255)
                                                                                                                                            aaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc           >  1:2806992/1‑150 (MQ=255)
                                                                                                                                            aaaGGCCTCGCAAGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGc           >  1:7333397/1‑150 (MQ=255)
                                                                                                                                                  ctctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCg     <  2:4457016/150‑1 (MQ=255)
                                                                                                                                                  ctctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCg     >  2:285772/1‑150 (MQ=255)
                                                                                                                                                   tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCgg    >  2:5079732/1‑150 (MQ=255)
                                                                                                                                                   tctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCgg    >  2:2939026/1‑150 (MQ=255)
                                                                                                                                                    ctTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGt   <  1:50263/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCATCACGGGATATCGGAGGGGTTGGTGCGTCTGTCGGTGGGGCTGGAGGATGTGGAGGACCTGCTGGCAGATGTCGAACAGGCATTGCAGGCCTGCTGTTAGGTTGCCCGCGAAGAGGCCATTGAATCCAGCCCACAAAAAAGGCCTCTTGCGAGGCCTTTGTCACATTCAGATCAACGCTTGAGCCCGTAGCTTTCGTCGAGCATGCCCGGCGAGTTAGGCGCTTTCGGTGCATAGTCGCGCGGCACTTCGGCGGTGTCTTTAGGTGGGGTCAGGCGGTCACGCGGGCCTTGCGGTG  >  NC_002947/1501622‑1501920

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