Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,609,8960GA59.6% 28.3 / 55.3 53R224C (CGT→TGT) PP_4929LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base A (17/14);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCA  >  NC_002947/5609749‑5610045
                                                                                                                                                   |                                                                                                                                                      
ctGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGcc                                                                                                                                                      >  2:275653/1‑150 (MQ=255)
    caccaAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATAc                                                                                                                                                  <  2:441212/150‑1 (MQ=255)
     accaAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATAcc                                                                                                                                                 >  2:229844/1‑150 (MQ=255)
      ccaAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCt                                                                                                                                                >  2:216842/1‑150 (MQ=255)
       caAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTg                                                                                                                                               >  1:961918/1‑150 (MQ=255)
                   ccgccAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGAtcgc                                                                                                                                   <  2:961918/150‑1 (MQ=255)
                                tCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCgggg                                                                                                                      >  1:228210/1‑150 (MQ=255)
                                tCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCgggg                                                                                                                      >  1:692143/1‑150 (MQ=255)
                                tCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAACAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCg                                                                                                                                             >  1:179774/1‑127 (MQ=255)
                                tCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAACAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCg                                                                                                                                             <  2:179774/127‑1 (MQ=255)
                                   aGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTaa                                                                                                                   >  2:519672/1‑150 (MQ=255)
                                    gTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAAc                                                                                                                  >  2:293468/1‑150 (MQ=255)
                                       cGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGc                                                                                                               <  1:275653/150‑1 (MQ=255)
                                          ccATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCttt                                                                                                            >  2:915171/1‑150 (MQ=255)
                                                  cccTGGTAGGTCGGGTACTGCACCACCTGCCGGGGTAACCAGGACCAGCCAAGGCCGCGCATCAGCAGTTCGGGCATGGTGTAGAAACTGTCGGCACGCCATACATGCGGGCTGATCGCCTCGACGCCGGGGTAACCGCTTTGTTGTGGc                                                                                                    <  1:312511/149‑1 (MQ=255)
                                                         aGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCa                                                                                             <  1:229844/150‑1 (MQ=255)
                                                                  aCTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGcc                                                                                    >  2:530978/1‑150 (MQ=255)
                                                                    tGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCgg                                                                                  >  2:606350/1‑150 (MQ=255)
                                                                       accacATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTgg                                                                               >  1:64349/1‑150 (MQ=255)
                                                                         cacaTGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCCCCCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCg                                                                             <  1:73895/150‑1 (MQ=255)
                                                                         cacaTGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCg                                                                             >  1:158753/1‑150 (MQ=255)
                                                                         cacaTGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCtt                                                                                                             <  1:596191/118‑1 (MQ=255)
                                                                         cacaTGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCtt                                                                                                             >  2:596191/1‑118 (MQ=255)
                                                                                    gTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCttt                                                                                                            >  1:1046708/1‑108 (MQ=255)
                                                                                    gTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCttt                                                                                                            <  2:1046708/108‑1 (MQ=255)
                                                                                     tAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAAttgttg                                                                 <  2:228210/150‑1 (MQ=255)
                                                                                         cAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCgggg                                                                                                                      <  2:514630/93‑1 (MQ=255)
                                                                                         cAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCgggg                                                                                                                      >  1:514630/1‑93 (MQ=255)
                                                                                         cAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑Cgg                                                            <  1:293468/150‑1 (MQ=255)
                                                                                           ggCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGt                                                          <  1:216842/150‑1 (MQ=255)
                                                                                             cccAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGa                                                        >  1:389646/1‑150 (MQ=255)
                                                                                                 gCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAAtt                                                                     >  2:735362/1‑134 (MQ=255)
                                                                                                 gCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAAtt                                                                     <  1:735362/134‑1 (MQ=255)
                                                                                                  ccAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCgg                                                   >  2:347715/1‑150 (MQ=255)
                                                                                                   cAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGc                                                  <  2:349980/150‑1 (MQ=255)
                                                                                                        ccGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCCGGGTCACCCgccaaca                                              <  1:915171/150‑7 (MQ=255)
                                                                                                        ccGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGATGATCGCCTCGCCGGCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTcg                                             >  1:411265/1‑150 (MQ=255)
                                                                                                          gcgcATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTACCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTcgcc                                           >  2:629391/1‑150 (MQ=255)
                                                                                                          gcgcATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑Cgg                                                            <  1:500365/133‑1 (MQ=255)
                                                                                                          gcgcATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑Cgg                                                            >  2:500365/1‑133 (MQ=255)
                                                                                                           cgcATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTcgccc                                          <  1:492857/150‑1 (MQ=255)
                                                                                                            gcATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑c                                                              >  1:656007/1‑129 (MQ=255)
                                                                                                            gcATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑c                                                              <  2:656007/129‑1 (MQ=255)
                                                                                                               tCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTAGCCCGCCa                                      >  1:47117/1‑150 (MQ=255)
                                                                                                                cagcagTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAg                                     >  2:108540/1‑150 (MQ=255)
                                                                                                                  gcagTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGc                                                                          >  2:17815/1‑112 (MQ=255)
                                                                                                                  gcagTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGc                                                                          <  1:17815/112‑1 (MQ=255)
                                                                                                                   cagTTCGGCCATGGCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAGCgg                                  <  1:519672/150‑1 (MQ=255)
                                                                                                                          gCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAGCGGGt                                >  2:527036/1‑145 (MQ=255)
                                                                                                                          gCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAGCGGGt                                <  1:527036/145‑1 (MQ=255)
                                                                                                                               ggCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACcgc                      >  2:678152/1‑150 (MQ=255)
                                                                                                                                gCGTAGAAACTGTCGGCACACCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACcgcg                     <  1:1008448/150‑1 (MQ=255)
                                                                                                                                gCGTAGAAACTGTCGGCACACAATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACGGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACcgcg                     <  1:321110/150‑1 (MQ=255)
                                                                                                                                                   gCCGTACCTGCTGGCTGATCGCCTCGCCGCCGGGGTAACCGCTGTGTTGTGGCGTGATCACCAACTGGCGGTGGGGGGCCAATTGTTGG‑CGGGTGACCCGGGCCTCGCCCGCAAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCa  >  2:257743/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CTGCCACCAAGGCCGGCGGCCGCCAGTCGCTGTCCAGTTCGACCATGTGCGCCTGGTAGGTCGGGTACTGCACCACATGCCGGGGTAACCAGGCCCAGCCAAGGCCGCGCATCAGCAGTTCGGCCATGGCGTAGAAACTGTCGGCACGCCATACCTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGG‑CGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCA  >  NC_002947/5609749‑5610045

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: ‑