Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,111,747	0	G	T	44.7%	48.3 / 18.3	38	N156K (AAC→AAA)	yiiQ	DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (15/6);  new base T (0/17);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 4.44e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTA  >  NC_000913/4111602‑4111889                                                                                                                                                   |                                                                                                                                                gctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                                                                                            <  1:395335/151‑1 (MQ=255)    ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑cc                                                                                                                                         <  2:506770/150‑1 (MQ=255)       ctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCccgc                                                                                                                                     >  2:257026/1‑151 (MQ=255)              tCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTtagt                                                                                                                              >  1:6796/1‑151 (MQ=255)                        aaccctcaaaGCGATTCAGATGGCGCCAGCTTAATCGTTTCAACAGCGAAGGCCGCCCCCTTTT‑CGCCTTTGCCCGCAACAACATAACGCGGTGCCCTTTCTGTCTCGGTGAAATAACGTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCa                                                                                                                    <  1:471248/143‑1 (MQ=255)                        aacccTGCAAGCGATTCATATGGCGCCAGATTAAGCGGTCCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGAAAACATAACGCAGTGCACCTGCTGCCTCGGTGTAATAACTTTTTTTTTC‑CCCCGCCGTTAGTAGc                                                                                                                           <  2:162290/141‑1 (MQ=255)                                aggcaaTTCAGATGCCCCCATCTTAATGGGTACAACAGCGAGGGCCACCCCTTTTC‑CGCCGTTTTCCGCGAAAACATACCGCAGTGCCCTTTCTTTCTGGTTTTAAAACCTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGttt                                                                                                                <  1:465371/142‑1 (MQ=255)                                     aTTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                       >  2:101606/1‑151 (MQ=255)                                        cAGATAGCCCCATCTTAATCGTTCCAACAGCGAAGTTCAGCCCCTTTT‑CGCTGTTGTCCGCGACAACATACCGCATTGCCCTTTTTGCCGCGGTGTAAGACATTTTTTTTTT‑CTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                                                                                    <  2:635646/151‑1 (MQ=255)                                               cgcCATCTTAATCGTTTCAACAGCGAGGGTCACCCCCTTTT‑CGCTGTTGCCCGCGAAAACATACCGCAGTGCCCCTTCTGTCTGGGTGAAATATCTTTTGTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                                  <  2:639928/146‑1 (MQ=255)                                                gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                                                                                            >  1:489353/1‑151 (MQ=255)                                                  cAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTc                                                                                          >  1:642692/1‑151 (MQ=255)                                                      taaaTCGGTTCAACAGCGAAGGCCACCCCCTTTC‑CGCCGTTGCCCGCAACAACAAACCGCAGTCCCCCTTCTGTCTGGTTGTAATAACGTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                       <  2:633407/148‑1 (MQ=255)                                                        aaTCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                                                                                    >  2:571347/1‑151 (MQ=255)                                                          tcgttcaaaCAGCGAAGGCCACCCCTTTTT‑CGCCGTTGTCCGCGAAAACATACCGCAGTGCCCCTTTTGTCTCGGTGTAATACCGTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATg                                                                                   <  2:285641/143‑1 (MQ=255)                                                              ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTg                                                                              >  2:162701/1‑151 (MQ=255)                                                                   cAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtgt                                                                         >  1:448138/1‑151 (MQ=255)                                                                         aGGCCACCCCTTTTC‑CGCCGTTGTCCGCAAAAACATACCGCATTCCCCCTTCTGTCTCGGTTTATTATTTTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGa                                                                    <  1:109924/150‑1 (MQ=255)                                                                          ggccAGCCCCTTTT‑CGCCTTTGCCCGCGACAACATACCGCAGTGCACCTTCTGTCGCGGTTTAATAAGGTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGAt                                                                   <  1:407478/147‑1 (MQ=255)                                                                          ggTCACCCCCTTTT‑CGCCGTTTCCCGCAACAACATAACGCAGTGCCCCTTCTGTCTCGGTTTAAT‑ACGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtgt                                                                         <  2:102637/144‑1 (MQ=255)                                                                           ggcccccccTTTTCCGCCGTTGCCCGCGCAAACATACCCCGGTCCCCTTTCTTTCTCGGTTTAATACCTTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                                  <  2:199957/146‑1 (MQ=255)                                                                                 ccttttc‑CGCCTTTGCCCGCAAAAACAAAACGCGGTGCCCTTTCTGTCTCGGTGTAATAATGTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                           <  1:433966/148‑1 (MQ=255)                                                                                      tt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                      >  2:190483/1‑151 (MQ=255)                                                                                                  ccGCAACAACATACCGCTTTCCCCCTTCTTTCTGGGTGTAATATCTTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCt                                               <  1:595265/147‑1 (MQ=255)                                                                                                       acaacaTACCGCAGTGCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttc                                       <  2:312112/150‑1 (MQ=255)                                                                                                       acaacaTAACGCAGTGCACCTTCTGTCTGGGTGTAATAATTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCt                                               <  2:528930/142‑1 (MQ=255)                                                                                                         aacaTACCGCTGTGCCCTTTCTGTCTCGGTTTAATATGTTTTGTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgct                                    <  1:380390/151‑1 (MQ=255)                                                                                                         aacaTACCGCAGTGCCCCTTCTGTCTCGTTTTAATATCTTTTGTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcg                                      <  2:232483/149‑1 (MQ=255)                                                                                                          acaTAACGCAGTGCCCCTTCTGTCTCGGTTTAATATTGTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                   <  1:302874/151‑1 (MQ=255)                                                                                                            aTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctttc                                 >  2:52871/1‑151 (MQ=255)                                                                                                              accGCAGTGCACCTTCTGTCTCGGTGTAAAAATTTTTGTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                                <  2:663586/148‑1 (MQ=255)                                                                                                                  cAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCt                           >  1:307360/1‑151 (MQ=255)                                                                                                                       cACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTg                      >  2:534763/1‑151 (MQ=255)                                                                                                                           ttctttctGGGTGTAATAACTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                            <  2:442117/141‑1 (MQ=255)                                                                                                                            tCTGTCTCGGTGTACTAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTg                 >  1:111783/1‑151 (MQ=255)                                                                                                                               gTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc              >  2:210728/1‑151 (MQ=255)                                                                                                                                 ctcGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcccc            >  1:177076/1‑151 (MQ=255)                                                                                                                                            accGTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTa  <  1:558244/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTA  >  NC_000913/4111602‑4111889

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