Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,271	0	T	G	44.7%	56.7 / 10.4	39	G165G (GGT→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base T (2/19);  new base G (17/0);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGC  >  NC_000913/1157122‑1157416                                                                                                                                                       |                                                                                                                                                   catATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGt                                                                                                                                                    <  1:18144/150‑1 (MQ=255)  atatCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATCCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTaa                                                                                                                                                  <  1:28818/151‑1 (MQ=255)    atCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGGGGCGGACTACACTGTTTTACAGGGGACAGAGAAACGGC‑GGGaaa                                                                                                                                                  >  1:73109/1‑149 (MQ=255)      ccAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAAtt                                                                                                                                               <  1:9154/150‑1 (MQ=255)      ccAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGaaaat                                                                                                                                               >  1:176689/1‑148 (MQ=255)        aGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACt                                                                                                                                            <  1:296471/151‑1 (MQ=255)                   gaactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGa                                                                                                                                  <  1:351440/150‑1 (MQ=255)                            aaGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATAATGGGTTTCCGGtttt                                                                                                                              >  1:405133/1‑145 (MQ=255)                             aGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACAtctc                                                                                                                       <  1:255054/151‑1 (MQ=255)                                   ttATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAATTACTGGGTCTCGGGTTTTTCATCTCCttttt                                                                                                                 >  2:405571/1‑150 (MQ=255)                                    tATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTc                                                                                                                 <  2:60142/150‑1 (MQ=255)                                        gTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATAATGGATATCGGGTTTTATATCTCCtctttttg                                                                                                              >  1:406826/1‑145 (MQ=255)                                                       gCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGATCTTGGGTTTTTCATCTCCATTTTCGGCCTTGTGACCTTTCGt                                                                                              >  2:312033/1‑150 (MQ=255)                                                              ccGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGa                                                                                      <  1:147241/151‑1 (MQ=255)                                                              ccGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGATTTTGGATTTTAAATATTCTTTTTCCGCCTTTTGGCCTTTCGTAATGCgg                                                                                       >  2:258204/1‑150 (MQ=255)                                                                    cGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACt                                                                                <  1:576031/151‑1 (MQ=255)                                                                     gCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTTATTGAACTCTGAATTTTCAAATCCCTTTCCGTGCTTGGGACCTTCCGTCATGCAGAGGAAca                                                                                >  2:87677/1‑149 (MQ=255)                                                                        ggCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTgcg                                                                             <  2:493509/150‑1 (MQ=255)                                                                                 gcgcGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTgcg                                                                   <  2:612376/151‑1 (MQ=255)                                                                                 gcgcGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTgcg                                                                   <  2:134579/151‑1 (MQ=255)                                                                                  cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAAATATTGGGTCTCGGGTTTTACATCTCCTTTTTCGGCATTGTGAACTTCCGTAATTCGGGGCAACCGCGCCAtgct                                                                      >  2:272167/1‑147 (MQ=255)                                                                                    cgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATTGATCTCGGATTTTTCAACTCCTTTTTCGGCATTGTGTCCCTCCCTAATTCGGAGCAACTGCGCGATGCTgcg                                                                   >  2:124724/1‑148 (MQ=255)                                                                                       agaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGTAATTACTGGATCTCGGGTTTTATATCTCCTTTTTTCGCATTTTGACCTTTCGTAATACGGAGCAACTGCGGGGTGCTGCGCGGTAt                                                             >  1:106584/1‑151 (MQ=255)                                                                                        gaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAATTAATTGCTGTTGGCTTTTAAAAATCCTTTTTCGGCCATGTGGCCCTTCGTAAATCGGGGCAAACTCGCCACTCTGCGCGt                                                                >  2:209745/1‑147 (MQ=255)                                                                                                gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCt                                                    <  1:178633/150‑1 (MQ=255)                                                                                                gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCATg                                                   <  1:217088/151‑1 (MQ=255)                                                                                                    gATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCTGGATc                                               <  2:243534/151‑1 (MQ=255)                                                                                                     aTTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCTGGATCg                                              <  1:127488/151‑1 (MQ=255)                                                                                                          ggTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGATCTCGGGTTTTATATCTCCTTTTCCGGCATTGTGGCCTTCCGTAATGCGGAGCAACTGCGCGGTGCTGCGCGTTATG‑TCCCCCCGGGTCgggt                                           >  1:581939/1‑147 (MQ=255)                                                                                                            tGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTAGGACTTGGATTTTAAATCTTCTTTTTCGGGATTGTGAACCTTCGTAATGCGGGGGAAATTCGCGATGCTGCGCGTTTTG‑TTCCCCCGGATCGGGTAc                                         >  2:365284/1‑149 (MQ=255)                                                                                                               cGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCTGGATCGGTTACtggtg                                    <  1:327526/151‑1 (MQ=255)                                                                                                               cGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAAATACTGGGTCTCGGGTTTTAAATCTCCTTTTCCCGCCTTTTGACCTTCCGTAATTCGGGGCAAACGCGCGATGCTGCGCGGTATG‑TCCCCCCGGGTCGGGTACtgggg                                    >  1:358478/1‑149 (MQ=255)                                                                                                                 tactacACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGGTCTCGGGTTTTACATTTCCTTTTTCGGCATTGTGACCTTCCGTAATGCGGGGCAAATGCGCGGTGCTGCGCGTTATG‑TCCCCCCGGGTCGGGTACtggggg                                   >  1:464032/1‑150 (MQ=255)                                                                                                                     acacTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCTGGATCGGTTACTGGTGGAAACt                              <  1:405707/151‑1 (MQ=255)                                                                                                                        cTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTGTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGAGCGTTATG‑TCCCCCTGGATCGGTTACTGTTGGAAACTGa                            <  2:143330/150‑1 (MQ=255)                                                                                                                         tGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATTATTGGGTTTTGGTTTTTAAACCCCCTTTT‑CGGCAATGTGAACTTCCGTAAATCGGGGCAACTGCGCGGTTCTTCGCGTTTTGTTCCCCCTGG‑TCGGGTACTGGTGGgaa                                >  2:52122/1‑145 (MQ=255)                                                                                                                            tttACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGATTTTCGAATTTAACTTTTCTTTTTCGGCATTGTGACCTTTCCTAATGCGGGGGAACTACGCGAAGCTTCGCGTTTTT‑TTCCCCCCGAGCGGGTACTGGGGGAAAATGACTccc                       >  2:440770/1‑149 (MQ=255)                                                                                                                               aCAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCTGGATCGGTTACTGGTGGAAACTGACTCacctt                    <  1:84980/151‑1 (MQ=255)                                                                                                                               aCAGAGGACAGAGAAACGGC‑GGGGAAAATACTGGATCTCCGGTTTTAAATCTCCTTTTTCCGCCTTTTGGCCTTCCGTAATGGGGGGGAACTTCGCGGTGCTGCGCGGTTTT‑TTCCCCCGGGTCGGTTACTTGGGGGAAAAGACTCccctt                    >  1:156641/1‑151 (MQ=255)                                                                                                                                  gagGACAGAGAAACGGC‑GGGGAAATTACTTGGTCTCGGGTTTTACATCTCCCTTTTCGGCATTGTGAACCTCCGTAATGCGGAGCAACTGCCCGGTGCTGCGCGTTTTG‑TCCCCCCGGGTCGGTTACTGGTGGGAACTGACTCacctt                    >  2:532906/1‑148 (MQ=255)                                                                                                                                                 gc‑ggGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGc  >  2:567479/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATG‑TCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGC  >  NC_000913/1157122‑1157416

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