Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	44.7%	52.4 / 13.6	38	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (5/16);  new base A (17/0);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 1.04e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508514‑3508811                                                                                                                                                       |                                                                                                                                                         aCGTTGTTTTTAAAGCCGGAGAGATCAATATGGAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCTAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAt                                                                                                                                                          <  1:608442/149‑1 (MQ=255)        ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCacacac                                                                                                                                                  >  2:179972/1‑150 (MQ=255)        ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCacacac                                                                                                                                                  >  1:534703/1‑150 (MQ=255)                ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGGAAAATCAAAAAC‑GGGAGcccc                                                                                                                                        >  1:469857/1‑151 (MQ=255)                   gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCcac                                                                                                                                     <  1:616275/151‑1 (MQ=255)                       gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAC‑CGTAGCCCCcac                                                                                                                                     >  2:472903/1‑147 (MQ=255)                          cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGa                                                                                                                              <  1:721284/150‑1 (MQ=255)                           aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGaaa                                                                                                                            <  2:253479/151‑1 (MQ=255)                                 cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACAC‑CGCAGGCCCCACAA‑CCGAAAGCgag                                                                                                                       >  2:333295/1‑148 (MQ=255)                                       cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTgcgc                                                                                                                 <  1:379600/150‑1 (MQ=255)                                         cATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGcaca                                                                                                              <  1:544899/151‑1 (MQ=255)                                             ggCCAGGGGTTAACCCGTTGGGATAAAACTCCATAGGGTGT‑TTCTGCAAAACGGAGGGAAGGAAGAAAAAAAGCGCCAATTTGGGGCCGGGGGGGCCCGAAAAAATTC‑CCCCC‑GGGTGGCCCCACAT‑CAGAAAGCGTGGTGCGCACACCg                                                                                                           <  2:588643/150‑1 (MQ=255)                                                gAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGc                                                                                                        <  1:530233/150‑1 (MQ=255)                                                        aaCCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACAC‑GGGAGCCCCCCCAT‑CAGAAAACGGGGGGGGCACACCGAGCTGCGTAg                                                                                                 >  2:661759/1‑149 (MQ=255)                                                         aCCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACAC‑GGGAGCCCCCACAT‑CAAAAAGCGGGGTGGGCACACCCCGCCGCGTAGCg                                                                                               >  2:165010/1‑150 (MQ=255)                                                            cGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat‑ca                                                                                          <  2:646109/151‑1 (MQ=255)                                                               cGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACAC‑GGTAGCCCCCCCAT‑CAGAAAGCGGGGGGCGCCCACCCAGCTGCGGAGCGAT‑CATcc                                                                                       >  1:17015/1‑151 (MQ=255)                                                                     aaaaCT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTCGGGCGGGGGAAAAATCACACAA‑CGAGGCCCCCCCATCCTCAAAGC‑TGGTGCGCACACCCACGAGCATAGGGAA‑CATccc                                                                                      >  2:261442/1‑146 (MQ=255)                                                                          t‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAC‑GGGAGCCCCCACAT‑CAGAAAGCGGGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTTCTGGTa                                                                             >  1:634105/1‑150 (MQ=255)                                                                            cATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATa                                                                             <  1:477164/149‑1 (MQ=255)                                                                               aGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAAcc                                                                         <  1:662738/150‑1 (MQ=255)                                                                                        ccTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTg                                                                  <  2:444744/148‑1 (MQ=255)                                                                                          tGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCAAACAA‑GGGAGCCCCCACAC‑CAGAAAGCGGGGGGGGCACACCAAGCGGCGTAGCGAG‑CATCCCCTTCTGGTAAAACAGCACCGGcccc                                                             >  1:530450/1‑151 (MQ=255)                                                                                            cAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGGCGGGAAAAATTAAAAAA‑GGGAGCACCCACAA‑AAGAAAAGGGGGGGGGCACAACGGGGAGCGGAGGGGA‑CAACCCTTTCTGGTAAAACAGCACTGGcacc                                                             >  1:596695/1‑149 (MQ=255)                                                                                              aaCCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACtt                                                         <  1:48569/151‑1 (MQ=255)                                                                                                 cGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACAC‑CGGAGCCCCCCCAT‑CAGAAAGCGGGGGGGGCACACCCAGCTGCGTAGCGAT‑CATCCCTTTCTGATAAAACAACACTGGCACCACTTcgc                                                      >  2:283003/1‑151 (MQ=255)                                                                                                 cGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACAC‑CGTAGCCCCCCCCT‑CAGAAAACGGGGGGCGCCCACCCAGCTGCGTAGCGAT‑CATCCCCTTCTGATAAACCCGCACCGGCACCCCCTcg                                                       >  2:172261/1‑150 (MQ=255)                                                                                                  ggaggg‑agGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACAC‑GGTAGCCCCCACAT‑CAGAAAGAGGGGGGCGCACCCCCAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGGCCCACTTcgc                                                      >  1:166084/1‑150 (MQ=255)                                                                                                     ggg‑aggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                                  <  1:228472/151‑1 (MQ=255)                                                                                                         aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACAC‑GGTAGCCCCCCCAT‑CAGAAAACGGGGTGCGCACACCCAGCTGCGGAGCGAT‑CATCCCCTTCTGATAAAACAGCACTGGCACCACCTCGCGGGGAtg                                               >  2:525425/1‑151 (MQ=255)                                                                                                          ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                              <  2:679176/151‑1 (MQ=255)                                                                                                           gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                             <  2:514689/151‑1 (MQ=255)                                                                                                             aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACAC‑GGGAGCCCCCCCAT‑CAGAAAACGGGGGGCGCACACCCCGCTGCGTAGCGAT‑CATCCCTTTCTGATAAAACAGCACTGGCACCCCCTCGCGGGGAtggggg                                           >  2:249745/1‑151 (MQ=255)                                                                                                                aaaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACAC‑CGTAGCCCCCCCAT‑CAGAAAGCGTGGGGGGCACACCCAGCCGCGTAGCGAT‑CATCCCCTTCTGATAAAACAACACCGGCCCCACTTCGCGGGGAtggtgggt                                         >  2:388104/1‑148 (MQ=255)                                                                                                                           cATTTTGTGGCGGTGGCGCGCCGGAAAATTTACACACGGGGAACCCCCCCA‑‑CACGAAACGGGGTGCGCCCACCCAGCTGCGGAGCGGTCCTTCCCTTCCGGTAAAACCGCCAC‑GGCACCACCTTGCCGGTATGGGGGCggtgcgcaaaagt                             >  1:8692/1‑142 (MQ=255)                                                                                                                                  tggcggtggcgCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                      >  2:44271/1‑151 (MQ=255)                                                                                                                                  tggcggtggcgCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                      >  1:287048/1‑151 (MQ=255)                                                                                                                                           cgcgcgGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGa             <  1:60059/151‑1 (MQ=255)                                                                                                                                                      aTTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  >  2:13869/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                         ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACT‑CATAGGGTGTCTCCTGCAAACCGGAGGG‑AGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAC‑GGTAGCCCCCACAT‑CAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGAT‑CATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508514‑3508811

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