Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,739,7640GT50.0% ‑6.5 / 22.5 16A12E (GCG→GAG) ydhBLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base T (5/3);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGTCAGCTCCACATCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCGCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAATTGAGAGATTTAGAGAATATAC  >  NC_000913/1739620‑1739910
                                                                                                                                                |                                                                                                                                                  
gcGGTCAGCTCCACATCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACTCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCTCATc                                                                                                                                                <  3:112319/149‑1 (MQ=255)
           cacaTCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCTCATCAACAACTTc                                                                                                                                       >  5:136534/1‑147 (MQ=255)
             caTCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCGCATCAACAACTTCGAGt                                                                                                                                   <  3:117732/149‑1 (MQ=255)
             caTCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCGCATCAACAACTTCGAGt                                                                                                                                   <  4:1685/149‑1 (MQ=255)
                                 aaGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCTCATCAACAACTTCGAGTGAATATTCTGACCACATAGt                                                                                                               <  6:136534/149‑1 (MQ=255)
                                                  ccaCTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCGCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAAttttt                                                                                              <  6:86856/149‑1 (MQ=255)
                                                                   gCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCTCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCaaa                                                                              >  5:61800/1‑148 (MQ=255)
                                                                                            aaCGCGATGCAGCTCATGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCGCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAat                                                                             >  5:62464/1‑124 (MQ=255)
                                                                                            aaCGCGATGCAGCTCATGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCGCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAat                                                                             <  6:62464/124‑1 (MQ=255)
                                                                                                                            aaCTACCATTACGCGCTACCTCATCAACACCTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAACCCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAAtt                      >  4:83505/1‑147 (MQ=255)
                                                                                                                            aaCTACCATTACGCGCTACCTCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAATTga                    <  6:61800/149‑1 (MQ=255)
                                                                                                                            aaCTACCATTACGCGCTACCGCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAAtt                      >  1:68715/1‑147 (MQ=255)
                                                                                                                                 ccATTACGCGCTACCGCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAATTGAGAGAt               >  2:5237/1‑149 (MQ=255)
                                                                                                                                  cATTACGCGCTACCTCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAATTGAGAGAtt              >  3:99929/1‑149 (MQ=255)
                                                                                                                                            tACCTCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGACTGCAAAATTTTTGAAACCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAATTGAGAGATTTAGAGAata     >  6:140015/1‑148 (MQ=255)
                                                                                                                                              ccGCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAATAGAGAGATTTAGAGAATATAc  <  3:60852/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCGGTCAGCTCCACATCACGGTGCCGTCGTTCAAAGAGCGGCACCGCCAGCCACTCTTCCAGCTGACGCACGGTATAGCTGACCGCAGAAGGAACGCGATGCAGCTCCTGTGCCGCAGCGCTAAAACTACCATTACGCGCTACCGCATCAACAACTTCGAGTGAATATTCTGACCACATAGTCTGCCTGCAAAATTTTTGAAACCAGTCATCAAATATTACCGTTTCACAACACTAATTTCACTCCCTACACTTTGCGGCGGTGTTTAATTGAGAGATTTAGAGAATATAC  >  NC_000913/1739620‑1739910

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