Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	827,499	0	T	G	44.4%	51.7 / 13.7	48	T495P (ACC→CCC)	ybhF	putative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (6/19);  new base G (20/0);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 7.46e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTC  >  NC_000913/827352‑827632                                                                                                                                                    |                                                                                                                                       ccGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAg                                                                                                                                     <  2:597156/151‑1 (MQ=255)                  aTTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAg                                                                                                                    <  1:1012848/151‑1 (MQ=255)                    tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTc                                                                                                                  <  1:907618/151‑1 (MQ=255)                                gTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAg                                                                                                                                     <  1:480336/120‑1 (MQ=255)                                gTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCCGGCGTCGGTttgtcccga                                                                                                      >  2:317726/1‑147 (MQ=255)                                  acacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGGGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCCACGCCGGGAGTggggttttccggaaa                                                                                                    >  1:169091/1‑139 (MQ=255)                                   cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGGCGAGaaaa                                                                                                   >  1:78974/1‑150 (MQ=255)                                   cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGCGaaaa                                                                                                   >  1:357343/1‑150 (MQ=255)                                                      cAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCa                                                                                <  2:1148356/151‑1 (MQ=255)                                                        atatTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTGAACGGCGGGAGGCGGGTTGTCGAGACAAAGAAAATCCGGtgcatgcac                                                                               >  2:740697/1‑149 (MQ=255)                                                        atatTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCCCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCAt                                                                               >  1:1099745/1‑150 (MQ=255)                                                          tttCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGTGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                                                                            <  2:1120046/150‑1 (MQ=255)                                                                     tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCCGAAGTCGGTTCGGCGGGAAAAAGAATGGCCGGTTCATGCATCAGCGAAAA‑gggc                                                                 >  1:242523/1‑148 (MQ=255)                                                                      ccATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTTATGCATCAGCGGAAA‑Ggcccg                                                               >  1:310708/1‑149 (MQ=255)                                                                         tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑Ggccagcgc                                                            <  1:343253/151‑1 (MQ=255)                                                                         tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑Ggccagcgc                                                            <  2:824295/151‑1 (MQ=255)                                                                         tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATGCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑Ggccagcgc                                                            <  2:564184/151‑1 (MQ=255)                                                                         tAAAGAGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTAAACGCCCGGAGTCCGTTTGTCGAGGAAAAGAATGTTCGGGTCTTGGATTAGCGAAAA‑Agcccgcg                                                             >  1:575037/1‑150 (MQ=255)                                                                                     tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGGGCCAGCCGCTGttt                                                <  1:932720/151‑1 (MQ=255)                                                                                      gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTCCTGCATCAGCGAACA‑GGCCAGCGCCCGCCCCCGTTTa                                               >  1:244453/1‑151 (MQ=255)                                                                                      gACCATCACCGTGACGCCTTCCTCTACCATGCTGTTGATTTGCAGCCAAAATTCACGGCGGGCGGGGGGGTTCATGCCCGGACTCGCGTCGTTTGGGAACACCAAGTGCGGTTCAATCATCAGCGCACA‑AGCCAGCGCCAGcccaggtttt                                               >  2:687777/1‑143 (MQ=255)                                                                                          aacACCGTGACGCCTGTCTCTACCATGCGGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATACATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAAc                                           <  2:61697/149‑1 (MQ=255)                                                                                             aCCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTGCGGGTCATGCCTCAGCGAAAA‑GGCCAGCGCCCGCCCCCGTTTTAAACCTa                                        >  2:926446/1‑151 (MQ=255)                                                                                              ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGTGGGGGGCCACGCCCGGAGGCGGTTCGTCGGGAAACAGAAAGTCCGGTTCATGGATCAGCGAACA‑AGCCCGCGCCAACCGCCGTCTTTAACCtta                                       >  2:253589/1‑149 (MQ=255)                                                                                                   aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                   <  1:838906/150‑1 (MQ=255)                                                                                                   aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                   <  1:764275/150‑1 (MQ=255)                                                                                                      ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGTGGAGTCGTCCAGAAACAGAAAGTTCGGGTTATGCATCAACGAACA‑GGCCAACGGCCGACGCTGTTTAAAACCTAATGTCCGtt                               >  2:818616/1‑151 (MQ=255)                                                                                                      ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGGTCGTCGGGAAAAAGAATGTCCGGGTCATGCCTCAGCGAACA‑GGCCCGCGCCCGCCGCCGTTTTAAAACTAATGGGAggt                               >  1:371503/1‑149 (MQ=255)                                                                                                      ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTCATTCATCAACGAAAA‑AGCCAGCGCCCGGCCCTGTTTTAAAACCAATTGCAggt                               >  1:121796/1‑149 (MQ=255)                                                                                                        tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCTGGAGTCGGTTCGTCGAGGAAAAGAATGTCCCGGTCCTGGCTCCGCGAAAA‑GGCCCGCGCCCGCCGCCGTTTAAAACCCAAACGCAggtta                             >  2:75150/1‑147 (MQ=255)                                                                                                        tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTTGGGTTTTCGAGAAAAAGAATGTCCGGGTCATGTCTCCGCGAACA‑AGCCAGCGCCAGCCGCTGTTTTAAAACTAATGGccggt                               >  1:540062/1‑145 (MQ=255)                                                                                                        tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGGGGGGGGCGACGCCCGAAGCCCGTTCTTCCGGAAAAAGAGTGGCGGGTGCCTCCATCTGCGAACA‑GGCCCGCGCCCCCCGCCGTTCAAAACCTAATGCCAGtgtt                             >  2:82671/1‑148 (MQ=255)                                                                                                          tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGTAACGCCCGAAGGGGGTTTCGCTAGAAAAAGAGTGTTCGGGTCATTTGACAGCGGAAA‑GGACCGCGCGAGCCCGCGTTTTTAACACAATGGCAgcttatg                           >  2:766095/1‑145 (MQ=255)                                                                                                          tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAAAA‑AGCCAGCGCCAGCCGCCGTTTAAAAACTAATGGCAggttat                            >  1:198089/1‑145 (MQ=255)                                                                                                             ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTCAACGCCCGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCAACAGCGAAAA‑GGCCCGCGCCCGCCCCTGGTTAAAAACTAAAGGGAGGTCATCggg                        >  1:1066461/1‑150 (MQ=255)                                                                                                             ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAAGGGGGCAACTCGGGCCCACGGTTCCGTGGGAAACAGAAAGTTCGGTTAAAGCACCAGCCAACA‑CGCCAGGCCCCACCGCTTTTTTAAACCCAATCGGAGTTTAtaagt                        >  1:822920/1‑147 (MQ=255)                                                                                                               aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTTGTCGGGAAAAAGAATGTCCGGTTCATGCCTCAGCGAACA‑GGCCCGCGCCCCCCCCCGTTTTAAACCCCAAGGCCGTTCATCCGTgg                      >  1:572381/1‑151 (MQ=255)                                                                                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAGTCCGGGTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAAACCAATGGGCGGTCATCGGTGGCgg                   >  2:523805/1‑150 (MQ=255)                                                                                                                   tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGGCCGGGTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCCGTTTAAAACCCAAAGGCAggttctcggggggggg                  >  2:483212/1‑136 (MQ=255)                                                                                                                     cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAACCCGGGAATCGGGGTGGGGGAAAAAAGAAAGGCCGGTTCATGGATTAGCAAACACGCCCAGGGCCA‑CCGCCGTTTTAAAACTAAAGGCAGTTTAACGgggggggggg                >  1:586747/1‑151 (MQ=255)                                                                                                                     cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGAAGGCGGTTTGTCGAGAAAAAGAATGTCCGGGTCATGCATCAACGAACA‑GGCCCGCGCCCGCCGCTGGTTTAAACCTAAAGGGAGGTTATTGGGGGCGgggg                >  1:1141383/1‑151 (MQ=255)                                                                                                                       gTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGTAACGCCGGGAGGCGGGTTGTCGAGAAAAAAAAAGTCCGGGTCATGCAACAGCGAAAA‑GGCCCGCGCCCGCCGCTGTTTTAAACCCAAAGGCCGGTCAAcgggggcggggggg              >  1:146492/1‑149 (MQ=255)                                                                                                                            tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGAt         <  1:842453/151‑1 (MQ=255)                                                                                                                              tgCAGCCAAAATTCACGGCGGGTGAGGGTGTCAACGCCGGAAGTCGGTTCGTCGATAACCAGAATGTCCGGTTCATGCATCAGCGACCA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTTGCGTGGGAGGCGATAc       <  2:801238/151‑1 (MQ=255)                                                                                                                              tgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATAc       <  1:199515/151‑1 (MQ=255)                                                                                                                              tgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATAc       <  2:318577/151‑1 (MQ=255)                                                                                                                                   ccAAAATTCACGGCGGGTGAGGGGGTCACCGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTc  <  2:614911/151‑1 (MQ=255)                                                                                                                                   ccAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTc  <  2:284696/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                       CCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTC  >  NC_000913/827352‑827632

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