Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,352,7540AC50.0% 3.0 / 15.4 18V91G (GTG→GGG) arcBsensory histidine kinase ArcB
Reads supporting (aligned to +/- strand):  ref base A (7/2);  new base C (0/9);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGA  >  NC_000913/3352693‑3352817
                                                             |                                                               
attatCTTTTAACTGAACGTTGAGGCTCAAATCGCTCTCGCGCATCTCCTCCAGTTTTTGCACCa                                                              >  1:658350/1‑65 (MQ=255)
attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCa                                                              >  1:2424366/1‑65 (MQ=255)
attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCa                                                              >  1:2112099/1‑65 (MQ=255)
attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCa                                                              >  1:3457417/1‑65 (MQ=255)
     cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt                                                         >  1:3614306/1‑65 (MQ=255)
     cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt                                                         >  1:2049604/1‑65 (MQ=255)
     cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt                                                         >  1:2081439/1‑65 (MQ=255)
     cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt                                                         >  1:2281488/1‑65 (MQ=255)
      ttttAGCTGAGCGTGGAGGCTCAAGCCGCGCCCGCGTATCTCCTCCAGTTTTTGCACCAGCCGTg                                                        <  1:3179562/65‑1 (MQ=255)
      ttttAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTg                                                        <  1:1928123/65‑1 (MQ=255)
      ttttAACAGAACTTTGAGGCTCAAACCGCGCTCGCGCATCCCCTCCAGTTTTTGCCCCAGCCGTg                                                        <  1:2509312/65‑1 (MQ=255)
          aaCTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAg                                                    >  1:303456/1‑65 (MQ=255)
            cTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGAc                                                  <  1:3442555/65‑1 (MQ=255)
                 cGTTGGGGCTCAAACCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt                                             <  1:1954647/65‑1 (MQ=255)
                 cGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt                                             <  1:2740609/65‑1 (MQ=255)
                 cGTTGAGGCTCAAATCGCGCTCGCGCATCCCCTCCATTTTTTCCCCCAGCCGTGACAGACGTTGt                                             <  1:2900193/65‑1 (MQ=255)
                 cGTTGAGGCTCAAATCGCCCTCGCGCTTCTCCTCCATTTTTTGCACCAGCCGTGACAGACGTTGt                                             <  1:3796255/65‑1 (MQ=255)
                 cGTTGAGGCTAAAACCGCGCTCGCCCATCTCCGCCAGTTTTGGCCCCAGCCGTGACAGACGTTGt                                             <  1:2609433/65‑1 (MQ=255)
                 cGTTGAGCCACAAATCGCGCTCGCGCATCCCCTCCATTTTTTGCCCCAGCCGTGACAGACGTTGt                                             <  1:1036062/65‑1 (MQ=255)
                      aGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGa                                        <  1:1348979/65‑1 (MQ=255)
                        gCTCAAATCGCGCTCGCGCATCTCCTCCATTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct                                      <  1:2501809/65‑1 (MQ=255)
                        gCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct                                      <  1:3821145/65‑1 (MQ=255)
                        gCTAAAATTGCCCTCTCGCATCTCCCCCATTTTTTCCACCAGCCGTGACAGACGTTGTCGTGAct                                      <  1:912267/65‑1 (MQ=255)
                        gCTAAAACCGCGCTGGCGCATCTCCTCAAGTTTTTCCACCAGCCGTGACAGACGTTGTCGTGAct                                      <  1:1568914/65‑1 (MQ=255)
                        ccTCAAACCGCGCTCGCGCATCCCCCCCACTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGAct                                      <  1:3643829/64‑1 (MQ=255)
                         cTCAAACCGCGCTCGCCCATCCCCTCCATTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGActc                                     <  1:3224019/65‑1 (MQ=255)
                          tCAATTCGCGCCCGCGCATCTCCTCCAGTTTTTCCACCAGCCGTGACAGACGTTGTCGTGActcc                                    <  1:3178003/65‑1 (MQ=255)
                          tCAACTCGCCCTCGCGCACCCCCCCCTGTTTTGCCCCCAGCCGTGACAGACGTTGTCGTGActcc                                    <  1:1921330/65‑1 (MQ=255)
                           cAAAGCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGActcct                                   <  1:843137/65‑1 (MQ=255)
                              aTCGCACTCGCGCCTCTCCCCCATTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCa                                <  1:2661051/65‑1 (MQ=255)
                                       ccgcTTCTCCCCCAGTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTcga                       <  1:235618/64‑1 (MQ=255)
                                             ctcctcCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGacca                 >  1:2294676/1‑65 (MQ=255)
                                               cttccAGTTTTTGGCCCAGCCGTGACAGACGTTGTCGTGCCTCCTCCAGTTGCTCGACGaccacc               <  1:1748005/63‑1 (MQ=255)
                                               cctccGTTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGcccacc               <  1:2332601/65‑1 (MQ=255)
                                                  ccAGTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGAt            <  1:2219321/65‑1 (MQ=255)
                                                            cACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGa  >  1:900303/1‑65 (MQ=255)
                                                             |                                                               
ATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGA  >  NC_000913/3352693‑3352817

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: