Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,589	0	A	T	40.6%	32.6 / 20.3	33	L135F (TTA→TTT)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (6/13);  new base T (12/1);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 9.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

CATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCC  >  NC_000913/2151442‑2151735                                                                                                                                                      |                                                                                                                                                   cATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑TT‑AACAATGAATCtttttt                                                                                                                                                  >  1:237710/1‑150 (MQ=255) cATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTg                                                                                                                                                 <  1:519938/151‑1 (MQ=255)       tATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTTATCCTTTAttttttacc                                                                                                                                           >  2:298665/1‑144 (MQ=255)           ggCGGAATCAGCAATGACCCATATATTGGTACTAACGTAGAATAGGGCCAGAAGAGGGAAATCAAGAAAGAAGATATTTCAGACTGGATGGCTACTCGCGGGGCTAAAATTTACGGCGGTT‑AT‑CCCATTGAGCCTTTTTTGTTAAccg                                                                                                                                          <  1:154077/148‑2 (MQ=255)                             ccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑TT‑AACATTTATTCTTTTTTTGTCACCTATCCaaaaaaaaaa                                                                                                                            >  2:521416/1‑144 (MQ=255)                                         aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCa                                                                                                          <  1:421205/150‑1 (MQ=255)                                           aaaCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑TT‑AACCTTTATCCTTTTTTTGTCTCCCCTCCCAAAAAAAAAACCGACGAGGTAAAAGaaaa                                                                                                        >  2:707570/1‑150 (MQ=255)                                                      atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑AACATTTATTCTTTTTTTGTCACCCATTCAAAAAAAAAAAACGACGAGCCAAGAGCaaaa                                                                                                       >  2:124815/1‑140 (MQ=255)                                                       tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑AACATTTATCCTTTTTTTGTTACCCTTCCCAAAAAAGAAAACGACGAGGCAAAAACAAAAATAAtt                                                                                                 >  2:420238/1‑145 (MQ=255)                                                           cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGCT‑AT‑AACAATGTTTCTTTTTTTGTTAACCATCCAAAAAAAGAAAACGAAGAGCTTAAAGAAAAAATAGTTTCCTaa                                                                                           >  1:208322/1‑147 (MQ=255)                                                             ggAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑AACATTTATTCTTTTTTTGTTACCTCTTCCAAAAAAGAAGCCGACGAGCTAAGAGCAAAAATAGTTTGCtcatca                                                                                        >  2:682265/1‑148 (MQ=255)                                                               aGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑TT‑ACCATTTATTCTTTTTTTTTTACCCATCCCAAAAAAAAAAACGAAGAGCTAAAAGGAAAAATAATTTGCTTAttaacc                                                                                     >  2:221936/1‑144 (MQ=255)                                                                gTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑TT‑TCCATTTTACTCTTTTTTGTCTCCTATCCAAAAAAAGAAACCCAAGAGGTAAAAGGAAAAATAAATTTCTAATCAGcc                                                                                     >  2:480673/1‑148 (MQ=255)                                                                gTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGTTCCTTTTTTGGGTACCTATCCCAAAAAAAAAGACAACGAGCTAAGAGAAAAATTAATTTTCTAATCaaacgt                                                                                   >  1:649615/1‑145 (MQ=255)                                                                 tGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑TTAACAATTATTCCTTTTTTGGTAACC‑ATTCAAAAAAAAAAGACGACGAGCTAAGAGCAAAAATAATTTGCTAATCAGcccct                                                                                  >  1:100628/1‑147 (MQ=255)                                                                    aaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTAAA‑ACAATTTTTCCTTT‑TTTGTCACCCATCCAAAAAAAAAAACCGACGAGCTAAAAGAAAAATTAATTTGCTAAACAACCGTTatat                                                                              >  1:66237/1‑151 (MQ=255)                                                                     aaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTTTTTACCCATCCCAAAAAAAAAGGCGACGAGCTAAAAGGAAAAATAATTTGCTAATCAGCCGTTAtttt                                                                             >  2:253588/1‑151 (MQ=255)                                                                         cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttatt                                                                         <  2:386517/151‑1 (MQ=255)                                                                         cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttat                                                                          <  1:502404/150‑1 (MQ=255)                                                                         cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttat                                                                          <  2:512379/150‑1 (MQ=255)                                                                              aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTa                                                                    <  2:494500/151‑1 (MQ=255)                                                                                gaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTaaa                                                                  <  1:126139/151‑1 (MQ=255)                                                                                    atatTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTgg                                                              <  1:343484/151‑1 (MQ=255)                                                                                      atTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCg                                                            <  1:560141/151‑1 (MQ=255)                                                                                         tCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGt                                                          <  2:323184/150‑1 (MQ=255)                                                                                           aGACCGAATGTATACGCGCGGGGATAAAAGTTACGGCGGTT‑AT‑ACAATTGATCCATTATGGGTCCCCTAGCCAAAAGAAGAAGCCGAAGAGCTAAGAGCAAAATTAGTGCGCTAATCAGCCGTTATATTTATTGTCTAAACAGACgga                                                          <  1:138227/148‑2 (MQ=255)                                                                                            gACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ATCAATTGTTTTTTTTTGGTTACCTATCCCAAAAAAAAAGAAGACCAGGTAAAAGCAAAAATAGTTTGCTCAATTCCCGTTTTATTTATTTTTTataaagacggggtt                                                       >  1:561424/1‑137 (MQ=255)                                                                                            gACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGAT‑TT‑AACACTTAATCTTTTTTTGTGACCTATCCCAAAAAAAAAAACGGCCAGCTAAAAGGAAAAATAGTTTGCTAATAAACCGCTATTTTTTTTTTTTTAAATAGCgg                                                           >  2:256623/1‑146 (MQ=255)                                                                                                    gTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTc                                              <  2:206981/151‑1 (MQ=255)                                                                                                                              gcgtTT‑AT‑ACCATTGATCTTTTATGGGTCCCCTATCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTCCGCTAATCAGCCGTTATATTTATTGTTAAAACGGGCGTTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCActct                    <  2:420440/151‑1 (MQ=255)                                                                                                                                ggTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCActct                    >  2:432202/1‑149 (MQ=255)                                                                                                                                               aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCGCCATTAGAAAATCACTCTGCCCACAGCGCCAg      >  1:165183/1‑150 (MQ=255)                                                                                                                                               aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCCCTCTGCCCCCAGCGCCAGt     >  1:133491/1‑151 (MQ=255)                                                                                                                                                  cTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGcc  >  1:141964/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CATGTTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTT‑AT‑ACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCC  >  NC_000913/2151442‑2151735

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