Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,637,349 A→G 7.9% intergenic (+16/+50) pdxJ → / ← PP_1437 pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,637,3490AG 7.9% 317.0 / 11.1 112intergenic (+16/+50)pdxJ/PP_1437pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (44/58);  new base G (4/5);  total (48/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

Only 100 of 115 total aligned reads displayed.

GGTGGCGGCGATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATC  >  NC_002947/1637207‑1637498
                                                                                                                                                |                                                                                                                                                     
ggTGGCGGCGATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCt                                                                                                                                                <  1:2067634/150‑1 (MQ=255)
      ggcgATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTcgc                                                                                                                                          >  1:1218138/1‑150 (MQ=255)
        cgATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCgg                                                                                                                                        <  1:362559/150‑1 (MQ=255)
         gATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCggg                                                                                                                                       >  1:1398882/1‑150 (MQ=255)
               gggTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCC‑‑GTTGAAGTCACTCTCATGTTTGCGGCTTCGCGATGCTCGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGcc                                                                                                                                 >  1:1862592/1‑150 (MQ=255)
                     cAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGaa                                                                                                                           <  1:1367467/150‑1 (MQ=255)
                       aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCGAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCCTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGcc                                                                                                                                 >  2:1849379/1‑142 (MQ=255)
                       aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                         <  1:730132/150‑1 (MQ=255)
                       aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                         <  1:1038958/150‑1 (MQ=255)
                       aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGaaaa                                                                                                                         <  2:753733/150‑3 (MQ=255)
                       aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCCGCCTTTTCCCGGGCAAGCCCGCGCaga                                                                                                                         >  2:1038958/1‑150 (MQ=255)
                        tGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAagag                                                                                                                        >  1:960530/1‑150 (MQ=255)
                         gAGTTGAACATTGGCCATGCGCTGGTGGCGCCTGCCTTGTTTGTCGGTTCCAGCACGGCTGGGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAgg                                                                                                                       <  2:1244676/150‑1 (MQ=255)
                         gAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAgg                                                                                                                       <  1:504279/150‑1 (MQ=255)
                         gAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAgg                                                                                                                       <  2:2440935/150‑1 (MQ=38)
                          aGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGc                                                                                                                      <  1:2137658/150‑1 (MQ=255)
                            ttGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCa                                                                                                                    <  1:483740/150‑1 (MQ=255)
                            ttGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCa                                                                                                                    <  1:254890/150‑1 (MQ=255)
                             tGAACATTGGCCATGCGCTGGTGGCGCATGCGTGGTGTGTCGGCGTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGGGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGCGGCCAg                                                                                                                   <  2:493187/150‑1 (MQ=255)
                               aaCATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGca                                                                                                                 <  2:960530/150‑1 (MQ=255)
                                  aTTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGGTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAg                                                                                                              >  2:1046645/1‑150 (MQ=25)
                                  aTTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAAAAgg                                                                                                                       <  1:1401021/141‑1 (MQ=255)
                                       ccATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcggc                                                                                                         <  2:1791401/150‑1 (MQ=255)
                                       ccATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggaggc                                                                                                         >  1:403064/1‑150 (MQ=255)
                                       ccATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGCCGCCTGTGCTGGTCTCTTCGCTTTCTTGCCCGCGAAAAGGCCGGTACAggcgac                                                                                                         <  1:1046645/150‑3 (MQ=25)
                                        cATGCGCTGGTGGCGCGTGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGACGCCAGCACTGGCGGcc                                                                                                        >  2:1224117/1‑150 (MQ=255)
                                            cgcTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCNATc                                                                                                    <  2:1503247/150‑1 (MQ=255)
                                               tggtggCGCATGCCTTGTTTGTCGGCTAAATGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACt                                                                                                 >  1:1024024/1‑150 (MQ=255)
                                                  tggCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCgg                                                                                              <  1:1154811/150‑1 (MQ=255)
                                                   ggCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCggg                                                                                                                                       <  1:123559/108‑1 (MQ=255)
                                                   ggCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCggg                                                                                                                                       >  2:123559/1‑108 (MQ=255)
                                                    gcgcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCa                                                                                            <  2:2347205/150‑1 (MQ=255)
                                                     cgcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAg                                                                                           <  1:1160394/150‑1 (MQ=255)
                                                      gcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAaaag                                                                                                                        >  1:1830016/1‑120 (MQ=39)
                                                      gcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAaaag                                                                                                                        <  2:1830016/120‑1 (MQ=39)
                                                       cATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAggc                                                                                         <  2:403064/150‑2 (MQ=255)
                                                         tGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGaaa                                                                                       >  1:1112380/1‑150 (MQ=255)
                                                         tGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGaaa                                                                                       >  2:2109854/1‑150 (MQ=255)
                                                            cTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAacac                                                                                    <  1:2267949/150‑1 (MQ=255)
                                                            cTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAacac                                                                                    >  1:607784/1‑150 (MQ=255)
                                                             ttgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcgc                                                                                                                              <  1:893581/107‑1 (MQ=255)
                                                             ttgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcgc                                                                                                                              >  2:893581/1‑107 (MQ=255)
                                                             ttgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGc                                                                                                                      <  1:1998342/115‑1 (MQ=255)
                                                             ttgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGc                                                                                                                      >  2:1998342/1‑115 (MQ=255)
                                                               gtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAg                                                                                 <  2:2070435/150‑1 (MQ=255)
                                                               gtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAg                                                                                 <  1:1703909/150‑1 (MQ=255)
                                                                    tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                              >  2:1406197/1‑100 (MQ=25)
                                                                     cGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTaaa                                                                           <  2:689299/150‑1 (MQ=255)
                                                                      ggCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAAc                                                                          >  2:2066427/1‑150 (MQ=255)
                                                                         ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCa                                                                                            <  1:1397648/129‑1 (MQ=255)
                                                                         ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCa                                                                                            >  2:1397648/1‑129 (MQ=255)
                                                                         ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCgg                                                                       <  1:74856/150‑1 (MQ=255)
                                                                         ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCgg                                                                       >  1:587200/1‑150 (MQ=255)
                                                                         ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCgg                                                                       >  1:1043203/1‑150 (MQ=255)
                                                                         ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCgg                                                                       >  1:1680739/1‑150 (MQ=255)
                                                                          tCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                      <  1:2109854/150‑1 (MQ=255)
                                                                          tCAAGGCGGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                      <  2:1566056/150‑1 (MQ=255)
                                                                            aaGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCggg                                                                                                                                       <  1:1627269/83‑1 (MQ=255)
                                                                            aaGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCggg                                                                                                                                       >  2:1627269/1‑83 (MQ=255)
                                                                            aaGGCTGCTGTGGCCGAGATGAAC‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTa                                                                    <  1:1052165/150‑1 (MQ=255)
                                                                               gctgctGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGc                                                                 <  1:2066427/150‑1 (MQ=255)
                                                                                     gtgGCCGAGATGAAG‑‑GCGCTGGTCGTGGCGGCTTCGCGCTTATTCTTGGTCGCCTGTAGCAGTGCGTCCACACGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAgg                                                           <  1:1224117/150‑1 (MQ=255)
                                                                                        gCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGa                                                                                                                            >  1:1605778/1‑82 (MQ=255)
                                                                                        gCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGa                                                                                                                            <  2:1605778/82‑1 (MQ=255)
                                                                                        gCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcgg                                                                                                          <  1:799740/100‑1 (MQ=255)
                                                                                        gCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcgg                                                                                                          >  2:799740/1‑100 (MQ=255)
                                                                                         ccGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGc                                                       >  2:220442/1‑150 (MQ=255)
                                                                                           gagaTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTa                                                                             <  2:1814279/126‑1 (MQ=255)
                                                                                           gagaTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTa                                                                             >  1:1814279/1‑126 (MQ=255)
                                                                                             gaTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgct                                                   >  2:2278890/1‑150 (MQ=255)
                                                                                             gaTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgct                                                   >  2:2278314/1‑150 (MQ=255)
                                                                                              aTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctg                                                  <  1:316000/150‑1 (MQ=255)
                                                                                               tGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                 <  2:1680739/150‑1 (MQ=255)
                                                                                                 aaG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGctgctg                                               >  2:139502/1‑150 (MQ=255)
                                                                                                 aaG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGctgctg                                               <  1:2278890/150‑1 (MQ=255)
                                                                                                 aaG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGctgctg                                               <  1:2278314/150‑1 (MQ=255)
                                                                                                   g‑‑gCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAc                                             >  2:734531/1‑150 (MQ=255)
                                                                                                   g‑‑gCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAc                                             >  1:2053200/1‑150 (MQ=255)
                                                                                                       cgcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCt                                           <  1:1067883/150‑1 (MQ=255)
                                                                                                          tggtCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCGGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGa                                        <  1:194025/147‑1 (MQ=255)
                                                                                                             tCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGa                                                                                                      >  1:1338902/1‑85 (MQ=255)
                                                                                                             tCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGa                                                                                                      <  2:1338902/85‑1 (MQ=255)
                                                                                                                tGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgcc                                  <  2:1849223/150‑1 (MQ=255)
                                                                                                                tGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgcc                                  <  2:1398882/150‑1 (MQ=255)
                                                                                                                  gcggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTTTTTTGTTTTCCCGCGAAGAGGCCAGCCCAGGCGGCCAATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGCcccgcccg                                >  2:341326/1‑150 (MQ=255)
                                                                                                                    ggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCcgcccgcc                              >  1:1389823/1‑150 (MQ=255)
                                                                                                                                aTTGTTGGCCCCCTGTTCCGGCTTCTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCAt                  <  1:2250699/150‑1 (MQ=255)
                                                                                                                                 ttgttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAgg                                                                                          >  2:389975/1‑77 (MQ=255)
                                                                                                                                 ttgttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAgg                                                                                          <  1:389975/77‑1 (MQ=255)
                                                                                                                                   gttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGAt               >  1:1757146/1‑150 (MQ=255)
                                                                                                                                    ttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAGTTGCATggcg              <  2:1043203/150‑4 (MQ=255)
                                                                                                                                    tcggccGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGACAGCTGCTGACCTCGACCCGCCCGCCa                             <  2:1036244/130‑1 (MQ=18)
                                                                                                                                    tcggccGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGACAGCTGCTGACCTCGACCCGCCCGCCa                             >  1:1036244/6‑135 (MQ=18)
                                                                                                                                          gCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                      >  1:1135473/1‑136 (MQ=25)
                                                                                                                                          gCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                      <  2:1135473/136‑1 (MQ=25)
                                                                                                                                          gCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTc        >  1:2210909/1‑150 (MQ=255)
                                                                                                                                           ccTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCa       >  2:1789699/1‑150 (MQ=255)
                                                                                                                                            cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCAc      >  2:561778/1‑150 (MQ=255)
                                                                                                                                                aCCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc  <  2:1112380/150‑1 (MQ=255)
                                                                                                                                                aCCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc  <  2:1584577/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGTGGCGGCGATCAAGGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATC  >  NC_002947/1637207‑1637498

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