Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,129,303	0	A	G	50.0%	9.5 / 18.5	26	I489V (ATT→GTT)	murJ	putative lipid II flippase
Reads supporting (aligned to +/- strand):  ref base A (2/11);  new base G (13/0);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 2.02e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATGGCGTTTCTGTTGCGTCTGGTGGTGGCGGTACTGGTGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGT  >  NC_000913/1129159‑1129446                                                                                                                                                 |                                                                                                                                                      tGGATGGCGTTTCTGTTGCGTCTGGTGGTGGCGGTACTGGTGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGGTTGTcg                                                                                                                                                 >  2:488920/1‑151 (MQ=255)  ggATGGCGTTTCTGTTGCGTCTGGTGGTGGCGGTACTGGTGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGGGTGCcg                                                                                                                                                 >  2:456942/1‑150 (MQ=255)                               gCGGTACTGGTGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGc                                                                                                              <  2:210168/150‑1 (MQ=255)                                         tGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGGGGGGGGTTCCGCG‑GACTTCCCCGCA‑CC‑GGGGGGACC‑T‑GGGTTTTAAAATTa                                                                                                    >  2:619056/1‑151 (MQ=255)                                           aTGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAg                                                                                                 <  1:448770/151‑1 (MQ=255)                                              tCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGaa                                                                                               <  2:428863/150‑1 (MQ=255)                                                  gCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGGTGGCGGGGGGTGGCGCG‑GGCTTTCCTGGG‑CTGGGGGGGG‑T‑T‑GGGGTTCAAAGTTAAAAAAtttt                                                                                           >  1:458369/1‑150 (MQ=255)                                                  gCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGGGTGCCGCG‑GACTTCGCTGCA‑CC‑GGCGGGG‑C‑T‑GGGGGTCaaaaattaaaaaaattt                                                                                          >  2:507139/1‑137 (MQ=255)                                                              gTATGTTACATATCATTCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTAGTGGTGGGGGGGGGTGTCCCGTTACTTCTCTGCA‑TG‑GGCGGGA‑T‑A‑GGGGTATAAAGTTAAA‑AATTTTTCCGCCggggggg                                                                              >  1:162307/1‑146 (MQ=255)                                                                  gTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGGGTGTCGCG‑TACCTTTCTGCACCG‑GGGGGGG‑C‑G‑GGGGTTAAAAATAAAAAAATTTGGCCCCCGGGCGGGGTa                                                                           >  2:655539/1‑151 (MQ=255)                                                                  gTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGGGTGTCGCG‑GACTTTGCTGCA‑CC‑GGGGGGA‑C‑TGGGGGTTTAAATTTAAAAAATTTGCCCCCCggcgggggga                                                                           >  2:170672/1‑143 (MQ=255)                                                                     aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACaa                                                                       <  2:463275/151‑1 (MQ=255)                                                                     aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACaa                                                                       <  2:153689/151‑1 (MQ=255)                                                                     aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACaa                                                                       <  2:519888/151‑1 (MQ=255)                                                                     aCATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACaa                                                                       <  1:133757/151‑1 (MQ=255)                                                                         atCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGGTGGCGGGGGGTTCCGCG‑GTATTCGCTGCA‑CT‑GGCGGGG‑G‑TGGGGGTTCAAAGTAAAAAAATTTTCCCCccggccgggggaaaaaac                                                                     >  1:186077/1‑140 (MQ=255)                                                                            aTGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGGGGGGGGTTCCGCC‑TGATTTGCTGCACCG‑GGGGGGA‑G‑T‑GGGTTTCAAATTTAAAAAATTTTCCCCCCGGCGGGGGTAaaaagaattt                                                                 >  2:364210/1‑145 (MQ=255)                                                                               ccGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCgg                                                             <  1:236717/151‑1 (MQ=255)                                                                                  gAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGcc                                                           <  2:537962/150‑1 (MQ=255)                                                                                    gTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGGTTTCCCCG‑GACTTCGCTGCC‑CT‑GGCGGGG‑CTT‑GGGGTTCAAATTTAAAAAATTTTCCCCCCGGCGGGGGGAAAAATGCTTTCCGGccg                                                          >  2:199106/1‑149 (MQ=255)                                                                                     tGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgc                                                        <  1:116603/150‑1 (MQ=255)                                                                                         cATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGGTTGCCGCG‑GTCTTTGCTGCA‑AT‑GGGGGGA‑CTT‑GGGGTTCAAAGTTAAAAAATTTTCCCCCCGGCGGGGGGAAAAATGCTTTCCGGCCtcagggg                                                    >  2:546435/1‑145 (MQ=255)                                                                                                 aCCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGGTTTGCGCG‑CGCCTTGCTGCA‑CC‑TGCGGGG‑G‑TGGGGGTTCAAAATTAAAAAATTTGCCCCCCGGCGGGGGTAAAAAAGCTTTCCCGcccccggggggggccgg                                            >  1:486238/1‑136 (MQ=255)                                                                                                        ccTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATc                                    <  1:67010/151‑1 (MQ=255)                                                                                                                          acgtagagggcGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGGA‑C‑T‑GGGGTTCCAAGTTAAAGAAATTGCCCGCCCGGCGGGGGAACAATGGATTCCGGGCTGCAGTGCAGGCCGGAGATAATCTTTCGATCGGAATCtttt                  >  1:49231/11‑151 (MQ=255)                                                                                                                                          ggCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATgggg  >  1:463961/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TGGATGGCGTTTCTGTTGCGTCTGGTGGTGGCGGTACTGGTGATGTCTGGCGTGCTTTTAGGTATGTTACATATCATGCCGGAGTGGTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCG‑TACTTCGCTGCA‑CT‑GGCGGTA‑C‑T‑GGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGT  >  NC_000913/1129159‑1129446

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