Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	915,190	0	C	G	24.1%	117.3 / 23.6	54	intergenic (‑333/+162)	ybjE/aqpZ	putative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base C (25/16);  new base G (0/13);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 6.59e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCT  >  NC_000913/915052‑915324                                                                                                                                           |                                                                                                                                        gTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCg                                                                                                                                        <  2:863895/140‑1 (MQ=32) gTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCg                                                                                                                                        <  1:956111/140‑1 (MQ=32) gTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCg                                                                                                                                        <  1:697062/140‑1 (MQ=32)  ttCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgg                                                                                                                                       <  1:1360591/140‑1 (MQ=32)   tCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCggg                                                                                                                                      <  1:722417/140‑1 (MQ=32)    cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgggg                                                                                                                                     <  1:982796/140‑1 (MQ=32)         aaGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCa                                                                                                                                <  1:1554755/140‑1 (MQ=32)                 ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                        <  2:1468369/140‑1 (MQ=32)                 ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                        >  1:643013/1‑140 (MQ=32)                 ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                        <  1:231674/140‑1 (MQ=32)                 ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                        <  2:271059/140‑1 (MQ=32)                 ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                        <  1:231290/140‑1 (MQ=32)                 ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcgg                                                                                                                        >  1:204244/1‑140 (MQ=32)                      cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGgcggcggcg                                                                                                                   >  1:1232205/1‑140 (MQ=32)                            tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGgc                                                                                                            <  2:1491177/140‑1 (MQ=32)                            gggATGTGGGCCGAGGGCGCGTATGACAGGGACGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGgc                                                                                                            <  2:1784081/139‑1 (MQ=255)                             ggttgtgGGCTGAGGGCGCGTTTTACAGGTATGTTACCTCGCGCCCGACCCGGTTGCCGTAAGGGATAAGTTGCGGGCACCGCGAAGCGGCGATTTTGTTCGCCAGAGCCGGGGGGTGCAGGGGGCGGCGGCG‑ATTGgcc                                                                                                           <  2:1275513/137‑1 (MQ=11)                                tgtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGgccgcc                                                                                                        >  1:1379551/1‑140 (MQ=32)                                tgtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGgccgcc                                                                                                        >  1:1379545/1‑140 (MQ=32)                                tgtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGgccgcc                                                                                                        >  1:1379399/1‑140 (MQ=32)                                        gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTgcgc                                                                                                >  2:376220/1‑140 (MQ=32)                                        gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTgcgc                                                                                                >  2:48550/1‑140 (MQ=32)                                        gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTgcgc                                                                                                >  2:48588/1‑140 (MQ=32)                                         cgggCGCTTTTTACTGGGATGTTACCACGCGCCCGCCCGGGTAGCCGTACGGTATAAGTGGTGGGCACGGCGCAGCGGCTATTTTGTTCGCATGTGCCGGGGGGTGCAGGGGGGGGCGGCT‑ATTGGCCGCCCCCTgcgcg                                                                                               <  2:1682857/139‑1 (MQ=11)                                                      ctgggATGTTCCCTCGCGCCCGCCCGGTTTGCCGTATGGTTTATGTCGTGGGCACCGCGCTCGGGCTATTTTGTTCGCCAGAGCCGGGGGGGGCTGGGGGGGGCGGCT‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                  <  2:1557687/138‑1 (MQ=2)                                                           ttttacccTCGCGCCCGACCCGGTCGCCTTATGGTATCTGTGTTGGGCACCGCGCTGGGGGTTTTTTTTTGGCCTTCGCCGGGGGGTGCTGGGGGGGGCGGCT‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                             <  2:194679/134‑1 (MQ=2)                                                           ttGTTACCTCGCCCCCCACCGGGTCGCCTTATGGTATACGACGTGGGCCCCGCGCTGCGGCTTTTTTTTTCGCAAGAGCCGGGGGGTGCCGGGGGCGG‑GGCGTTTTGCCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCa                                                                             <  2:485970/139‑1 (MQ=2)                                                              ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGGGGCCata                                                                          >  1:1478858/1‑140 (MQ=18)                                                                     gcgcCCGACCCGGTCGCCGTATGGTATACGTCGCGTGCCCGGCGCTGCTGGTATTTTGTTCGCCACAGCCCGGGGGTGCAGGGGGCGGCGGCG‑TATGCCCGCCCCTCGCGCGCTCCTTGCGACATTGGCAATATGTTGCt                                                                   <  2:1682321/140‑1 (MQ=9)                                                                     gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCt                                                                   >  1:1047990/1‑140 (MQ=255)                                                                      cgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTTGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                  >  2:864497/1‑140 (MQ=255)                                                                      cgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                  >  1:1087600/1‑140 (MQ=255)                                                                      cgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                  >  2:163011/1‑140 (MQ=255)                                                                      cgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCtt                                                                  >  1:1038493/1‑140 (MQ=255)                                                                       gcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTa                                                                 >  1:293239/1‑140 (MQ=255)                                                                       gcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTa                                                                 >  1:456277/1‑140 (MQ=255)                                                                        cccGACCCGGTAGCCGAATGGTTTATGCCGAGGGCACCGCGCAGCGGCTATTTTGTCCGCCAGACCCCGGGGTTGCTGGGGGCGGCGCCT‑TTTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                <  2:1720269/140‑1 (MQ=14)                                                                                               tAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGGGGCAATATGTTGCTTAGCGCATGAAATGAGCGCAACAAga                                         >  1:90161/1‑140 (MQ=255)                                                                                               tAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                         >  1:1083745/1‑140 (MQ=255)                                                                                                  ttCGTGGCCACCGCGCAGCGGCTATTTTGTTCGCCAGAGCCCGGGGTTGCTGGGGGCGGCGGCG‑TTGGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                      <  2:652009/139‑1 (MQ=18)                                                                                                  gtggtgggCACGGCGATGGGGCTTTTTTGTTCGCCTGTGCGGGGGGGTGCTGGGGGCGGGGGCT‑TTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                      <  1:1778534/135‑1 (MQ=18)                                                                                                   tCGAGGGCACCGCGCAGCGGCGTTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATGGACCGCCCCCTGCGCGCTACTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatgat                                     >  2:1484624/1‑140 (MQ=255)                                                                                                    ggtgggCCCGGCGCTGGGGGGTTTTTGTTCGCCTGCGCCGGGGGGTGGTGGGGGGGCCGGCG‑TTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGatgatg                                    <  2:1417299/137‑1 (MQ=18)                                                                                                    ggtgggCACCGGGAGGCGCCTATTTTTTGCCCCTTAGCCGGGGGGTGCTGGGGGGGGGGGCT‑TTGGGCGTCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGatgatg                                    <  2:402339/137‑1 (MQ=16)                                                                                                    ggAGGGCCCGGCCCAGGGCCTATTGTGTTCGCCTTCGCCGGGGGGTGCCGGTGGCGGCGGCG‑AGGGGCCGCCCCGAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGatgatg                                    <  2:1417227/139‑1 (MQ=18)                                                                                                      aGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGaa                                  >  2:566507/1‑140 (MQ=255)                                                                                                       gggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAAAATTTTGCTTAGCTCATGAAAGGAGCGCAACAATATTATTAAt                                 >  1:15879/1‑140 (MQ=18)                                                                                                           aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAAc                             >  2:996810/1‑140 (MQ=255)                                                                                                             cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACaa                           >  1:555838/1‑139 (MQ=255)                                                                                                               cgcAGCGCCGATTTAGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCG‑ATTGGCCGCCCCCATCGAGCTCCTTGCGCCAGGGGCAATAATTTGCTTATCTCATGAAAGTATCGCAACAAGATGATGAAACAACaaat                         >  1:1446110/1‑140 (MQ=18)                                                                                                                                tCGCCAGCGCGGGGGGGTGGGGGGGGGGGGGGGT‑TTGGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaa        <  1:471075/140‑1 (MQ=25)                                                                                                                                     tgcgCCGGGGGTTGCGGGGGGGGGGGGCG‑ATTGGCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc   <  1:1646258/137‑1 (MQ=37)                                                                                                                                     tgcgCCGGGGGGTGCTGGGGGGGGGGGTT‑TTGGGCGCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc   <  1:1464959/137‑1 (MQ=37)                                                                                                                                      gagCCGGGGGGTGCCGGGGGCGGCGGCG‑ATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCt  <  2:1047990/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        GTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCT  >  NC_000913/915052‑915324

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