Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,271	0	T	G	47.1%	76.7 / 14.1	52	G165G (GGT→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base T (6/21);  new base G (24/0);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 5.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATC  >  NC_000913/1157126‑1157419                                                                                                                                                  |                                                                                                                                                       tCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGaaa                                                                                                                                                     >  1:123344/1‑149 (MQ=255) tCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAtt                                                                                                                                                   >  2:654631/1‑151 (MQ=255)  ccAGATTGGTCGTGAACTGACCAAGCCGGTTATCGTCCATACCCGTGATGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAtt                                                                                                                                                   <  2:998493/150‑1 (MQ=255)  ccAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAtt                                                                                                                                                   <  2:233281/150‑1 (MQ=255)            cGTGATCTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGAtctc                                                                                                                                         <  2:338311/150‑1 (MQ=255)            cGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCg                                                                                                                                        <  1:420935/151‑1 (MQ=255)             gTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCgg                                                                                                                                       <  1:783713/151‑1 (MQ=255)                    gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATATTTGGTCCCGGGtttt                                                                                                                                  >  2:330119/1‑149 (MQ=255)                       cAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATATTTGGTCTTGGGtttt                                                                                                                                  >  1:129342/1‑146 (MQ=255)                       cAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGATCTCTGGtttt                                                                                                                                  >  1:202304/1‑146 (MQ=255)                          gCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGCAAATAATCGATCTCCGGTTTTACATCT‑cc                                                                                                                         >  1:21759/1‑151 (MQ=255)                           ccGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCt                                                                                                                        <  1:1067841/151‑1 (MQ=255)                           ccGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCt                                                                                                                        <  2:556819/151‑1 (MQ=255)                                 aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCg                                                                                                                  <  2:379267/151‑1 (MQ=255)                                 aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCg                                                                                                                  <  2:197022/151‑1 (MQ=255)                                 aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTTCATCT‑CCtttt                                                                                                                     >  2:312154/1‑148 (MQ=255)                                   cGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCgg                                                                                                                 <  1:498114/150‑1 (MQ=255)                                     tCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTAATGGGTTTCGGGTTTTATTTCT‑CCtttttt                                                                                                                   >  1:606714/1‑144 (MQ=255)                                      ccATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGGTCTCGGGGTTTATATTT‑CCTTTTTCggggtt                                                                                                             >  1:25883/1‑147 (MQ=255)                                          aCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATtgtg                                                                                                          <  2:575719/150‑1 (MQ=255)                                                       gTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTTCATCT‑CCTTTTTCCGCATTGTGACCTTCC‑GTaa                                                                                              >  2:711425/1‑148 (MQ=255)                                                         gCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGATCTCGGGTTTTTCATTT‑CCTTTTTCCGCATTTTGACCTTTC‑CTAATTCgg                                                                                         >  2:842033/1‑151 (MQ=255)                                                           cGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTTAATCT‑CCTTTTTCGGCATTGTGACCTTTC‑GTAATGCgg                                                                                         >  2:755102/1‑149 (MQ=255)                                                                cGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATACTTGGTTTCCGGTTTTACATCT‑CCTTTTTCCGCCTTTTGACCTTCC‑CTAATGCGGAGCa                                                                                     >  2:303148/1‑148 (MQ=255)                                                                 gCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTg                                                                                 <  2:878337/151‑1 (MQ=255)                                                                        aTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGTCTTGGATTTTTCATCT‑CCTTTTTCGGCATTTTGGCCTTCC‑GTAATGCGGAGCAACTgcgc                                                                              >  2:1034515/1‑147 (MQ=255)                                                                          tCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGAtgc                                                                         <  2:965585/150‑1 (MQ=255)                                                                           cTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACATGGTTTTGGGTTTTTCAATT‑CCTTTTTCCGCCTTGTGACCTTCC‑GGAAAATGGGGCAACTGCGCGAtactg                                                                       >  2:164985/1‑151 (MQ=255)                                                                            tGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGAAAAATATTGGGTATCCGATTTTATATCC‑CCTTTTCCGGCATTGTGGCCTTCC‑GCAATGCGGGGCAAccaagcgaag                                                                          >  2:535693/1‑138 (MQ=255)                                                                               gcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGATTTCGGGTTTTTAATCT‑CCTTTTTCCGCCTTGTGGACTTCC‑GTAATGCGGAGGCACTGCGGGCTGCTgcgc                                                                    >  1:666976/1‑150 (MQ=255)                                                                                cgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGATCTCGGATTTTTAATAC‑CCTTTTTCGGCAATGTGAGCCTTC‑GTAAAGCGGGGCCAACCCGCGATGCCTcgcg                                                                   >  1:281815/1‑150 (MQ=255)                                                                                 gaagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGtt                                                                 <  2:202496/151‑1 (MQ=255)                                                                                  aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTa                                                                <  2:702654/151‑1 (MQ=255)                                                                                  aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGTCTCGGATTTTTCATCT‑CCTTTTTCCGCATTGTGACCTTCC‑GTAATGCGGAGCAAACGCGCGATGCTGcgcggt                                                                 >  2:452570/1‑148 (MQ=255)                                                                                       aaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGAGCTCGGGTTTTATATCT‑CCTTTTTCGGCCTTGTGGCCTTCC‑GGTATTGGGAGCAACCGCGCCAGTCTGCGGGTTTTTTcc                                                           >  1:1106892/1‑151 (MQ=255)                                                                                        aaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTTCTTGAGCTTGGGTTTTACATCTCCCTTTTCCGGC‑TTGTGGACTTCC‑CTACTGGGGGGGAACTGCGCGGTGCTGcgcggtgttttcc                                                          >  1:723761/1‑142 (MQ=255)                                                                                         aGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATACTGGGACTCCGGTTTTATATCT‑CCCTTTTCGGGATTGTGACCTTCC‑CTAATACGGGGCCACTGCGCGCTGCTGCGCGGTATGTcccc                                                         >  2:773381/1‑151 (MQ=255)                                                                                               ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGAt                                                   <  2:317948/151‑1 (MQ=255)                                                                                                  ttGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAATTAATGGGTCTCGGATTTTAACTCT‑CCTTTTTCCGCATTGTTAGCTTCC‑GTAATGCGGGGCAAACACGCGGTGGTGCGCGTTATGTTCCCCCGGATCgg                                                >  1:1032618/1‑151 (MQ=255)                                                                                                       gTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCTCGGGTTTTACATCC‑CCTTTTTCCGCCTTGTGACCTTCC‑CTAATGCGGGGCAAATGCGCGATGCTGCGCGGTATTTTCCCCCGGATCGGTTAc                                            >  2:1028109/1‑150 (MQ=255)                                                                                                        tGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGTAAATATTGGAACTCCGAGTTTAAATCT‑CCTTTTTCGGCATTGTGAACTTCCGGAAATCGGGGGCAAAT‑CGCGGTGCTGCGCGGTTTGTTCCCCTCGGTCg                                                 >  1:995048/1‑144 (MQ=255)                                                                                                           cGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACtggtg                                       <  2:1111071/151‑1 (MQ=255)                                                                                                            gTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTGGGGCTCCGGGTTTACATCT‑CCTTTTTCCGCATTGTGGCCTTCC‑GTAAATGGGGGGAACAGCGCGGTGCTGCGCGCTATGTTCCCCCCGGTCGGGTACtggggg                                      >  2:72862/1‑151 (MQ=255)                                                                                                             tactacACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATACTTGGTCTCGGGTTTTAAATCT‑CCTTTTTCGGCCTTGTGACCTTCC‑CTAATGCGGAGCAACTGCGCGGTGCTGCGCGGTATTTCCCCCCGGATCGGTTACtggtggg                                     >  2:8304/1‑150 (MQ=255)                                                                                                             tactacACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTTGATTTCGGTTTTTTAATCT‑CCTTTTTCGGCCTTGTGGACTTCC‑GCAATGCGGGGCAACTTCGCGGTGCTGCGGGTTATGTCCCCCCGGAGCGGGTATtggtggg                                     >  1:998803/1‑150 (MQ=255)                                                                                                               ctacACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTATTGGGTCTCGGGTTTTTCAACT‑CCCTTTTCCGCCTTGTGACCCTCC‑GGAATTCGGGGCAAACGCGCGATGGTGCGCGGTTTTTCCCCCCGGATCGGtttattggggaa                                    >  2:838035/1‑140 (MQ=255)                                                                                                                        tttACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTATCGGGTTTTAAATCT‑CCTTTTTCGGCCTTTTGGCCTTCC‑GTAATGCGGGGGAACCGCGCGATGCTGCGCGTTTTGTCCCCCCGGGTCGGTTACTTGGGGAAACTGACTCac                          >  2:49009/1‑151 (MQ=255)                                                                                                                               agGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCaccttacct                   <  2:920989/151‑1 (MQ=255)                                                                                                                                      agaAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCg            <  2:684035/151‑1 (MQ=255)                                                                                                                                      agaAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCg            <  2:886035/151‑1 (MQ=255)                                                                                                                                               gggTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAt   >  1:805090/1‑151 (MQ=255)                                                                                                                                                ggTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATc  <  1:7435/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                       TCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCTCGGATTTTACATCT‑CCTTTTCCGGCATTGTGACCTTCC‑GTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATC  >  NC_000913/1157126‑1157419

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