Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,187 C→T 5.1% intergenic (+14/+105) 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,967,1870CT 5.1% 304.1 / 4.2 97intergenic (+14/+105)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (37/55);  new base T (2/3);  total (39/58)
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.82e-01

CCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAA  >  NC_002947/1967040‑1967306
                                                                                                                                                   |                                                                                                                       
cccTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCG                                                                                                                       <  1:2321283‑M1/150‑1 (MQ=255)
cccTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCAAGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCG                                                                                                                       <  2:3144905‑M1/150‑1 (MQ=255)
cccTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCAAGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCG                                                                                                                       <  2:3144268‑M1/150‑1 (MQ=255)
  cTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGC                                                                                                                     <  2:5310923‑M1/150‑1 (MQ=255)
  cTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGC                                                                                                                     <  2:2679610‑M1/150‑1 (MQ=255)
  cTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGC                                                                                                                     <  1:5203852‑M1/150‑1 (MQ=255)
       gTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTT                                                                                                                >  1:4225428‑M1/1‑150 (MQ=255)
            ggACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGG                                                                                                           <  2:3968478‑M1/150‑1 (MQ=255)
            ggACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGG                                                                                                           <  2:3803771‑M1/150‑1 (MQ=255)
             gACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGG                                                                                                          >  2:2274134‑M1/1‑150 (MQ=255)
             gACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGG                                                                                                          <  1:3146483‑M1/150‑1 (MQ=255)
             gACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCAGG                                                                                                          <  2:945530‑M1/150‑1 (MQ=255)
              cccgAGGTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGC                                                                                                         <  2:3518898‑M1/149‑1 (MQ=255)
              aCCGAGTTGCGCAGGGCGATCTGCAAGAGTTTCGACTGCATGGTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGC                                                                                                         <  2:3519041‑M1/150‑1 (MQ=255)
                   gTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCC                                                                                                    <  2:149888‑M1/150‑1 (MQ=255)
                    ttGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCC                                                                                                   <  1:302439‑M1/150‑1 (MQ=255)
                    ttGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCC                                                                                                   <  1:302850‑M1/150‑1 (MQ=255)
                        gcAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                               >  2:3123692‑M1/1‑150 (MQ=255)
                        gcAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                               >  1:3140530‑M1/1‑150 (MQ=255)
                        gcAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGCATCGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                               >  2:1280934‑M1/1‑150 (MQ=255)
                        gcAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGCATCGGCCTCTTCGCGGGCTTGCCCGCTC                                                                                               >  2:1280915‑M1/1‑150 (MQ=255)
                          aGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCgaag                                                                                             <  2:640573‑M1/150‑5 (MQ=255)
                             gCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACA                                                                                          <  1:2812399‑M1/150‑1 (MQ=255)
                              cGATCAGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAG                                                                                         <  2:1781214‑M1/150‑1 (MQ=255)
                                 tCTGCAAGAGTTACGACTGCATGTCCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGAt                                                                                      <  1:3847937/150‑1 (MQ=255)
                                   tGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCa                                                                                    <  1:2669942/150‑1 (MQ=255)
                                    gCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCAc                                                                                   >  1:677424/1‑150 (MQ=255)
                                    gCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCAc                                                                                   <  1:2378731/150‑1 (MQ=255)
                                               gACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCcaca                                                                        >  2:1749011/1‑150 (MQ=255)
                                                aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                       <  2:4894649/150‑1 (MQ=255)
                                                aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                       <  1:5373755/150‑1 (MQ=255)
                                                aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                       <  2:2802188/150‑1 (MQ=255)
                                                aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGCATCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttga                                                                       <  2:2861315‑M1/150‑18 (MQ=255)
                                                  tGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAATCCttgagg                                                                     <  2:1183292/150‑6 (MQ=255)
                                                   gCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGa                                                                    <  1:4827935/150‑1 (MQ=255)
                                                    cATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGG                                                                                                           <  2:3906469‑M1/110‑1 (MQ=255)
                                                    cATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGG                                                                                                           >  1:3906469‑M1/1‑110 (MQ=255)
                                                      tgttcctgttccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTg                                                                 <  1:4515495/150‑1 (MQ=255)
                                                      tgttcctgttccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTg                                                                 <  1:4515583/150‑1 (MQ=255)
                                                      tgttcctgttccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTg                                                                 <  1:4704831/150‑1 (MQ=255)
                                                       gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                      <  1:491378/144‑1 (MQ=255)
                                                       gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                      >  2:3219105/1‑144 (MQ=255)
                                                       gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                      <  1:3219105/144‑1 (MQ=255)
                                                       gttcctgttcctGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAg                                                                      >  2:491378/1‑144 (MQ=255)
                                                        ttcctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcg                                                               <  1:1434924‑M1/150‑26 (MQ=255)
                                                        ttcctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcg                                                               <  2:704810‑M1/150‑26 (MQ=255)
                                                            tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                           <  1:3136174/150‑1 (MQ=255)
                                                            tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcct                                                           <  1:2159497/150‑1 (MQ=255)
                                                                ccTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcct                                                       >  1:4141020‑M2/1‑147 (MQ=255)
                                                                  tgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgt                                                     <  2:1690262‑M1/150‑36 (MQ=255)
                                                                  tgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgt                                                     <  1:2020968‑M1/150‑36 (MQ=255)
                                                                  tgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgt                                                     <  1:1794367‑M1/150‑36 (MQ=255)
                                                                  tgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgt                                                     <  1:3945098‑M1/150‑36 (MQ=255)
                                                                  tgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGatcaccgcaatccttgagggttgcggtgatcctgt                                                     <  1:3123692‑M1/150‑36 (MQ=255)
                                                                  tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                     <  2:129472‑M2/150‑6 (MQ=255)
                                                                  tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                     >  1:2919329‑M2/1‑145 (MQ=255)
                                                                  tGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                     >  1:3160343‑M2/1‑145 (MQ=255)
                                                                    gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTTTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgg                                                   >  1:1713441‑M2/1‑143 (MQ=255)
                                                                    gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgg                                                   >  1:1672148‑M2/1‑143 (MQ=255)
                                                                    gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGCATCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgg                                                   <  1:1280915‑M1/150‑38 (MQ=255)
                                                                    gAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGCATCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtgg                                                   <  1:1280934‑M1/150‑38 (MQ=255)
                                                                     aGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtggg                                                  <  2:4974738‑M1/150‑39 (MQ=255)
                                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtggga                                                 <  1:4102123‑M1/150‑40 (MQ=255)
                                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtggga                                                 <  1:5172459‑M1/150‑40 (MQ=255)
                                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                 >  2:677424‑M2/1‑141 (MQ=255)
                                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                 >  1:582655‑M2/1‑141 (MQ=255)
                                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                 >  2:3140530‑M2/1‑141 (MQ=255)
                                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                 >  1:1006619‑M2/1‑141 (MQ=255)
                                                                      gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                 >  1:1423016‑M2/1‑141 (MQ=255)
                                                                       cgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCAGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggag                                                >  2:2528114‑M2/1‑140 (MQ=255)
                                                                         cTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGC                                                                                                     >  1:77903‑M1/1‑95 (MQ=255)
                                                                         cTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGC                                                                                                     <  2:77903‑M1/95‑1 (MQ=255)
                                                                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                             >  1:510943‑M2/1‑137 (MQ=255)
                                                                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                             >  1:4562751‑M2/1‑137 (MQ=255)
                                                                          tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgg                                             >  1:129472‑M2/1‑137 (MQ=255)
                                                                           gcatgcaGCTGGACGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCTAGGATTGCGCCCTGcctgtgggagcggg                                            >  1:4305442‑M2/1‑136 (MQ=255)
                                                                           gcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcggg                                            >  1:4305268‑M2/1‑136 (MQ=255)
                                                                           gcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcggg                                            >  1:4305437‑M2/1‑136 (MQ=255)
                                                                              tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATcctgtgggagcgggcaa                                         >  1:1292428‑M2/1‑133 (MQ=255)
                                                                                 aGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtg                                                    <  2:4009602‑M1/136‑37 (MQ=255)
                                                                                 aGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtg                                                    >  1:4009602‑M1/1‑100 (MQ=255)
                                                                                  gCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGGACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtgggagcgggcaagccc                                     >  2:2245878‑M2/1‑129 (MQ=255)
                                                                                    tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                 <  2:2768206‑M2/136‑10 (MQ=255)
                                                                                    tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgtggga                                                 >  1:2768206‑M2/1‑127 (MQ=255)
                                                                                      gAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGGGTGGCCGCTACCCCCCGCGCCCCGCAACCCAAAAGGCATTCGCCCCGCCGGGGGGGGAGGGCCAGCCcgcgc                                 >  2:3808223/1‑149 (MQ=255)
                                                                                      gAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCGGccggtgggtgcgggcaagcccgcga                                 <  2:3808260‑M2/150‑26 (MQ=255)
                                                                                       aaGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtggga                                                 <  1:4992872‑M1/133‑40 (MQ=255)
                                                                                       aaGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGcagggcgcaatccttgagggttgcggtgatcctgtggga                                                 >  2:4992872‑M1/1‑94 (MQ=255)
                                                                                           gcgcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGC                                                                                                                     <  2:3119147‑M1/61‑1 (MQ=255)
                                                                                           gcgcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGC                                                                                                                     >  1:3119147‑M1/1‑61 (MQ=255)
                                                                                           gcgcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGC                                                                                                                     >  1:3119241‑M1/1‑61 (MQ=255)
                                                                                           gcgcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGC                                                                                                                     <  2:3119241‑M1/61‑1 (MQ=255)
                                                                                             gcCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGCCCACAGGATCACCGCCACCCGCACGGATTGCGCCCGGcctgggggcgcgggccagcccgcggcgaggcc                          <  2:3480730‑M2/150‑33 (MQ=255)
                                                                                                                  cggcggGGGGTGCGGGCTGAGCGCTGGCCGGGCCCGGCCTCTGCGCGGGCTGGCCCGCCCCCACCGGAGCCCCGCACCCCCCACGGATTGCGCCCGGCCTGGGGGAGCGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGtt     <  2:2978498/150‑1 (MQ=255)
                                                                                                                     cGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCATTCCCCCACCCCCAGGATGAACCCCACCCAAAAGGCTTTGGCCCGGCCCGGGGGGGCGGGCAAGCCCGCGGACAGGCCCGTGGCGGGGTATTGGGCGGTGac  >  2:4577478/1‑149 (MQ=255)
                                                                                                                                   tGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                     <  2:1853041‑M2/85‑6 (MQ=255)
                                                                                                                                   tGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGcctgt                                                     >  1:1853041‑M2/1‑80 (MQ=255)
                                                                                                                                                   |                                                                                                                       
CCCTGGGGTTGCGGACCGAGTTGCGCAGGGCGATCTGCAAGAGTTACGACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAA  >  NC_002947/1967040‑1967306

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: