Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,348,0230GA57.5% 17.1 / 38.0 40A20T (GCA→ACA) nprphosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr)
Reads supporting (aligned to +/- strand):  ref base G (8/9);  new base A (9/14);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTCTGCCAAAGGACGGCAGATTGAAG  >  NC_000913/3347885‑3348158
                                                                                                                                          |                                                                                                                                       
tGCAGAGCACATGGCAGACTACTTCCGCTCGCGCGGTAATAACGTCCAGTCACGCCATCCTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGTGCATGCATGCCCGGCCTGCAAt                                                                                                                                     <  1:135874/143‑1 (MQ=255)
  cAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCGGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCACACTGTTGAAATCACACCCAAGCTGGGCATGCATGCCCGGCCTGCAATGa                                                                                                                                   >  2:631975/1‑143 (MQ=255)
  cAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAACGa                                                                                                                                   >  2:157296/1‑143 (MQ=255)
  cAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACCATGa                                                                                                                                   >  1:19489/1‑143 (MQ=255)
   gggccAACTGGCAGACTCCTTCCGCTCGCGGCGTACCAGTAACCGGTCGCGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAATCTGGGCATGCATGCCCGGCCTGCAATGaa                                                                                                                                  <  2:352548/139‑1 (MQ=255)
   agagCAACCGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACCGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGaa                                                                                                                                  <  2:403474/143‑1 (MQ=255)
             gCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGaa                                                                                                                        <  1:760268/143‑1 (MQ=255)
             gCACACGACGTCCGCTCGCCCGGTAAAGACGTCCAGCCACGCCATCGTACGCTGGAACAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGaa                                                                                                                        <  2:2920/143‑1 (MQ=255)
                 acgactTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTaa                                                                                                                    <  2:358938/140‑1 (MQ=255)
                       gccGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAgg                                                                                                              <  1:35189/142‑1 (MQ=255)
                              gcgcgGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGa                                                                                                       >  1:228118/1‑143 (MQ=255)
                              gcgcgGCAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGa                                                                                                       <  1:91958/143‑1 (MQ=255)
                                gcgGTAAAAACGTCCAGTCACGCCATCGTACGCGAGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACg                                                                                                     <  2:359191/143‑1 (MQ=255)
                                gcgGTAAAAACGTCCAGTCACGCCACCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTCGGCATGCCTGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACg                                                                                                     <  1:630622/143‑1 (MQ=255)
                                    tAAAAACGTCCAGGCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGa                                                                                                 <  1:314879/143‑1 (MQ=255)
                                     aaaagCGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAACCTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTAAAATGAATCTCTGTTAATTAATGCACGTTTTTGACGCTGaa                                                                                                <  2:173987/143‑1 (MQ=255)
                                     aaaaaCGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATCACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGaa                                                                                                <  1:255365/143‑1 (MQ=255)
                                               aGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTAc                                                                                      <  2:369046/143‑1 (MQ=255)
                                                      ccATCGTACGCTGGAAAAACGTAACCCATGGCCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTTAtta                                                                               <  2:535616/143‑3 (MQ=255)
                                                            tACGCTGGAAAAACGTAAACGATGACAGTCAAGCAAACTGATGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGACTTAATGCAGGGTTTTGACGCTGAAGTGTGCTTACGGAATGAGGAggg                                                                         >  1:163484/1‑143 (MQ=255)
                                                                  ggAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGa                                                                   <  2:112818/143‑1 (MQ=255)
                                                                       aaCGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCAGTTTGAATTAATGCAGGGTATTGACGCAGAAGTGCTCTTACGTAATGACGAAGACACCGAGGCTg                                                              >  2:261366/1‑143 (MQ=255)
                                                                          gTAAACCATGACCGTCAAGCAGACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAg                                                           <  2:183895/143‑1 (MQ=255)
                                                                          gTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGCCGCTGACGTGCTCTTACGTAATTACGAAGGCACCGAGGCTAAAg                                                           >  1:518550/1‑143 (MQ=255)
                                                                            aaaCCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGCAGCTGTTCGAATTAATGCAGGGTTTTGACTCTGAAGTGCTCTTACGGCGTGACGAAGGCCGCGAGGCtcgaacc                                                         >  1:208087/1‑137 (MQ=255)
                                                                               ccATGACCTTCACGCAAACTGCGGAAGTCACAAAGAAGCTTGCCATGCATGCCCGGCCTACACTGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTTCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAAc                                                      <  1:57933/143‑1 (MQ=255)
                                                                                    aCCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAAAGAACCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGt                                                 <  1:71731/143‑1 (MQ=255)
                                                                                           aGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTgcg                                          >  2:774380/1‑143 (MQ=255)
                                                                                           aGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTgcg                                          <  2:263506/143‑1 (MQ=255)
                                                                                             cAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTCTGGCgct                                        >  1:491396/1‑143 (MQ=255)
                                                                                                      gAACTCACGAACAAGCGGGGCCTGCATGCCCGGCCTACAATGAGTCTGTGTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACTAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGtt                               <  2:105721/143‑1 (MQ=255)
                                                                                                            acacacagGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTGTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTc                         <  1:464146/143‑1 (MQ=255)
                                                                                                              aaacaaGCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTCTg                       >  1:348735/1‑143 (MQ=255)
                                                                                                                acaaGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGAGTTATTGCGCTGCTGACGTTGGATTCTGcc                     >  2:500927/1‑143 (MQ=255)
                                                                                                                   aGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTCTGCCaaa                  >  2:681385/1‑143 (MQ=255)
                                                                                                                    gccgggCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTAGGATTCTGCCAAAg                 >  2:345240/4‑143 (MQ=255)
                                                                                                                    gCTGGGCATGCATGCCCGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTCATGTTGGATTCTGCCAAAg                 >  2:711626/1‑143 (MQ=255)
                                                                                                                        ggCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTCTGCCAAAGGACg             <  1:530419/143‑1 (MQ=255)
                                                                                                                                  cccGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTCTGCCAAAGGACGGCAGAttgaa   >  1:66103/1‑143 (MQ=255)
                                                                                                                                   ccGGCCTACAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTCTGCCAAAGGACGGCAGAttgaag  >  1:658229/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCATGACCGTCAAGCAAACTGTTGAAATCACAAACAAGCTGGGCATGCATGCCCGGCCTGCAATGAAGCTGTTTGAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTCTGCCAAAGGACGGCAGATTGAAG  >  NC_000913/3347885‑3348158

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