Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,111,747	0	G	T	21.1%	103.7 / 11.6	38	N156K (AAC→AAA)	yiiQ	DUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base G (28/2);  new base T (0/8);  total (28/10)
Fisher's exact test for biased strand distribution p-value = 9.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111615‑4111885                                                                                                                                      |                                                                                                                                            ttttttATTAAAACCTAAAAGCAATCAAGAAGCCGCCAGTATAAGCGGTAAAACAGAGAAGGTCAGCCCCTTTGCGCGGTTGCCCGCGAAAAAATAACGCAGTGCACC‑TTCTGTCTGGGTGTAATAA‑TTTTTTTTTT‑ccc                                                                                                                                    <  2:1610925/137‑1 (MQ=255)     ttATTCAACCCTTAAAGCGTTTCAGATAGCGCCAGTTTAAGCGGTTCAACAGCGAAGGTCACCCCCTTTTCGCTGTTGTCCGCGAAAACATAGCGCAGTGCCCC‑TTCTGTCGCGGTGTAATAACTTTTGTTTTT‑CCccgc                                                                                                                                 <  2:1335859/140‑1 (MQ=255)       tttaaacccTTCAAGCGATTAAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCACCCCTTTTCCGCGGTTGTCCGCAACAACATACCGCGGTGCCCC‑TTCTGTCTGGGTGTAATAACTTTTTTTTTC‑CCccgccg                                                                                                                               <  2:1050303/133‑1 (MQ=255)                  cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGAc                                                                                                                    >  1:877217/1‑140 (MQ=255)                   aaGTGATTGAGATAGCGCCAGCTTAATCGTTAAAACAGCGAGGGTCACCCCTTTTGCGCTTTTGTCCGCAACAACAAACCGCGGTGCCCC‑TTCTGTCTCGGTGTAATAACTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACt                                                                                                                   <  2:1090588/140‑1 (MQ=255)                     gCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGc                                                                                                                 >  1:773902/1‑140 (MQ=255)                       gATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAg                                                                                                               >  2:1048734/1‑140 (MQ=255)                         ttCAGATAGCGCCGCTTTAACCGTTCCAACAGAGAAGGACACCCCTTTTTCCCTTTTGTCCCCAACAAAATACCCCAGTCCCCC‑TTCTCTCTGGGTAAAAAACCGTTTGTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                             <  2:420651/140‑1 (MQ=255)                                  cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                                                                                    >  1:966214/1‑140 (MQ=255)                                  cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                                                                                    >  1:423581/1‑140 (MQ=255)                                  cgcCAACTTAATCGGTTCAACATCGCAGGTCACCGCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGCCTGCAGCTTTTACCGGc                                                                                                    >  1:1505550/1‑140 (MQ=255)                                      aGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                                                                                >  1:770011/1‑140 (MQ=255)                                        cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                              >  2:1131846/1‑140 (MQ=255)                                            aTCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                                                                                          >  2:1549481/1‑140 (MQ=255)                                                 ttCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     >  1:1056970/1‑140 (MQ=255)                                                 ttCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     >  1:526163/1‑140 (MQ=255)                                                 tccAACAGGGAAGGTCCCCCCTTTTTCGCCGTCGCCCCCAACAACAAAACGCGTCGCACC‑TTCTGCCGGGGTGTAATATCGTTTGTTTCC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                     <  2:778095/138‑1 (MQ=255)                                                  tCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                                    >  1:1149902/1‑140 (MQ=255)                                                   aaaCAGCGAAGGCCCCCCCTTTTTCGCCTTTGCCCGCAACAACAAACCGCAGTCCCCC‑TTCTTTCTGGGTTTAATAATTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                   <  2:62447/139‑1 (MQ=255)                                                     aCAGCGAGGGCCACCCCCTTTCCGCCTTTCCCCGCGAAAACA‑AACCCCGTTCCCCTTTCTGTTTGGGTAAAATAATTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                                 <  2:1490736/140‑1 (MQ=255)                                                           aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           >  1:379718/1‑140 (MQ=255)                                                                        tttCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                              >  2:790586/1‑140 (MQ=255)                                                                           cGCCGTTGCCCGCGACAACATAACGCGGTGCCCT‑TTCTGTCTGGGTGTAAT‑ACGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAcc                                                           <  2:210871/140‑1 (MQ=255)                                                                               gTTGCCCGCAAAAACATACCGCAGTCCCCC‑TTCTGTCTGGGTGTAATATTTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       <  2:379718/140‑1 (MQ=255)                                                                                 tgcccGCAACAACAAACCGCAGTCCCCC‑TTCTGTCTGGGTGTAATATTGTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                     <  2:966214/137‑1 (MQ=255)                                                                                 tgcccGCAACAACAAAACGCGGTCCACC‑TTTTGTCTGGTTGTAATAATTTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                     <  2:526163/137‑1 (MQ=255)                                                                                    ccGCGAAAACAAACCGCGTTGCCCC‑TTCGTTCTGGGTATATTAATTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                  <  2:936943/140‑1 (MQ=255)                                                                                     cgcgACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGt                                                 >  1:665148/1‑140 (MQ=255)                                                                                       caaaaacaTACCGCGGTCCCCC‑TTCTGTCTGGGTTTAATACTTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTAt                                               <  2:1056970/136‑1 (MQ=255)                                                                                             caTAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                         >  1:1039621/1‑140 (MQ=255)                                                                                                       tGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  2:466673/1‑140 (MQ=255)                                                                                                       tGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:351308/1‑140 (MQ=255)                                                                                                        gCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  1:1566803/1‑140 (MQ=255)                                                                                                        gCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  1:703222/1‑140 (MQ=255)                                                                                                        gCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  2:136880/1‑140 (MQ=255)                                                                                                        gCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  2:1596607/1‑140 (MQ=255)                                                                                                         cccc‑TTCTGTCGCGTTGTAATACCTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctttc                             <  1:1584218/138‑1 (MQ=255)                                                                                                          ccc‑TTCTGTCTGGTTTTAATATCTTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                            <  2:770011/139‑1 (MQ=255)                                                                                                                  tcctgggTGTAATACCTTTTTTTTCC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     <  1:1131846/135‑1 (MQ=255)                                                                                                                  gTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     >  1:94753/1‑140 (MQ=255)                                                                                                                    ctcGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttttt                   >  2:64994/1‑140 (MQ=255)                                                                                                                     tgggTGTAATATCTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTg                  <  2:694016/138‑1 (MQ=255)                                                                                                                     tcGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTg                  >  1:316691/1‑140 (MQ=255)                                                                                                                       ggTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                >  1:424257/1‑140 (MQ=255)                                                                                                                         tgtAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              >  1:881107/1‑140 (MQ=255)                                                                                                                         tgaaatatCTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  1:1180831/137‑1 (MQ=255)                                                                                                                            aataatTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc           <  2:423581/140‑1 (MQ=255)                                                                                                                             ataagtttttttttC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          <  1:152957/134‑1 (MQ=255)                                                                                                                                 tttttgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  2:238531/138‑1 (MQ=255)                                                                                                                                 cGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      >  2:527581/1‑140 (MQ=255)                                                                                                                                   tttgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:403734/1‑140 (MQ=255)                                                                                                                                   tttgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:1106454/1‑140 (MQ=255)                                                                                                                                   tttgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:1170965/1‑140 (MQ=255)                                                                                                                                   tttgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGAGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTggt    >  2:977794/1‑139 (MQ=255)                                                                                                                                   tttgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGAGTTTTGGTCATTAATGGGGGGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:856902/1‑140 (MQ=255)                                                                                                                                     tgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:253152/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                            TCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACC‑TTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111615‑4111885

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