Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,478,324	0	T	G	43.9%	66.7 / 10.0	41	G33G (GGT→GGG)	ynbB	putative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (3/20);  new base G (18/0);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 7.52e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTT  >  NC_000913/1478184‑1478474                                                                                                                                               |                                                                                                                                                        tAATTGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACa                                                                                                                                                <  1:794640/150‑1 (MQ=255)     cGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATC‑GGGGA‑GAAAAAGGGGGGGGGGTACAtttct                                                                                                                                           >  1:27148/1‑148 (MQ=255)     cGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGA‑GAAAAGGGGGGGGGGGTAAATTACg                                                                                                                                           >  2:712341/1‑151 (MQ=255)           cGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCg                                                                                                                                     <  1:993696/151‑1 (MQ=255)             gAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑a                                                                                                                                  <  1:205524/151‑1 (MQ=255)                tgtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCC‑GGAGA‑GAAAAGGGGGGGAGGGTAGATTTCGTATTTGttact                                                                                                                                >  1:275856/1‑147 (MQ=255)                     ttcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATC‑GGAGA‑GAAAAGGGGGGGAGGTTACATTTCGTATTTGT‑ACCtggggg                                                                                                                          >  1:288175/1‑151 (MQ=255)                          aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCAt                                                                                                                     <  2:200317/151‑1 (MQ=255)                               cgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAAC‑GGAGA‑GGAAAGGGGGGGGGGGTGCATTTAGGATTTGTAACCTGggggggatcg                                                                                                                    >  1:1028422/1‑144 (MQ=255)                                  tGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGA‑GAAAAGGGGGGGGGGGTACATTATGTTTTTGT‑ACCCGGTGGGGGGTCGTCCttt                                                                                                             >  1:295058/1‑151 (MQ=255)                                     aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGtt                                                                                                          <  1:435049/151‑1 (MQ=255)                                         taATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑CGAGG‑GAAAAGGGGGGGGGGGTACATTACGTATTCGT‑ACCTGGGGGGGGATGGGCATTTGttttt                                                                                                       >  2:999005/1‑150 (MQ=255)                                          aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTc                                                                                                      <  1:763003/150‑1 (MQ=255)                                          aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAAC‑GGGGG‑GAAAAGGGGGGGAGGGTACATTACGTATTTGT‑TCCTTGTGGGGGATCGTCATTTGttttt                                                                                                       >  1:893504/1‑149 (MQ=255)                                                ggAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTg                                                                                               <  1:1033595/151‑1 (MQ=255)                                                   aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGTGTGAGGTTACATTACGTATTCGT‑ACCTGGTCGGGCATCGTCATTTGTTTTTCAATGATGAtt                                                                                            <  1:73707/151‑1 (MQ=255)                                                      aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGA‑GAAAAGGGGGGGGGGGTCCATTACGTATTTGT‑TCCCGGGGGGGGCTCGACACTTTTTTTTTAATGGTGAttt                                                                                           >  2:493100/1‑149 (MQ=255)                                                          ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGG‑GAAAAAGGGGGGGGGGTAAATTACGTATTTGT‑ACCTCGGGGGGGGGCGTCCTTTGTTTTTTTATGGTGGttt                                                                                           >  2:126317/1‑145 (MQ=255)                                                          ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCa                                                                                         <  2:586623/147‑1 (MQ=255)                                                                 cacTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATC‑CGGGA‑GAAAAGGGGGGGAGGGTACATTACGTATTCGT‑ACCCGGTGGGGGATCGTCATTTTTTTTTTTATTGTGGTTTTTGGTCCAc                                                                                  >  1:761400/1‑147 (MQ=255)                                                                  acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTTCGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTGCACg                                                                                 <  1:928497/147‑1 (MQ=255)                                                                  acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAAC‑TGGGA‑GAAAAGGGGGGGGGGGTACATTACGTATTCGT‑TCCTGGGGGGGGCTCGTCATTTGTTTTTTTATTGGGGTTTTTGGTCCACCCTgg                                                                             >  1:720442/1‑151 (MQ=255)                                                                     tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGG‑GAAAAGGGGGGGGGGTTA‑ATTTCGTATTTGT‑ACCTGGTGGGGGGTCGTCATTTTTTTTTTTATGGTGGTTTTAGGGCCCCGCTgg                                                                             >  2:358914/1‑148 (MQ=255)                                                                        cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacg                                                                       <  2:70297/151‑1 (MQ=255)                                                                        cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgac                                                                        <  1:217333/150‑1 (MQ=255)                                                                         gcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgt                                                                      <  1:532317/151‑1 (MQ=255)                                                                          cgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTAT‑GCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttg                                                                    <  2:971541/151‑1 (MQ=255)                                                                             tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATC‑TGAGA‑GAAAAGGGGGGGAGGGTACATTTCGTATTTGT‑TCCTGGGGGGGCATCGTCCTTTTTTTTTTTATTGGGGTTTTTGGGCCACGCTGGGtgacg                                                                       >  1:344481/1‑146 (MQ=255)                                                                                   ttttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATC‑GGGGA‑GGAAAAAGGGGGGGGGTAAATTATATATTCGT‑ACCTGGGGGGGCATCGTCATTTGTTTTTTTATGGGGATTTTTGGGCCACGCTGGGTGACGTTGGCGttttt                                                            >  2:1012703/1‑151 (MQ=255)                                                                                         cAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCC‑CGAGAGGAAAAGGGGGGGGGGTT‑CATTTCCTATTCGT‑TCCTCGGGGGGGAGCGGCATTTGTTTTTCTTTGGGGATTTTAGGGCCCCGCCGGGTGTCGTTGTTGttttttttgct                                                      >  1:70941/1‑146 (MQ=255)                                                                                               tGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGA‑GAAAAGGGGGGGGGGGTACATTTCGTATTTGT‑ACCTGGGGGGGGGTCGTCATTTGTTTTTTTATTGTGATTTTAGGGCCACGCTGGGTGACGTTGACGGTTTTTTCGCTGATaa                                                 >  2:256074/1‑150 (MQ=255)                                                                                                gATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGt                                               <  2:813812/151‑1 (MQ=255)                                                                                                   aaaCCGTTTTCTGCTATGGCGGTTTCC‑GGGGA‑GGAAAAGGGGGGAGGGTACATTATGTATTTGT‑AACTGGTGGGGGGTCGGCCTTTGTTTTTTTATTGTGGTTTTTGGTCCCCCCCGGGTGGCGTTTACGGTTTTTTTTCTGGTAAGtttt                                            >  1:467976/1‑151 (MQ=255)                                                                                                        gTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGc                                       <  1:915195/151‑1 (MQ=255)                                                                                                            tCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCTCGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCAtt                                    <  2:380713/150‑1 (MQ=255)                                                                                                                    ggCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAata                           <  1:554417/151‑1 (MQ=255)                                                                                                                    ggCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAata                           <  2:432841/151‑1 (MQ=255)                                                                                                                      cGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATAtt                         >  2:793873/1‑151 (MQ=255)                                                                                                                          tACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTGAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTAc                     >  1:991193/1‑151 (MQ=255)                                                                                                                               ggAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCtt                 <  1:228061/151‑1 (MQ=255)                                                                                                                                      aaaaGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATc            >  1:861649/1‑150 (MQ=255)                                                                                                                                               tGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAtttt  >  1:330364/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                        TAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGA‑GAAAAGGGGGTGAGGTTACATTACGTATTCGT‑ACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTT  >  NC_000913/1478184‑1478474

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: ‑