Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,611	0	G	A	41.9%	30.1 / 10.6	31	E143K (GAA→AAA)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (6/12);  new base A (8/5);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCC  >  NC_000913/2151461‑2151753                                                                                                                                                       |                                                                                                                                               aGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAaagaag                                                                                                                                                <  1:624861/151‑1 (MQ=255)  gCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAaagaaga                                                                                                                                               <  2:34688/151‑1 (MQ=255)   cAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAATATATTTCATACTGGATGTATACGTGCGGGGATAAAATTTACGGCGGTTATTCCATTGATCCTTTATTGTTCAGCTACCCAAaagaagaa                                                                                                                                              <  2:598833/151‑1 (MQ=255)       gACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTTAATCTTTTTTTGTGACCCATTCaaaaaaaaaa                                                                                                                                              >  2:185515/1‑147 (MQ=255)             ggTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAg                                                                                                                                     <  1:215862/150‑1 (MQ=255)               tATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTTATTCTTTTTTTGTGACCTCTCCAAaaaaaaaaga                                                                                                                                            >  1:418263/1‑140 (MQ=255)                     tACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGc                                                                                                                             <  2:432273/150‑1 (MQ=255)                      aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCa                                                                                                                            <  2:525766/150‑1 (MQ=255)                        aaaCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCaaaa                                                                                                                         <  1:188091/151‑1 (MQ=255)                        aaaCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTGATCCTTTTTTTGTCACCTATCCAAAAAAAAAAGGCGACGAGCTAAGAGCaaa                                                                                                                          >  1:492638/1‑150 (MQ=255)                              aGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATTTACCATTCATCGTTTATTGGGTTCATATCAACAAGAAAAAGCCGGAGAGGTAAGATGAAAATTa                                                                                                                      >  2:494865/1‑148 (MQ=255)                                  tatGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTGATTCTTTTTTTGTTACCTATCCAAAAGAAAAAACCGGCGAGCTAAGAGGAAAATTAGTTTGCt                                                                                                               >  2:144810/1‑151 (MQ=255)                                   atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCt                                                                                                               <  2:424316/150‑1 (MQ=255)                                    tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCAATTATCCTTTTTTTTTCACCTATCCCAAAAAAGAAGCCGACGAGCTAAGAACAAAATTTGttt                                                                                                                  >  2:151072/1‑145 (MQ=255)                                    tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATATATTATTTGTTCTTTTTTTGTCAACTAATCAAAAGAAAAAACCGACCAGCTAAGAGTAAAAATAGtt                                                                                                                   >  2:611663/1‑145 (MQ=255)                                     ggCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAAt                                                                                                            <  2:391218/151‑1 (MQ=255)                                         aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTTATCCTTTTTTTGTCACCCATCCCAAAAAAAAAAACGACGAGCCAAAAGCAAAAATAATTTGCTAAtaaa                                                                                                         >  1:465513/1‑147 (MQ=255)                                             gTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGtt                                                                                                    <  1:108706/151‑1 (MQ=255)                                             gTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACAATTATCCTTTTTTTGTCACCCATCCCAAAAAAAAAACCGACGAGCTAAGAGGAAAAATAATTTGCTAATCAACCggt                                                                                                    >  1:482161/1‑149 (MQ=255)                                                  aaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatt                                                                                               <  1:460314/151‑1 (MQ=255)                                                                    tttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCATTGATCCTTTATTTGTTACCTATCCCAAAAAAAAAGCCGGAGAGCTAAGAGCAAAAATAGTTTGCTAAACAACCGCTTTATTTTTTTTTTtaa                                                                                    >  1:502410/1‑144 (MQ=255)                                                                     ttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTTATCCTTTTTTTGTTACCTCTTCCAAAAAAGAAACCGCCGCGGTAAGAGGAAAAATAAATCGCTAATTAACCGCTTTATTTTTTTTTTTTAATGGCggg                                                                            >  1:644059/1‑150 (MQ=255)                                                                            tGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCATTTATCCTTTTTTTGTCACCCATCCCAAAAAAAAAAACGACGAGCTAAAAGCAAAAATAGTTTGCTAATCAACCGGTTTATTTTTTTTTTTAAATAGCgggt                                                                           >  1:377902/1‑143 (MQ=255)                                                                                     aCGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCa                                                            <  1:573310/151‑1 (MQ=255)                                                                                       gcgcgGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAc                                                          <  1:714769/151‑1 (MQ=255)                                                                                        cgcgGGGATAAAATTTACGGCGGTTATAACCTTTTTCCTTTTTTTGTTACCCATTCAAAAAAAGAAAACGAAGAGCTAAAAGCAAAAATAGTTTGCTAATCAGCCGCTATATTTTTTGTTTTTAATGGGGGGTGGAttgttttttttt                                                            >  1:520763/1‑141 (MQ=255)                                                                                                   atttgaCGGCGTTTATACTATTGATCTTTTATTGTTCCCCTATCAAAAAGAAGAACCCGAGGAGCAAAGAGCAAAATTATTTCGCAAATCAGCCGTTATATTTTTTGTTAAAACGGGGGGTGGAATTTTTGTCTTCAACCGCCATGAGaaa                                              <  2:321239/146‑1 (MQ=255)                                                                                                          ggggtttATACCTTTGTCCCTTTTTTGGCCCCCTCTCAAAAAAAAAAGGCCGACGGGCAAAGAAAAAATTTAGTCCCTTACCCACCCGTTATTTTTTTTGTTAAAATGGGCGGTTGATTTGTTTTCTTCAACCGCCATGAGAAactca                                          <  2:66870/143‑1 (MQ=255)                                                                                                           gggtttATACCATTGATCTTTTTTTATTCCCTTTTCAAAAAGAAAAAGCCGAGGTGAAAAGAGAAAAATTATTCCGCTACTCACCGTTTATTTTTATTTTTAAAACGGGCGGTTTAATTTTTTTTTTCAACCGCCATGAGAAACTCActct                                      <  1:611952/147‑1 (MQ=255)                                                                                                               ttATCCCTTTGATCTTTTATTGTTCCCCTATCAAAAAGAAGAAGCCGCCGAACTAAGAGAAAAATTATTTCGCTAATCAGCCTTTATATTTATTGTTAAACCTGGCGTTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTg                                     <  2:438556/148‑1 (MQ=255)                                                                                                                      aTTGTTCCTTTATGGGTCACCTATCAAAAAAAAAAAGCCGAAGAGTTAAGAGCAAATTTAGTGCGCTAATCACCCTTTTTTTTTTTTGTTTAAACGGGCGTTTTAGTTGTTGTCCTCACCCGCCATGAGAAACTCACTCTGCCca                                 <  2:15598/145‑1 (MQ=255)                                                                                                                       ttGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGcc                          >  2:504159/1‑151 (MQ=255)                                                                                                                                  tggttccccTACCAAAAAAAAAAAGCCGACGGTCTAAGAACAAAATTATTTCGCTAACCAGCCTTTATTTTTATTTTTAAACCGGGGGGTTGATTTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAAc                 <  2:216207/142‑1 (MQ=255)                                                                                                                                               cAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTTTCATCATcc  >  2:269887/1‑151 (MQ=255)                                                                                                                                               cAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTCAACCGCCATGCGAAACTCACTCTCCCCACCGCGCCCGTGCCAAAGCTTTATCATTATcc  >  1:652494/1‑151 (MQ=255)                                                                                                                                                aaaaGAAAAAGCCGACGAGTAAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTTTTAAAACGGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat       <  1:151097/145‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               AGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCC  >  NC_000913/2151461‑2151753

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