Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,972,180 C→T 7.6% intergenic (+58/+61) PP_1762 → / ← yciK hypothetical protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,972,1800CT 7.6% 557.7 / 25.1 183intergenic (+58/+61)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (87/82);  new base T (3/11);  total (90/93)
Fisher's exact test for biased strand distribution p-value = 4.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 184 total aligned reads displayed.

ACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCG  >  NC_002947/1972037‑1972329
                                                                                                                                               |                                                                                                                                                     
aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                                 <  1:6450663/150‑1 (MQ=255)
aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                                 <  1:567129/150‑1 (MQ=255)
aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                                 <  1:4919673/150‑1 (MQ=255)
aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                                 <  2:4910958/150‑1 (MQ=255)
 cTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                                                                                                <  2:5156828/150‑1 (MQ=255)
 cTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                                                                                                <  2:662475/150‑1 (MQ=255)
    gctgcccgaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGAt                                                                                                                                             <  1:4785376‑M1/54‑1 (MQ=255)
     aTGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGAtt                                                                                                                                            <  2:1670834/150‑1 (MQ=255)
       gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGTTTgc                                                                                                                                          >  2:6783251/1‑150 (MQ=255)
           gaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTTCCACAGGATCACCGCAACCCTCAAGGATTGCGccc                                                                                                                                      >  1:7223461‑M1/90‑150 (MQ=255)
           ctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGccc                                                                                                                                      <  1:669141/150‑1 (MQ=255)
           ctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGccc                                                                                                                                      <  1:2862576/150‑1 (MQ=255)
            tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                                                                                                     <  2:5168396/150‑1 (MQ=255)
            tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                                                                                                     <  1:3253055/150‑1 (MQ=255)
            aagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAACCCTCAAGGATTGCGCcct                                                                                                                                     >  1:1332461‑M1/89‑150 (MQ=255)
              gagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                                   >  1:6223121‑M1/87‑150 (MQ=255)
              agagatcatgccggtgtacctgtacctgatggggccgacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                                   >  1:1204382‑M1/87‑150 (MQ=255)
               agatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgcc                                                                                                                                  <  1:6308746‑M1/65‑1 (MQ=255)
                gatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgcct                                                                                                                                 >  2:442241‑M1/85‑150 (MQ=255)
                ccTGGAGCGCTGCATGCAGCCGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcct                                                                                                                                 >  2:4761570‑M2/1‑147 (MQ=255)
                  tcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                                               >  2:7342803‑M1/83‑150 (MQ=255)
                  tcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                                               >  1:698519‑M1/83‑150 (MQ=255)
                  tcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                                               >  2:3705024‑M1/83‑150 (MQ=255)
                  tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                                                               >  2:1430575‑M2/1‑145 (MQ=255)
                  tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                                                               >  1:3438349‑M2/1‑145 (MQ=255)
                  tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                                                                                               >  2:5547914‑M2/1‑145 (MQ=255)
                   catgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                                              >  2:262781‑M1/82‑150 (MQ=255)
                    atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                             >  2:4680864‑M1/81‑150 (MQ=255)
                    atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                             >  2:2104369‑M1/81‑150 (MQ=255)
                    atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                             >  2:6308746‑M1/81‑150 (MQ=255)
                     tgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTggg                                                                                                                            >  1:3272565‑M1/80‑150 (MQ=255)
                      gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                           >  1:1902327‑M1/79‑150 (MQ=255)
                      gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                           >  1:5367817‑M1/79‑150 (MQ=255)
                       ccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaaagcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAg                                                                                                                          >  1:5770741‑M1/78‑150 (MQ=255)
                       ccggtgtacctgtacctgagggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaatcaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAg                                                                                                                          >  2:5901677‑M1/78‑150 (MQ=255)
                        gcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagc                                                                                                                         >  1:6546342‑M2/1‑139 (MQ=255)
                        cggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGc                                                                                                                         >  2:6901145‑M1/77‑150 (MQ=255)
                        cggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGc                                                                                                                         >  2:4048487‑M1/77‑150 (MQ=255)
                        cggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCCTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGc                                                                                                                         >  1:223986‑M1/77‑150 (MQ=255)
                        cggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCCTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGAGTGCGCCCTGCCTGTGGGAGc                                                                                                                         >  2:2056320‑M1/77‑150 (MQ=255)
                         ggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCg                                                                                                                        >  2:4785376‑M1/76‑150 (MQ=255)
                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgg                                                                                                                       >  1:3586385‑M2/1‑137 (MQ=255)
                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgg                                                                                                                       >  2:2564707‑M2/1‑137 (MQ=255)
                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                                                                                                       >  2:2468407‑M2/1‑137 (MQ=255)
                          gtgtacctgtacctgatggggccggacaacacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                                                       >  2:6947531‑M1/75‑150 (MQ=255)
                            catgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCGTGGCCGCGCCCCCCCGGTCCCCCCCACCCCCACGGCTTGCGCCCGGCCGGGGGgggcgggc                                                                                                                     >  1:1041735/1‑150 (MQ=255)
                                caGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGAATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagc                                                                                                                 >  2:810346‑M2/1‑131 (MQ=255)
                                  gCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagccc                                                                                                               >  1:6108293‑M2/1‑129 (MQ=255)
                                  gCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagccc                                                                                                               >  2:5768636‑M2/1‑129 (MQ=255)
                                  gCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagccc                                                                                                               >  1:3246767‑M2/1‑129 (MQ=255)
                                         gggCGCCAGCACGGCCGTATGCGGCCGGCGTTGGGGGCGGGCTGATCGCTTTCCTGTACCGGCCTCTGCGCGGGCTTGCCCGCTCCCCCAGGCTGGGCGAAATCCTTGTGGGTTGCGTTGTTCGTGTGGTAGCTGGAAAGCCCGCGAaga                                                                                                        <  2:3756800/150‑1 (MQ=255)
                                              gccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                           >  1:1159993‑M1/55‑126 (MQ=255)
                                                         gTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggctt                                                                                                                   <  2:1348725‑M2/123‑18 (MQ=255)
                                                                        aggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCTGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                                   >  2:5866172‑M1/29‑92 (MQ=255)
                                                                        aggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                                   <  1:5866172‑M1/64‑1 (MQ=255)
                                                                                                   cgggCCCCTGCGCGGGCTGGCCCGCACCCCCCGTATCCCCGCAACCCTCACGTTGTGCGCCGGGCTGGGGGGGGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCgg                                              <  1:3509795/148‑1 (MQ=255)
                                                                                                           ttCGCGGGCGTGCCCGCGCCCCCCGGATCACCCCAACCCGCAAGGCTTGCGCCCGGCCCGTGGGCGCGGGCCAGCCCCCGCACAGGCCGCTGCAGGTTTATTGGGCGTTGAACCCCAGCCCGTGCACCCCGGGGCTGTCCGGcccccgcc                                      >  2:5122831/1‑146 (MQ=255)
                                                                                                                  GCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacag                               <  2:3099938‑M1/150‑1 (MQ=255)
                                                                                                                   CTTCCCCGCTCCCACACGAACACCCCAACCCACACGCATTGCCCCCCCCCTGTGGGAGCGGCAAACCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacagg                              >  1:6783251‑M1/1‑150 (MQ=255)
                                                                                                                    TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGCTGCAGGGGGATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                             <  2:4986931‑M1/150‑1 (MQ=255)
                                                                                                                    TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                             <  2:1052131‑M1/150‑1 (MQ=255)
                                                                                                                    TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                             <  2:2708537‑M1/150‑1 (MQ=255)
                                                                                                                    TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                             <  1:4508728‑M1/150‑1 (MQ=255)
                                                                                                                    TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                             <  1:3097626‑M1/150‑1 (MQ=255)
                                                                                                                    TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                             <  2:2179183‑M1/150‑1 (MQ=255)
                                                                                                                    TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTGTATTGGGCGTTGAACGCCTGCCCGGTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                             <  2:332553‑M1/150‑1 (MQ=255)
                                                                                                                    TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATccggtccaggt                             <  2:1248039‑M1/150‑1 (MQ=255)
                                                                                                                        ccGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccgatacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgca                         <  2:1274426‑M2/150‑108 (MQ=255)
                                                                                                                        CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                         <  1:4315384‑M1/150‑1 (MQ=255)
                                                                                                                        CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                         <  2:6015058‑M1/150‑1 (MQ=255)
                                                                                                                        CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                         <  1:6779468‑M1/150‑1 (MQ=255)
                                                                                                                            tCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgc                     <  2:6726053‑M2/150‑112 (MQ=255)
                                                                                                                            TCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc                     <  2:3577605‑M1/150‑1 (MQ=255)
                                                                                                                            TCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc                     <  1:2832869‑M1/150‑1 (MQ=255)
                                                                                                                             cccACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccgatacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgct                    <  1:1712675‑M2/150‑113 (MQ=255)
                                                                                                                              CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt                   <  1:6625135‑M1/150‑1 (MQ=255)
                                                                                                                               cacaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcc                  <  1:1023074‑M2/150‑115 (MQ=255)
                                                                                                                               cacaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcc                  <  2:1922622‑M2/150‑115 (MQ=255)
                                                                                                                               CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg                  <  1:2436305‑M1/150‑1 (MQ=255)
                                                                                                                               CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg                  <  1:2932991‑M1/150‑1 (MQ=255)
                                                                                                                                acaGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagcccgcgaagaggccggtacaggaaagcgatcagcccgcagccaacgccggccgcatacggccttgctggcgcccttccagctgcatgcagcgctcca                 <  1:6505134‑M2/150‑116 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTc                                                               >  1:2060682‑M1/1‑104 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                 <  1:264708‑M1/150‑1 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                 <  2:5364966‑M1/150‑1 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                 <  2:5732971‑M1/150‑1 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                 <  1:5012159‑M1/150‑1 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                 <  1:3385582‑M1/150‑1 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                 <  1:6129008‑M1/150‑1 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGa                 <  2:6588988‑M1/150‑1 (MQ=255)
                                                                                                                                ACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTGTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTCCAGGTACACCGGCATGa                 <  2:4607888‑M1/150‑1 (MQ=255)
                                                                                                                                 CAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAt                >  1:4765556‑M1/1‑150 (MQ=255)
                                                                                                                                    GATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAtctc             <  1:4983645‑M1/150‑1 (MQ=255)
                                                                                                                                      tCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGTGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCTTCAGGTACAGGTACACCGGCATGATCTCtt           >  2:6440786/1‑150 (MQ=255)
                                                                                                                                      tCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCtt           >  2:6193872/1‑150 (MQ=255)
                                                                                                                                       cACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTc          >  2:964717/1‑150 (MQ=255)
                                                                                                                                         ccGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCgg        >  2:3235303/1‑150 (MQ=255)
                                                                                                                                         ccGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCgg        >  2:3407039/1‑150 (MQ=255)
                                                                                                                                         ccGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCgg        >  1:452941/1‑150 (MQ=255)
                                                                                                                                             aaCCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAg    >  2:6342091/1‑150 (MQ=255)
                                                                                                                                               cccTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCg  >  2:7253561/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
ACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCG  >  NC_002947/1972037‑1972329

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