Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,637,367 A→T 23.5% intergenic (+34/+32) 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,3670AT23.5% 282.5 / 86.0 140intergenic (+34/+32)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 = 9.39e-01

Only 100 of 142 total aligned reads displayed.

GGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCC  >  NC_002947/1637222‑1637507
                                                                                                                                                   |                                                                                                                                            
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‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGaaaa                                                                                                                                  <  2:753733/150‑3 (MQ=255)
        aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCCGCCTTTTCCCGGGCAAGCCCGCGCaga                                                                                                                                  >  2:1038958/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)
              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                                                                                                                                >  2:1401021/1‑141 (MQ=255)
                   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)
                             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)
                                      cgcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAg                                                                                                    <  1:1160394/150‑1 (MQ=255)
                                       gcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAaaag                                                                                                                                 <  2:1830016/120‑1 (MQ=39)
                                        cATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAggc                                                                                                  <  2:403064/150‑2 (MQ=255)
                                         aTGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGAt                                                                                                              >  1:573472/1‑137 (MQ=255)
                                         aTGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGAt                                                                                                              <  2:573472/137‑1 (MQ=255)
                                          tGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGaaa                                                                                                >  1:1112380/1‑150 (MQ=255)
                                          tGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGaaa                                                                                                >  2:2109854/1‑150 (MQ=255)
                                             cTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAacac                                                                                             >  1:607784/1‑150 (MQ=255)
                                              ttgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcgc                                                                                                                                       >  2:893581/1‑107 (MQ=255)
                                              ttgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCcgc                                                                                                                                       <  1:893581/107‑1 (MQ=255)
                                              ttgtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGc                                                                                                                               >  2:1998342/1‑115 (MQ=255)
                                                gtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAg                                                                                          >  2:1703909/1‑150 (MQ=255)
                                                gtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAg                                                                                          <  2:2070435/150‑1 (MQ=255)
                                                gtttgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAg                                                                                          <  1:1703909/150‑1 (MQ=255)
                                                     tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGACAGTACTAGCGGCCGATCACTCGGCAGGAAACACCAGGGTaa                                                                                     >  1:1307540/1‑150 (MQ=255)
                                                     tCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                       <  1:1406197/100‑1 (MQ=37)
                                                      cGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTaaa                                                                                    <  2:689299/150‑1 (MQ=255)
                                                          ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCa                                                                                                     >  2:1397648/1‑129 (MQ=255)
                                                          ttCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCgg                                                                                <  1:74856/150‑1 (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)
                                                           tCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                               <  2:2319398/150‑1 (MQ=255)
                                                           tCAAGGCGGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                               <  2:1566056/150‑1 (MQ=255)
                                                                gctgctGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGc                                                                          <  1:2066427/150‑1 (MQ=255)
                                                                      gtgGCCGAGATGAAG‑‑GCGCTGGTCGTGGCGGCTTCGCGCTTATTCTTGGTCGCCTGTAGCAGTGCGTCCACACGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAgg                                                                    <  1:1224117/150‑1 (MQ=255)
                                                                       tgGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggg                                                                   <  2:1218138/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                                                                                      >  1:1814279/1‑126 (MQ=255)
                                                                              gaTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgct                                                            >  2:2278890/1‑150 (MQ=255)
                                                                               aTGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctg                                                           <  1:316000/150‑1 (MQ=255)
                                                                                tGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                          <  2:1680739/150‑1 (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)
                                                                                        cgcTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCt                                                    <  1:1067883/150‑1 (MQ=255)
                                                                                           tggtCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCGGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGa                                                 <  1:194025/147‑1 (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)
                                                                                                                     ttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAGTTGCATggcg                       <  2:1043203/150‑4 (MQ=255)
                                                                                                                     tcggccGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGACAGCTGCTGACCTCGACCCGCCCGCCa                                      >  1:1036244/6‑135 (MQ=18)
                                                                                                                           gCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                               <  2:1135473/136‑1 (MQ=25)
                                                                                                                           gCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTc                 >  1:2210909/1‑150 (MQ=255)
                                                                                                                            ccTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCa                >  2:1789699/1‑150 (MQ=255)
                                                                                                                             ctgtgctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAGCCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCAACCCGCCCGCCATGCAACTCCATGATCGCCTTCAc               >  1:378325/8‑150 (MQ=21)
                                                                                                                             cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCAc               >  2:561778/1‑150 (MQ=255)
                                                                                                                                tgctggCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggg                                                                   <  1:721653/91‑1 (MQ=255)
                                                                                                                                 aCCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc           <  2:1024024/150‑1 (MQ=255)
                                                                                                                                  ctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGGAACACCAGGGTAAACCggg                                                                               >  1:2035720/3‑81 (MQ=17)
                                                                                                                                  ctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGGAACACCAGGGTAAACCggg                                                                               <  2:2035720/79‑1 (MQ=17)
                                                                                                                                  ctggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGAc                    <  1:557697/138‑1 (MQ=38)
                                                                                                                                   tggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                          <  2:1837740/100‑1 (MQ=25)
                                                                                                                                   tggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                          >  1:1837740/2‑101 (MQ=25)
                                                                                                                                   tggCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc         >  1:2359080/2‑150 (MQ=255)
                                                                                                                                   cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATc                       <  2:371722/136‑1 (MQ=255)
                                                                                                                                   cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc         >  1:379719/1‑150 (MQ=255)
                                                                                                                                     gCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACAAGGGTAAACCGGGTAGGCCCCAggg                                                                   <  1:263445/90‑1 (MQ=25)
                                                                                                                                      ccTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgc                                            >  1:2414955/1‑112 (MQ=255)
                                                                                                                                      ccTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgc                                            >  1:2414161/1‑112 (MQ=255)
                                                                                                                                      ccTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgc                                            <  2:2414161/112‑1 (MQ=255)
                                                                                                                                      ccTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcccgc                                            <  2:2414955/112‑1 (MQ=255)
                                                                                                                                      acTTTTCGCGGGCAAGCCCGAGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaa      >  2:1696275/2‑150 (MQ=255)
                                                                                                                                       ctcttCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc           >  1:1029037/4‑144 (MQ=38)
                                                                                                                                       ctcttCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATc           <  2:1029037/141‑1 (MQ=38)
                                                                                                                                       cTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaac     <  1:139502/150‑1 (MQ=255)
                                                                                                                                        ttttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaacc    >  1:163098/1‑150 (MQ=255)
                                                                                                                                         tttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaaccc   <  2:163098/150‑1 (MQ=255)
                                                                                                                                          ttCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaacccc  >  1:2024365/1‑150 (MQ=38)
                                                                                                                                               gggCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACtt                  <  1:178204/129‑1 (MQ=255)
                                                                                                                                                   aaGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAc                                                >  2:1960287/1‑95 (MQ=255)
                                                                                                                                                   aaGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAc                                                <  1:1960287/95‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GGGTATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAG‑‑GCGCTGATCGTGGCGGCTTCGCGCTGATTGTTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCC  >  NC_002947/1637222‑1637507

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