Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4820GC52.9% ‑0.6 / 16.2 17intergenic (+287/‑169)dapB/carA4‑hydroxy‑tetrahydrodipicolinate reductase/carbamoyl‑phosphate synthetase small subunit
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base C (0/9);  total (5/12)
Fisher's exact test for biased strand distribution p-value = 9.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.24e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGGTGTTATGTTTTTAAATTGTCTATAAGTGCCAAAAATTACATGTTTTGTCTTCTGTTTTTGTTGTTTTAATGTAAATTTT‑GACCATTTG‑GT‑CC‑ACTTTTTTCTGCTCGTTTTTATTTCATGCAAT‑CTTCTTGCTGCGCAAGCGTTTTCCAGAACA‑GGTTAGATGATCTTTTTGT‑CGCTTAAT‑GCCTGTAAAACATGCA‑T‑GAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGGTGTTTTGATTAAGTCAGCGCTATTGGTTCTGGAAGACGGAA  >  NC_000913/29235‑29687
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                             
atttaTTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGGAGTTATGTTTTTAAATTGTCTATAAGTGCCAAAAATTACATGTTTTGTCTTCTGTTTTTGTTGTTTTAATGTAAATTTT‑GACCATTTG‑GT‑CC‑ACTTTTTTCTGCTCGTTTTTATTTCATGCAAT‑CTTCTTGCTGCGCAAGCGTTTTCCAGAACA‑GGTTAGATGATCTTTTTGT‑CGCTTAAT‑GCCTGTAAAACATGCA‑T‑GAGcc                                                                                                                                                                                                             >  1:147794/1‑250 (MQ=255)
                                                                gTCTATAAGTGCCAAAAATTCAATGTTTTGTCTTCTGTTTTTGTTGTTTTAATGTAAATTTT‑GACCATTTG‑GT‑CC‑ACTTTTTTCTGCTCGTTTTTATTTCATGCAAT‑CTTCTTGCTGCGCAAGCGTTTTCCAGAACA‑GGTTAGATGATCTTTTTGT‑CGCTTAAT‑GCCTGTAAACCATGCA‑T‑GAGCCACACAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                                                                            >  2:200369/1‑251 (MQ=255)
                                                                                    aatttttttgttttttgttTTTGTTTTTTTAAT‑TAAATTTTTGACCTTTTG‑GT‑CC‑TTTTTTTTCTGCTTTTTTTTATTTAATGCATT‑TTTCTTGCTGCGCAAGCGTTTTCAAAAAAA‑GGTTAGATGTTTTTTTTGT‑CTCTTAAT‑GCTTGAAAAACTGGCA‑G‑GACCCACAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                                                                        <  1:88611/246‑1 (MQ=255)
                                                                                         ttttttttttttttttttttGTTTTATTGTAATTTTT‑GAC‑TTTTGTGT‑CC‑ATTTTTTTCTGCTTTTTTTTTTTACATGCATT‑TTTTTTGCTGCGAAGGCGTTTCCAAAAAAA‑GGTAAGAT‑ATTTTTTTGTTCG‑TTAATTGCTTGTAAAACA‑GCAAT‑GACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTaa                                                                                                                    <  2:10844/238‑1 (MQ=255)
                                                                                              tttttttttttttttttttTAAT‑TAAATTTTTGAC‑ATTTG‑GTTCC‑TCTTTTTTCTGCTCTTTTTTATT‑CAAGGAATTCTTTTTGCTGCGCAAGCTTTTCCCAGAACG‑GGTTGGATGATTTTTTTGT‑CGCTTAAT‑CCCTGTAAAGAATGCA‑T‑GACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                                                                               <  1:177029/243‑1 (MQ=255)
                                                                                                       tttttttttttATTGTAAATTTT‑GACCATTTG‑TT‑CC‑ATTTTTTTCTGCTCGTTTTTATTTCATGCAAT‑CTTCTTGCGGCGCAAGCTTTTCCCAAAACA‑GGTTAGTTGATTTTTTGGC‑CGCTTAAT‑GCCTGAAAAACATGCA‑T‑GAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                                                                      <  1:402121/243‑1 (MQ=255)
                                                                                                          ttttttttaaatgaaaTTTT‑GACCTTTTT‑GT‑CCCTCTTTTTTTTGCTTTTTTTTATTTCTTGAAAT‑CTTCTTGCTGCGAAAGCGTTTCCCAGACCA‑GTTTAGATGCTTTTTTGGC‑CGTTTGAC‑CCCTTAAAAATTTGCA‑T‑GACCCACAAAAAAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAt                                                                                                    <  1:170567/237‑1 (MQ=255)
                                                                                                                   atttAAATTTT‑GACCATTTG‑GT‑CC‑TTTTTTTTCTGCTGTTTTTTATTCCAGGCAAT‑CTTCTTGCTGCGCAAGCTTTTTCCAGAACA‑GGTGAGATGATTTTTTTGT‑CGTTAAAT‑GCCTGAAAAACATGCA‑T‑GACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                                                                         <  1:220642/248‑1 (MQ=255)
                                                                                                                                  gacttttggctccttttttttCTGCTTGTTTTTTTTTAATGAATT‑TTTTTTGCGCCGCAAGC‑TTTTCCAAAACAAGGTTAGATTATTTTTTTGT‑CGTTTAAT‑GCTTTAAAAACATGAA‑G‑GACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                                                               <  1:72604/235‑1 (MQ=255)
                                                                                                                                  cATTTG‑GC‑CC‑TCTTTTGTCTGCTTGTTTTTATTACATGCAAT‑CTCTTTGCTGCGAAGGCTTTTTCCAGAACA‑GGTTAGATGATCTTTTTGT‑CTCTAAAT‑GCCTTAAAAACAT‑CA‑TTGCGCCAAAAAAAAATAAAAAAATTTCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCTGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                                                                           <  1:268494/251‑1 (MQ=255)
                                                                                                                                   ttttG‑GT‑CC‑TTTTTTTTCGGCTTTTTTTTTTTTCATTCAAT‑CTTCTTGCGGCGCAGGCTTTTTCAAAAACA‑TGTTAGATGTTTTTTTTGT‑CTTTTAAT‑CCTTGAAAACCATGCA‑G‑GGCCCACAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                                                                           <  1:291654/249‑1 (MQ=255)
                                                                                                                                                      tCTGCTCGTTTTTATTTCATGCAAT‑ATTCTTGCTGCGCAAGCGTTTTCCAGAACA‑GGTTAGATGATCTTTTTAT‑CGCTTAAT‑GCCTGTAAAACATGCA‑T‑GAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTAATCATTGTGAATTAATATGCAAATAAagtgagtg                                                           >  2:194483/1‑250 (MQ=255)
                                                                                                                                                                     ttcattaaaT‑CTTCTTGCTGCGCAAGCTTTTTCAAGAACA‑GGTTAAATTATTTTTTTGT‑CGTTAAAT‑GCTTTAAAAACATGCA‑T‑GAGCCAAAAAATAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAgg                                            <  1:31076/243‑1 (MQ=255)
                                                                                                                                                                                                    ttCCAGAACA‑GGTTAGATGATCTTTTTGT‑CGCTTAAT‑ACCAGTACAACATGCA‑T‑GAGCCACAAAATAATATCAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGGTGTTTTGATTAAGTCAGCGCTATTGGTTCt             >  1:32843/1‑251 (MQ=255)
                                                                                                                                                                                                         gAACA‑GGTTAGGTGATCTTTTTGT‑CGCTTAAT‑GCCTGTAAAACATGCA‑T‑GAGCCACAAAATAATATAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGTATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGCGGGTGTTTTGATTAAGTCCGCGCTATTGGTTCTgtaac        >  2:315546/1‑247 (MQ=255)
                                                                                                                                                                                                                gttagagtattttttttT‑CTCTTAAT‑GCCTGAAAACCTTGAA‑T‑GACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGGTGTTTTGATTAAGTCAGCGCTATTGGTTCTGGAAGACGGaa  <  1:237339/240‑1 (MQ=255)
                                                                                                                                                                                                                 ttatatgatCTTTTTGT‑CGTTTAAG‑GCCGATAAAACAGGCA‑T‑GACCCAAAAAAAAATATAAAAAATTCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAAGAGTCTCTGGAGGGTGTTTTGATTAAGTCAGCGCTATTGGTTCTGGAAGACGGaa  <  1:200369/246‑1 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                             
ATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGGTGTTATGTTTTTAAATTGTCTATAAGTGCCAAAAATTACATGTTTTGTCTTCTGTTTTTGTTGTTTTAATGTAAATTTT‑GACCATTTG‑GT‑CC‑ACTTTTTTCTGCTCGTTTTTATTTCATGCAAT‑CTTCTTGCTGCGCAAGCGTTTTCCAGAACA‑GGTTAGATGATCTTTTTGT‑CGCTTAAT‑GCCTGTAAAACATGCA‑T‑GAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGGTGTTTTGATTAAGTCAGCGCTATTGGTTCTGGAAGACGGAA  >  NC_000913/29235‑29687

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