Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,956,334 G→T 19.7% intergenic (‑121/+40) oguA ← / ← PP_2585 8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,956,3340GT19.7% 673.5 / 164.6 284intergenic (‑121/+40)oguA/PP_25858‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (117/111);  new base T (32/24);  total (149/135)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 292 total aligned reads displayed.

AGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACA  >  NC_002947/2956185‑2956483
                                                                                                                                                     |                                                                                                                                                     
      tttGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgc                                                                                                                                                 <  2:4082788/150‑1 (MQ=255)
      tttGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgc                                                                                                                                                 <  1:570689/150‑1 (MQ=255)
      tttGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCcgc                                                                                                                                                 <  2:1315423/150‑1 (MQ=255)
         gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                              <  2:5856756/150‑1 (MQ=255)
         gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                              <  1:3889037/150‑1 (MQ=255)
         gCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCaa                                                                                                                                              <  2:1441652/150‑1 (MQ=255)
         gCAGGGTTTGGAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCaa                                                                                                                                              <  2:2175448/150‑1 (MQ=255)
         gCACGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCATTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGaa                                                                                                                                              <  2:200736/150‑1 (MQ=255)
           aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAaga                                                                                                                                            <  1:540128/150‑1 (MQ=255)
           aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAaga                                                                                                                                            <  1:5224262/150‑1 (MQ=255)
           aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAaga                                                                                                                                            <  1:5046860/150‑1 (MQ=255)
           aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAaga                                                                                                                                            <  1:3087901/150‑1 (MQ=255)
           aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                            <  1:1054064/150‑1 (MQ=255)
           aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                            <  1:2119774/150‑1 (MQ=255)
            gggTTTGTAGATGACAGGGTGCAAGCCGGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAagag                                                                                                                                           <  2:6078136/150‑1 (MQ=255)
             ggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAgg                                                                                                                                          <  2:5918860/150‑1 (MQ=255)
                     gATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACg                                                                                                                                  >  1:6544682/1‑150 (MQ=255)
                           aTGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATATGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGAGGCTCGAGCCGCGAAGAGGCCGGCACGGGTTGc                                                                                                                            >  2:1218777/1‑150 (MQ=255)
                           aTGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGc                                                                                                                            >  2:4958505/1‑150 (MQ=255)
                           aTGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGc                                                                                                                            >  1:117725/1‑150 (MQ=255)
                            tGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCa                                                                                                                           <  1:5952712/150‑1 (MQ=255)
                             ggTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAg                                                                                                                          <  2:1551174/150‑1 (MQ=255)
                                gCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCa                                                                                                                       >  1:1780566/1‑150 (MQ=255)
                                  aaGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATg                                                                                                                     >  1:5302989/1‑150 (MQ=255)
                                   aGCCTGCGCAAGACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                    <  2:2446577/150‑1 (MQ=255)
                                        gcgcAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGAt                                                                                                               <  2:6206213/150‑1 (MQ=255)
                                            aaTACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATcaca                                                                                                           <  1:6586910/150‑1 (MQ=255)
                                                         gATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgct                                                                                              <  1:4632998/150‑1 (MQ=255)
                                                              ccGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTca                                                                                         >  1:2445954/1‑150 (MQ=255)
                                                               agATGTCATAGGCGCTCGGTATTGTAAGCAAGATTTTAACTCATTGTATACAGCTCGAAGTGCCCGTGCAGGCCTATTGTCCGCTAGATCACCGAAGATGCCGGAACTGGTAGCATCAATTAAATAACAAAACACCTAGCGGATGCTcac                                                                                        >  1:1980777/2‑150 (MQ=255)
                                                                 aTGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCacac                                                                                      <  1:5165793/150‑1 (MQ=255)
                                                                  tGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACAcc                                                                                     >  1:102733/1‑150 (MQ=255)
                                                                     cATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTc                                                                                  <  1:3301069/150‑1 (MQ=255)
                                                                       tAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGATCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                <  2:2290479/150‑1 (MQ=255)
                                                                       tAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                >  2:1621993/1‑150 (MQ=255)
                                                                       tAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                <  2:742903/150‑1 (MQ=255)
                                                                       tAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                <  1:7429345/150‑1 (MQ=255)
                                                                       tAGGCGCTCGGTATTGTAAGCAAGATTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                <  1:4995415/150‑1 (MQ=255)
                                                                              tCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCg                                                                         >  1:4398692/1‑150 (MQ=255)
                                                                               cGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                                        <  2:4591035/150‑1 (MQ=255)
                                                                                    ttGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAAt                                                                   >  2:1122018/1‑150 (MQ=255)
                                                                                      gTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATcc                                                                 <  1:5241606/150‑1 (MQ=255)
                                                                                          gcaagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAg                                                             >  2:5013/1‑150 (MQ=255)
                                                                                           caagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGc                                                            >  2:2684372/1‑150 (MQ=255)
                                                                                             agcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCgg                                                          >  1:6048416/1‑150 (MQ=255)
                                                                                               cTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGt                                                        <  1:5522356/150‑1 (MQ=255)
                                                                                                tttttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGtt                                                       >  2:1696087/1‑150 (MQ=255)
                                                                                                tttttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGtt                                                       >  2:7331124/1‑150 (MQ=255)
                                                                                                tttttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCGAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCAGGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCCATCCCCAGCGGGtt                                                       >  2:3996968/1‑150 (MQ=255)
                                                                                                     aCTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGt                                                  <  2:3762358/150‑1 (MQ=255)
                                                                                                       tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttt                                                <  2:2583641/150‑1 (MQ=255)
                                                                                                        gATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGtttt                                               <  2:6395085/150‑1 (MQ=255)
                                                                                                          ttGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTc                                             >  1:6079458/1‑150 (MQ=255)
                                                                                                            gTATACCGGTCGATGTGCCCGTGCCGGCATCTTGGCGGCTAGAGCCGCGAAGTGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                           <  1:2725556/150‑1 (MQ=255)
                                                                                                             tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGgc                                                                                               >  2:1805929/1‑97 (MQ=255)
                                                                                                                aCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCGAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCCGCGGGTTAGCGTTTTTCAACGCa                                       >  2:2174230/1‑150 (MQ=255)
                                                                                                                  aGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATc                                     >  2:4342433/1‑150 (MQ=255)
                                                                                                                  aGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATc                                     >  1:3730112/1‑150 (MQ=255)
                                                                                                                     tCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCgg                                  >  2:1086739/1‑150 (MQ=255)
                                                                                                                       gAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgca                                <  2:4723243/150‑1 (MQ=255)
                                                                                                                              ccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGc                         <  2:25393/150‑1 (MQ=255)
                                                                                                                              ccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGc                         <  2:5105195/150‑1 (MQ=255)
                                                                                                                                 tGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTc                      >  2:6726042/1‑150 (MQ=255)
                                                                                                                                   ccGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAg                    >  2:3229872/1‑150 (MQ=255)
                                                                                                                                    cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                   >  2:6333536/1‑150 (MQ=255)
                                                                                                                                    cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                   <  1:7283796/150‑1 (MQ=255)
                                                                                                                                     ggCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg                  <  2:6851979/150‑1 (MQ=255)
                                                                                                                                     ggCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg                  >  1:4438709/1‑150 (MQ=255)
                                                                                                                                     ggCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg                  <  2:5138795/150‑1 (MQ=255)
                                                                                                                                      gCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGt                 >  1:1525850/1‑150 (MQ=255)
                                                                                                                                      gCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGt                 >  2:367550/1‑150 (MQ=255)
                                                                                                                                        ctctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                <  1:6989445/85‑1 (MQ=255)
                                                                                                                                        ctctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                >  2:6989349/1‑85 (MQ=255)
                                                                                                                                        ctctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg               >  2:2754993/1‑150 (MQ=255)
                                                                                                                                        ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg               >  1:3441639/1‑150 (MQ=255)
                                                                                                                                        ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg               >  1:4785/1‑150 (MQ=255)
                                                                                                                                         tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt              >  2:481581/1‑150 (MQ=255)
                                                                                                                                         tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt              >  2:5629454/1‑150 (MQ=255)
                                                                                                                                         tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt              >  1:920976/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTc                                                                                  <  1:3171928/82‑1 (MQ=25)
                                                                                                                                         tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTc                                                                                  >  2:3171928/1‑82 (MQ=25)
                                                                                                                                         tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                >  2:2978486/1‑84 (MQ=25)
                                                                                                                                         tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                                <  1:2978486/84‑1 (MQ=25)
                                                                                                                                         tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGt              >  2:2220725/1‑150 (MQ=255)
                                                                                                                                           ttcgcggcTCTAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGttg            >  2:3367077/4‑149 (MQ=255)
                                                                                                                                           ttcgcggcTCTAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc            >  1:4922657/4‑150 (MQ=255)
                                                                                                                                           ttcgcggcTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc            >  1:5842364/4‑150 (MQ=255)
                                                                                                                                           ttcgcggcTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc            >  1:6486265/4‑150 (MQ=255)
                                                                                                                                           ttcgcggcTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc            >  1:2219917/4‑150 (MQ=255)
                                                                                                                                           ttcgcggcTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc            >  1:1845007/4‑150 (MQ=255)
                                                                                                                                           ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc            >  2:5285736/1‑150 (MQ=255)
                                                                                                                                           ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc            >  2:1755043/1‑150 (MQ=255)
                                                                                                                                            tGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCt           <  2:6800941/150‑1 (MQ=255)
                                                                                                                                             ggcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTg          >  2:1875526/1‑150 (MQ=255)
                                                                                                                                              gcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTgg         >  1:695305/1‑150 (MQ=255)
                                                                                                                                              gcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTgg         <  1:6189612/150‑1 (MQ=255)
                                                                                                                                                ggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGaa       <  2:3333601/150‑1 (MQ=255)
                                                                                                                                                ggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGaa       <  1:67576/150‑1 (MQ=255)
                                                                                                                                                 gcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAg      <  2:4472487/150‑1 (MQ=255)
                                                                                                                                                     gagCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGcaca  >  2:657880/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACA  >  NC_002947/2956185‑2956483

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: