Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,637,368 A→T 23.5% intergenic (+35/+31) 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,3680AT23.5% 292.1 / 84.6 140intergenic (+35/+31)pdxJ/PP_1437pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (47/60);  new base T (19/14);  total (66/74)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

GGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCC  >  NC_002947/1637222‑1637517
                                                                                                                                                    |                                                                                                                                                      
gggTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCC‑‑GTTGAAGTCACTCTCATGTTTGCGGCTTCGCGATGCTCGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGcc                                                                                                                                                     >  1:1862592/1‑150 (MQ=255)
      cAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGaa                                                                                                                                               <  1:1367467/150‑1 (MQ=255)
        aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCGAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGcc                                                                                                                                                     <  1:1849379/142‑1 (MQ=255)
        aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACTGGCCTCTTCGCGGGCTTGCCCGCGaaaa                                                                                                                                             <  1:1085269/150‑3 (MQ=255)
        aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                                             <  1:730132/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‑‑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)
                   aTTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAAAAgg                                                                                                                                           >  2:1401021/1‑141 (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)
                               cTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGCTTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCAc                                                                                                                      <  2:1862592/150‑1 (MQ=255)
                                tggtggCGCATGCCTTGTTTGTCGGCTAAATGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACt                                                                                                                     >  1:1024024/1‑150 (MQ=255)
                                   tggCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCgg                                                                                                                  <  1:1154811/150‑1 (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)
                                        cATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAggc                                                                                                             <  2:403064/150‑2 (MQ=255)
                                         aTGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGAt                                                                                                                         <  2:573472/137‑1 (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                                                                                                                                                  >  2:893581/1‑107 (MQ=255)
                                              ttgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGc                                                                                                                                          <  1:1998342/115‑1 (MQ=255)
                                                gtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAg                                                                                                     >  2:1703909/1‑150 (MQ=255)
                                                     tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGACAGTACTAGCGGCCGATCACTCGGCAGGAAACACCAGGGTaa                                                                                                >  1:1307540/1‑150 (MQ=255)
                                                     tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                                  <  1:1406197/100‑1 (MQ=37)
                                                     tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                                  >  2:1406197/1‑100 (MQ=25)
                                                      cGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTaaa                                                                                               <  2:689299/150‑1 (MQ=255)
                                                       ggCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAAc                                                                                              >  2:2066427/1‑150 (MQ=255)
                                                         cTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCg                                                                                            >  1:2319398/1‑150 (MQ=255)
                                                          ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCa                                                                                                                >  2:1397648/1‑129 (MQ=255)
                                                          ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCa                                                                                                                <  1:1397648/129‑1 (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                                                                                          <  2:2319398/150‑1 (MQ=255)
                                                           tCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                                          <  1:2109854/150‑1 (MQ=255)
                                                                gctgctGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGc                                                                                     <  1:2066427/150‑1 (MQ=255)
                                                                      gtgGCCGAGATGAAG‑‑GCGCTGGTCGTGGCGGCTTCGCGCTTATTCTTGGTCGCCTGTAGCAGTGCGTCCACACGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAgg                                                                               <  1:1224117/150‑1 (MQ=255)
                                                                         gCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGa                                                                                                                                                <  2:1605778/82‑1 (MQ=255)
                                                                         gCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAggcgg                                                                                                                              <  1:799740/100‑1 (MQ=255)
                                                                          ccGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGc                                                                           >  2:220442/1‑150 (MQ=255)
                                                                            gagaTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTa                                                                                                 <  2:1814279/126‑1 (MQ=255)
                                                                              gaTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgct                                                                       >  2:2278890/1‑150 (MQ=255)
                                                                                tGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                                     <  2:1680739/150‑1 (MQ=255)
                                                                                  aaG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGctgctg                                                                  <  1:2278314/150‑1 (MQ=255)
                                                                                  aaG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGctgctg                                                                  >  2:139502/1‑150 (MQ=255)
                                                                                    g‑‑gCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAc                                                                >  1:2053200/1‑150 (MQ=255)
                                                                                    g‑‑gCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAc                                                                >  2:734531/1‑150 (MQ=255)
                                                                                        cgcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCt                                                              <  1:1067883/150‑1 (MQ=255)
                                                                                           tggtCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCGGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGa                                                           <  1:194025/147‑1 (MQ=255)
                                                                                                 tGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgcc                                                     <  2:1849223/150‑1 (MQ=255)
                                                                                                   gcggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTTTTTTGTTTTCCCGCGAAGAGGCCAGCCCAGGCGGCCAATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGCcccgcccg                                                   >  2:341326/1‑150 (MQ=255)
                                                                                                     ggcTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCcgcccgcc                                                 >  1:1389823/1‑150 (MQ=255)
                                                                                                                  ttgttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAgg                                                                                                              <  1:389975/77‑1 (MQ=255)
                                                                                                                    gttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGAt                                  >  1:1757146/1‑150 (MQ=255)
                                                                                                                     tcggccGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGACAGCTGCTGACCTCGACCCGCCCGCCa                                                >  1:1036244/6‑135 (MQ=18)
                                                                                                                     tcggccGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGACAGCTGCTGACCTCGACCCGCCCGCCa                                                <  2:1036244/130‑1 (MQ=18)
                                                                                                                           gCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                                         <  2:1135473/136‑1 (MQ=25)
                                                                                                                           gCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                                         >  1:1135473/1‑136 (MQ=25)
                                                                                                                           gCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTc                           >  1:2210909/1‑150 (MQ=255)
                                                                                                                             ctgtgctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAGCCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCAACCCGCCCGCCATGCAACTCCATGATCGCCTTCAc                         >  1:378325/8‑150 (MQ=21)
                                                                                                                                 aCCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc                     <  2:1584577/150‑1 (MQ=255)
                                                                                                                                 aCCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc                     <  2:1024024/150‑1 (MQ=255)
                                                                                                                                 aCCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc                     <  2:1112380/150‑1 (MQ=255)
                                                                                                                                  ctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGGAACACCAGGGTAAACCggg                                                                                          <  2:2035720/79‑1 (MQ=17)
                                                                                                                                  ctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGAc                              <  1:557697/138‑1 (MQ=38)
                                                                                                                                  ctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGAc                              >  2:557697/3‑140 (MQ=38)
                                                                                                                                   tggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                                     >  1:1837740/2‑101 (MQ=25)
                                                                                                                                   tggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                                     <  2:1837740/100‑1 (MQ=25)
                                                                                                                                   tggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc                   >  1:2359080/2‑150 (MQ=255)
                                                                                                                                   cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATc                                 >  1:371722/1‑136 (MQ=255)
                                                                                                                                   cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATc                                 <  2:371722/136‑1 (MQ=255)
                                                                                                                                     gCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACAAGGGTAAACCGGGTAGGCCCCAggg                                                                              <  1:263445/90‑1 (MQ=25)
                                                                                                                                      ccTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgc                                                      >  1:2414161/1‑112 (MQ=255)
                                                                                                                                      ccTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgc                                                      <  2:2414955/112‑1 (MQ=255)
                                                                                                                                      ccTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAGCTGCATGCTCTACTGCACGATCGgcaa                >  2:1960852/1‑150 (MQ=255)
                                                                                                                                       ctcttCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc                     >  1:1029037/4‑144 (MQ=38)
                                                                                                                                       ctcttCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAAACATAAATGCCCCAGGGGGCAAATAATATGAAAAGAAATGCCGCCATTCAACTGCATGATCGACTTCACGATCGccaac               >  1:533968/4‑150 (MQ=255)
                                                                                                                                       cTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaac               >  2:1689518/1‑150 (MQ=255)
                                                                                                                                        ttttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaacc              >  1:163098/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaacccc            >  1:2024365/1‑150 (MQ=38)
                                                                                                                                               gggCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACtt                            <  1:178204/129‑1 (MQ=255)
                                                                                                                                               gggCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACtt                            >  2:178204/1‑129 (MQ=255)
                                                                                                                                                   aaGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAc                                                          >  2:1960287/1‑95 (MQ=255)
                                                                                                                                                    aGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTcc  >  2:1768326/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCC  >  NC_002947/1637222‑1637517

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