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

GGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGC  >  NC_000913/1222446‑1222742                                                                                                                                                           |                                                                                                                                                    ggCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCg                                                                                                                                                     >  2:953509/1‑151 (MQ=255)           gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAaacaaaca                                                                                                                                           <  1:669377/151‑1 (MQ=255)            cattacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAg                                                                                                                                          <  2:385900/151‑1 (MQ=255)             attacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAgg                                                                                                                                         >  2:530444/1‑151 (MQ=255)             attacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAgg                                                                                                                                         >  2:484664/1‑151 (MQ=255)                   cattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGcc                                                                                                                                  <  2:928615/151‑1 (MQ=255)                      taccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAg                                                                                                                                <  2:293841/150‑1 (MQ=255)                        ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTa                                                                                                                              >  1:403693/1‑150 (MQ=255)                         cGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTaa                                                                                                                             >  2:789862/1‑150 (MQ=255)                            aacatTTTCTCCAACCTGGGCGCTGTCATTTTTTGCCTTTTCTTTGGTTTCTGCGAAAGTAATATAGTTTTCTGACGGCAAC‑CCTTTTTT‑AAT‑TTTGTAGTCCTGCGCCCAGG‑TCTTTGGTGCAAAAACAAACAGG‑CCCGCCAGTAAAg                                                                                                                           <  1:886762/149‑1 (MQ=255)                            aacaaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGc                                                                                                                          >  1:578411/1‑150 (MQ=255)                            aacaaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGc                                                                                                                          >  1:397720/1‑150 (MQ=255)                               aaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAAt                                                                                                                      >  1:585593/1‑151 (MQ=255)                                       cAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTtca                                                                                                               >  1:1012275/1‑150 (MQ=255)                                        acccTGTGCGCTGTCATTTTTTCCCTTTTCGTTGGCTTCTGCGAAAGTAATATAGTTTTCTGACGGCAAA‑CCTTTTTT‑AAT‑ATTGTAGTCCTGCCCCCAGT‑TCATTGGTGCAAAAACAAACGGG‑CCCGCCAGTAAAGCAAATTTTTtca                                                                                                               <  2:388373/147‑1 (MQ=255)                                         cccTGTGCGTTGTCATTTTTTCCCATTTCGTTGGTTTCTGCGCAAGAAATATAGTTTTCTGACGGCAAC‑CCTTTTTT‑AAT‑ATTGTATTCCTGCGCCGAGG‑TCTTTGGTCCAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                             <  2:540504/149‑1 (MQ=255)                                                     tCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAtt                                                                                                >  2:371264/1‑151 (MQ=255)                                                            ttGCCTTTTCGTGGGTTTCTGCGCAATTATTATAGTTTTATGCCGACAAC‑CCGTTTTT‑ATT‑ATTGTAGTCCTCCCCCAAGT‑TCTTTGTTGCAAAAACAAACAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAttt                                                                                          <  1:361415/150‑1 (MQ=255)                                                              gCCTTTTCGTGGGCTTCTGCGAAAGTAATATAGTTTTCTGAC‑GCAAA‑CCTTTTTTAAAT‑TTTGTAGCCCTGCGCCCAGG‑TCATTGGTGCAAAAAAAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtttt                                                                                         <  2:864831/149‑1 (MQ=255)                                                                         gtgttCGGCGAAAGTAATATAGTTTTCTGACGGCAAC‑CCTTTTTT‑AAT‑ATTGTAGTCCTCCCCCCAGG‑TCATTGGTGCAAAAACAAACAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                             <  2:60228/147‑1 (MQ=255)                                                                            ttcggcgaaataaaTATAGTTTTTTGACGGAAAC‑CCTTTTTT‑AAT‑ATGGTCGTCCGCCGCCATGT‑TCTTTGGTGCAAAAAAACACCGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                              <  1:659752/134‑1 (MQ=255)                                                                            ttCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacc                                                                          >  1:780278/1‑150 (MQ=255)                                                                                 cgcAAGAAATATAGGTTTCTGACGGAAAC‑CCTTTTTA‑AAT‑ATTGAAGTCCTGCGCCCAGG‑TCTTTGGTGCAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCaccacc                                                                       <  2:351260/148‑1 (MQ=255)                                                                                   cAAGTAATATAGGTTTCTGACGGCAAA‑CCTTTTTA‑AAT‑ATTGTAGTCCTGCCCCCAGG‑TCATTGGTGCAAAAAAAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCAccc                                                                      <  2:280648/147‑1 (MQ=255)                                                                                    aaGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAcc                                                                 >  1:26627/1‑151 (MQ=255)                                                                                      taaatttgtgTTTTTGACGGC‑AA‑GCTTTTTT‑AATGAAGGTAGTCCTGCCCCAAGG‑TCTTTGGTGGGAAAAGAACCGGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCg                                                                <  1:861144/141‑1 (MQ=255)                                                                                         atatAGTTTTCGGCCGGCAAC‑CCTTTTTT‑AAT‑TTTGTACCCCTGCCCCCAGT‑TCTTGGGTCCAAAAACAAACAGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCaa                                                            <  2:403840/151‑1 (MQ=255)                                                                                           ataGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAAt                                                           >  2:638163/1‑150 (MQ=255)                                                                                               gttttcgACCGCAAACCCTTTTTT‑AAT‑TTTGTAGTCCGGCCCCGAGT‑TCTTTGGTGCAAAAAAAAACGGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCaa                                                            <  1:164204/140‑1 (MQ=255)                                                                                                               ccttttta‑AAT‑TTTGTAGCCCTGC‑CCCAGGTTCTTGGGTGAAAAAACAAACGGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGccc                                       <  1:13631/148‑1 (MQ=255)                                                                                                                      t‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                  >  1:605760/1‑149 (MQ=255)                                                                                                                        aaT‑ATTTTAGCCCGGCCCCAAGG‑TCTTTGGTCCAAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGa                                 <  2:318916/149‑1 (MQ=255)                                                                                                                        aaT‑ATGGGAGTCCTGCGCCAAGG‑‑CATTGGGGGAAAAAAAAAAAGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                   <  1:735430/147‑1 (MQ=255)                                                                                                                             ttGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAt                           >  2:519685/1‑151 (MQ=255)                                                                                                                                     cTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                   >  1:927266/1‑151 (MQ=255)                                                                                                                                      cgcccccaagttttttgttgCAAAAACAAACCGC‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAtt                          <  1:121936/130‑1 (MQ=255)                                                                                                                                        cccccAGT‑TAATTGGGGCAAAAAAAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAg                 <  1:626410/148‑1 (MQ=255)                                                                                                                                              gg‑TTATTGGGGAAAAAAAAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCa          <  1:837868/151‑1 (MQ=255)                                                                                                                                                 tCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGt        >  1:71151/1‑151 (MQ=255)                                                                                                                                                       gTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGc  >  2:1036597/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                    GGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAA‑CCGTTTTT‑AAT‑ATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGG‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGC  >  NC_000913/1222446‑1222742

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