Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	42.2%	43.4 / 27.9	45	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (19/7);  new base A (0/19);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTG  >  NC_000913/1222447‑1222737                                                                                                                                                           |                                                                                                                                              gCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGa                                                                                                                                              <  1:826785/151‑1 (MQ=255)   aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGaa                                                                                                                                             >  1:494964/1‑150 (MQ=255)   aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGaa                                                                                                                                             >  2:727162/1‑150 (MQ=255)            attacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAgg                                                                                                                                   >  2:654775/1‑151 (MQ=255)              tacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGc                                                                                                                                  <  1:258092/150‑1 (MQ=255)                   attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCa                                                                                                                            <  1:44761/151‑1 (MQ=255)                         gCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTaa                                                                                                                        >  1:624607/1‑149 (MQ=255)                           aacaaTTTCTGCAACTTGTGCGCTGTCATTTTTTGCCT‑TTTCGTTGGCTTCTGCGCAAGTAATATAGTTTTCTGACGGCAAACCTTTTTT‑AA‑TATTGTAGTCCTGCGCCAGGGTCA‑TTGGTG‑CAAAAACAAACGGGCCCGCCAGTAAAGc                                                                                                                     <  2:347847/150‑1 (MQ=255)                              attttCTGCACCCTGTGCGCTGTCATTTTTTGCCA‑TTTCTTTGGTTTCTGCGCAAGAAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCAAGGTCA‑TTGGTG‑CAAAAAAAAACAGGCCCGCCAGTAAAGCaaa                                                                                                                  <  1:393527/148‑1 (MQ=255)                                tttCTGCAACCTGGGCGCTGTCAGTTTTTCCCT‑TTTCGTTGGTTTCTGCGAAAGTAATATAGTTTTCGGACGCCAACCTTTTTTA‑AA‑TATTGTACTCCGGCGCCCAGGTCT‑TTGGTG‑AGAAAAAAAACGGGCCCGCCAGTAAAGCAAAt                                                                                                                 <  1:193634/149‑1 (MQ=255)                                    tGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtc                                                                                                           >  2:300834/1‑151 (MQ=255)                                       aaCCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                        >  2:590885/1‑151 (MQ=255)                                        aCCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                        >  1:715124/1‑150 (MQ=255)                                         ccTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcat                                                                                                      >  2:748884/1‑151 (MQ=255)                                         accggggcgttgtcttttttGCCATTTTTGTGGGCTTTTGCGCAAGTAATATGGTTTTCTGACGCCAACCTTTTTTT‑AT‑TATTGTGGTCCTCCGCCAGGGTCT‑TTGTTG‑GAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                        <  1:763098/135‑1 (MQ=255)                                              gcgcTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTc                                                                                                 >  2:434550/1‑151 (MQ=255)                                                gcTGTCATTTTTTCCCA‑TTTCTTTGTCTTCTGCGCAAATAATATAGGTTTCAGGCGCCAAACTTTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGTTCA‑TTGGTG‑CAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTcc                                                                                                <  2:763089/150‑1 (MQ=255)                                                    tCATTTTTTGCAA‑TTTCTTTGGTTTCTGCGAAAGTAATATAGTTTTCTGACGGCAACCCTTTTTT‑AA‑TTTTGTAGTCCTGC‑CCCAAGTCATTTGGTA‑AGAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                            <  2:448033/150‑1 (MQ=255)                                                        tttttGCCT‑TTTCGTTGTTTTCGGCGAAATAAATATAGTTTTCTGACGGCAACCCGTTTTT‑AA‑TTTTGTATTCCTCCCCCCAGTTTT‑TTGGTG‑CAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAttt                                                                                          <  1:215893/148‑1 (MQ=255)                                                          tttGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAttt                                                                                     >  1:790395/1‑151 (MQ=255)                                                                  tttCGTTGGCTTTGGCGAAATTAATATAGTTTTCTGAC‑GCAAACCTTTTTTAAA‑TATTGTAGTCCTGC‑CCCAAGTTT‑TTTGTGACAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTa                                                                                   <  1:90745/146‑1 (MQ=255)                                                                  tttCGTGGGTTTCGGCGAAAGTAATATGGGTTTCGGACGAC‑AACCTTTTTT‑AATTTTTGTGGTCCGGCCCCAAGTTTT‑TGGGTG‑CAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCa                                                                               <  1:11714/150‑1 (MQ=255)                                                                       ttGGTTTCTGCGCAACAAATATATTTTTCGGACGGCAACCCTTTTTA‑AA‑TATTGTAGTCCTGCCCCCAGGTCT‑TTGGTG‑CAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                          <  2:243291/150‑1 (MQ=255)                                                                       tgggCTTCTGCGAAAATAATATATGTTTCTGCCGACAAACTGTTTAT‑AA‑TATTGTAGCCCTGCGCCAAGTTTT‑TTGGAG‑CAAAAAAAAAAAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAt                                                                           <  1:279817/147‑1 (MQ=255)                                                                         ggtttCTGCGCAAGTAATATAGGTTTCGGACGGAAAACTGTTTTT‑AA‑TTTTGTAGCCCGCCGCCCAGTTAT‑TTGGGG‑AAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                        <  2:396385/147‑1 (MQ=255)                                                                         ggCTTCTGCGAAAGTAATATAGGTTTCTGACGGCAAACCTTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGca                                                                       <  2:849327/151‑1 (MQ=255)                                                                             tCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcaccac                                                                   >  2:458782/1‑151 (MQ=255)                                                                                 cgaaaGAAATATAGTTTTCTGACGGCAACCCTTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCT‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcaccac                                                                   <  2:777231/144‑1 (MQ=255)                                                                                   cAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAc                                                             >  1:377071/1‑151 (MQ=255)                                                                                   aaaGTATGATAGGTTCTTGACGGCAACCCTTTTTT‑AA‑TTTTGAAGTCCTCCACCGAGGTTA‑GTGGGG‑GGAAACAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAc                                                             <  1:861394/150‑1 (MQ=255)                                                                                        aaTATAGTTTTCTGACGACAACCCTTTTTA‑AA‑TATTGTATCCCGGCCCCCGGTTCT‑TTGGTG‑AAAAAAAAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCg                                                           <  1:573694/148‑1 (MQ=255)                                                                                         atttaGTTTTCTGACGGCAACCTTTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGTTCA‑TTGGTG‑CAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCa                                                        <  2:495025/147‑1 (MQ=255)                                                                                           ataGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATc                                                     >  2:526934/1‑151 (MQ=255)                                                                                               ttttCTGCCGCAAAACTTTTTTT‑AA‑TATTGTATGCCGGCCCCAAGTTCT‑TTGGTG‑CAAAAAAAACCAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACg                                                  <  1:654689/149‑1 (MQ=255)                                                                                                tttttggcgcaaaaCCTTTTTT‑AA‑TTTTGTAGCCCTGCCCCCTGTTTT‑TTGGGG‑AAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACg                                                  <  1:534362/138‑1 (MQ=255)                                                                                                tttCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                                >  1:44898/1‑151 (MQ=255)                                                                                                        cgaaaaCCTTTTTTAAA‑TATTGTAGTCCCCCCCCCAGTTCT‑TTGGTG‑AAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                                <  1:434498/141‑1 (MQ=255)                                                                                                           aacccGTTTTT‑AA‑TATTGTATTCCTGCGCCCAGGTCA‑TGGGTG‑CAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcac                                                                      <  2:215899/115‑1 (MQ=255)                                                                                                           aaaCCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcac                                                                      >  1:215864/1‑118 (MQ=255)                                                                                                                       aa‑TATTGTAGCCCGGCCCCCAGTTCA‑TGGGGG‑AAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGa                            <  2:615265/149‑1 (MQ=255)                                                                                                                           atTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCAc                                              >  2:660168/1‑128 (MQ=255)                                                                                                                                aGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGcc                  >  1:260675/1‑151 (MQ=255)                                                                                                                                    cTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc              >  1:352810/1‑151 (MQ=255)                                                                                                                                          ctaaggttttgggTG‑CAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGCACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAtt                     <  1:537454/128‑1 (MQ=255)                                                                                                                                                cA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg  >  1:304175/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                              GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCA‑TTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTT‑AA‑TATTGTAGTCCTGCGCCCAGGTCA‑TTGGTG‑CGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTG  >  NC_000913/1222447‑1222737

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