Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,642,327 G→C 15.7% intergenic (‑113/+39) PP_1442 ← / ← lon‑I hypothetical protein/DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,642,3270GC15.7% 337.7 / 49.0 121intergenic (‑113/+39)PP_1442/lon‑Ihypothetical protein/DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base G (57/45);  new base C (11/8);  total (68/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 124 total aligned reads displayed.

AGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAG  >  NC_002947/1642181‑1642476
                                                                                                                                                  |                                                                                                                                                     
aGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCACCTGCACGGGTACGCTACCCCTGc                                                                                                                                                    <  2:2531809/150‑1 (MQ=39)
 gCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCa                                                                                                                                                   <  2:4524133/150‑1 (MQ=38)
 gCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGACTGTGCCGGTCCTTTCGCGAGAAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCa                                                                                                                                                   >  1:941341/1‑150 (MQ=255)
   tCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAgg                                                                                                                                                 >  1:2198244/1‑150 (MQ=255)
    cAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGt                                                                                                                                                >  1:1295855/1‑150 (MQ=255)
      ggggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGc                                                                                                                                              <  1:3983714/150‑1 (MQ=255)
      ggggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGCCCTCTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGgcgc                                                                                                                                              <  2:468026/150‑1 (MQ=21)
       gggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCg                                                                                                                                             >  1:3153143/1‑150 (MQ=255)
         gAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCTATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCCTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCggg                                                                                                                                           <  2:159810/150‑1 (MQ=25)
         gAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCggg                                                                                                                                           >  1:3723583/1‑150 (MQ=255)
         gAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCggg                                                                                                                                           >  2:1183581/1‑150 (MQ=255)
          aaCGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTACCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGc                                                                                                                                          <  1:592120/150‑1 (MQ=255)
             gAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTg                                                                                                                                       >  2:1367583/1‑150 (MQ=255)
                  gACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCcgc                                                                                                                                  >  2:4959550/1‑150 (MQ=255)
                     gaggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGaa                                                                                                                               >  1:2403342/1‑150 (MQ=255)
                      aggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGaaa                                                                                                                              <  2:2369297/150‑2 (MQ=21)
                      aggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCCTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGaaa                                                                                                                              <  1:4386561/150‑2 (MQ=21)
                       ggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAaga                                                                                                                             <  2:3937339/150‑1 (MQ=255)
                             gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGc                                                                                                                                                    >  1:5034824/1‑121 (MQ=255)
                             gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGc                                                                                                                                                    <  2:5034824/121‑1 (MQ=255)
                             gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCgg                                                                                                                       >  1:1293465/1‑150 (MQ=255)
                             gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTACCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCgg                                                                                                                       >  2:1101257/1‑150 (MQ=255)
                              tCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGt                                                                                                                      >  2:4866549/1‑150 (MQ=255)
                                      gCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCTTACCGGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTa                                                                                                              <  1:2732561/150‑1 (MQ=255)
                                       cTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAg                                                                                                             <  2:1293465/150‑1 (MQ=255)
                                          cGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAAc                                                                                                          >  1:4859273/1‑150 (MQ=255)
                                           gAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                         >  2:466024/1‑150 (MQ=255)
                                                 tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt                                                                                                   >  2:2186524/1‑150 (MQ=255)
                                                   gCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc                                                                                                 <  1:3505890/150‑1 (MQ=255)
                                                   gCGATGGCCTGAGGATACCCGATGATCGAGGAAATTAGTGCTGCCTGTGCCGGTCCTTTAGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc                                                                                                 >  2:4910293/1‑150 (MQ=255)
                                                       tGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCa                                                                                             <  1:4081226/150‑1 (MQ=255)
                                                            tGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGc                                                                                                                                          >  1:4582232/1‑100 (MQ=255)
                                                            tGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGc                                                                                                                                          <  2:4582232/100‑1 (MQ=255)
                                                                  tACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAg                                                                                  <  2:1384508/150‑1 (MQ=255)
                                                                     ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                               >  1:3698748/1‑150 (MQ=255)
                                                                     ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCAAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                               >  2:3561484/1‑150 (MQ=255)
                                                                     ccGATGATCGAGGAAATTTGTGCTGCCCGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGTTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                               >  2:1351163/1‑150 (MQ=255)
                                                                     ccGATGATCGAGGAAATGTGTGCTGCCTGTGCCGCTCCTTTCGCGAGTAAGCCCGCGCCTCCAGGGGTAGCGTCCCCGTCCAGGTGCGGGCGAGCCCGCGAAGAGGGCGGTGCACGCTAGAACAGCACTGACGCCACGACCCCAAAgcgc                                                                               >  2:2102370/1‑150 (MQ=255)
                                                                       gatgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCtt                                                                             >  1:439685/1‑150 (MQ=255)
                                                                       gatgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCtt                                                                             >  2:3340720/1‑150 (MQ=255)
                                                                        atgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTg                                                                            <  2:1756465/150‑1 (MQ=255)
                                                                         tgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTgg                                                                           <  1:697089/150‑1 (MQ=255)
                                                                                 gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGc                                                                   <  1:4835788/150‑1 (MQ=255)
                                                                                   aaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc                                                                 <  1:2677061/150‑1 (MQ=255)
                                                                                   aaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc                                                                 <  1:2677108/150‑1 (MQ=255)
                                                                                     tttGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTg                                                               >  1:628214/1‑150 (MQ=255)
                                                                                       tgtgCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTc                                                             >  2:3077553/1‑150 (MQ=255)
                                                                                        gtgCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCa                                                            <  2:960249/150‑1 (MQ=255)
                                                                                         tgctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAg                                                           >  2:2881765/1‑150 (MQ=255)
                                                                                          gctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAgg                                                          <  1:1637896/150‑1 (MQ=255)
                                                                                          gctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAgg                                                          <  1:1637992/150‑1 (MQ=255)
                                                                                             gcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCt                                                       >  1:744741/1‑150 (MQ=255)
                                                                                                tgtgCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcg                                                    >  1:2914734/1‑150 (MQ=255)
                                                                                                      ggTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAgg                                              >  1:3238260/1‑150 (MQ=255)
                                                                                                      ggTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAgg                                              <  2:3845355/150‑1 (MQ=255)
                                                                                                        tCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTa                                            >  1:4867641/1‑150 (MQ=255)
                                                                                                        tCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTa                                            <  1:3627836/150‑1 (MQ=255)
                                                                                                           tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGgtcgtc                                         >  2:3551093/1‑150 (MQ=255)
                                                                                                           tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc                                         >  2:2208301/1‑150 (MQ=255)
                                                                                                           tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc                                         >  2:818383/1‑150 (MQ=255)
                                                                                                            ttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTcc                                        >  2:4518818/1‑150 (MQ=255)
                                                                                                               gcgAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGc                                     <  1:2532963/150‑1 (MQ=255)
                                                                                                                cgAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa                                    <  2:2745271/150‑1 (MQ=255)
                                                                                                                 gggcaaGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAg                                   >  2:294256/5‑150 (MQ=37)
                                                                                                                 gAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAg                                   <  1:1512179/150‑1 (MQ=255)
                                                                                                                 gAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAg                                   >  2:3064702/1‑150 (MQ=255)
                                                                                                                       gCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                             <  2:439685/150‑1 (MQ=255)
                                                                                                                       gCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                             <  2:4467953/150‑1 (MQ=255)
                                                                                                                       gCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                             <  1:294256/150‑1 (MQ=38)
                                                                                                                        cccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTCGTCCGGCAGGTCTTCg                            >  2:282686/1‑150 (MQ=255)
                                                                                                                        cccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg                            >  1:2203253/1‑150 (MQ=255)
                                                                                                                        cccGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg                            >  1:1283329/1‑150 (MQ=37)
                                                                                                                         ccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGa                           >  1:2144286/1‑150 (MQ=255)
                                                                                                                         ccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGa                           >  2:393854/1‑150 (MQ=255)
                                                                                                                         ccGCACCTGCACGGGTACGCTACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGa                           >  1:1184058/1‑150 (MQ=38)
                                                                                                                          cgcgCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                          <  2:3153143/150‑1 (MQ=255)
                                                                                                                          cgcgCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                          <  1:622725/150‑1 (MQ=255)
                                                                                                                          cgcgCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                          <  2:2011465/150‑1 (MQ=255)
                                                                                                                          cgcaCCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                          >  1:4747621/1‑150 (MQ=37)
                                                                                                                          cgcaCCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                          >  2:99480/1‑150 (MQ=38)
                                                                                                                          cgcaCCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                          >  2:99406/1‑150 (MQ=38)
                                                                                                                           gcaccTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAg                         >  2:94957/4‑150 (MQ=38)
                                                                                                                            cgcCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGt                        <  2:3723583/150‑1 (MQ=255)
                                                                                                                               cTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTTGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgc                     >  1:2957269/1‑150 (MQ=255)
                                                                                                                                  cAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgccgc                  >  1:1053930/1‑150 (MQ=255)
                                                                                                                                    ggggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCTCTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCgg                <  2:2450923/150‑1 (MQ=255)
                                                                                                                                      ggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGtt              <  2:1295855/150‑1 (MQ=255)
                                                                                                                                       gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg             >  2:1050826/1‑150 (MQ=255)
                                                                                                                                       gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg             <  1:4576390/150‑1 (MQ=255)
                                                                                                                                         aGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTggg           <  1:816923/150‑1 (MQ=255)
                                                                                                                                           cGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCt         >  1:2685525/1‑150 (MQ=255)
                                                                                                                                            gTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTc        >  2:1842630/1‑150 (MQ=255)
                                                                                                                                             tACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                         >  1:3546361/1‑52 (MQ=17)
                                                                                                                                             tACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                         <  2:3546361/52‑1 (MQ=17)
                                                                                                                                              aCCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCg      <  1:798754/150‑1 (MQ=255)
                                                                                                                                              aCCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCg      >  2:2603525/1‑150 (MQ=255)
                                                                                                                                                ccGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGc    >  2:2588821/1‑150 (MQ=255)
                                                                                                                                                 cGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCa   <  2:3338834/150‑1 (MQ=255)
                                                                                                                                                  gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAg  >  1:3435272/1‑150 (MQ=255)
                                                                                                                                                  gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAg  >  2:2439921/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAG  >  NC_002947/1642181‑1642476

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: