Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,765	0	G	T	45.2%	45.9 / 27.5	43	Q249K (CAA→AAA)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base G (5/18);  new base T (19/0);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.06e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGA‑GCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCA  >  NC_000913/646617‑646908                                                                                                                                                        |                                                                                                                                                  aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCa                                                                                                                                                 <  2:611916/151‑1 (MQ=255)       ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑AATGACCCCTTTTTTCaata                                                                                                                                              >  1:619252/1‑148 (MQ=255)        tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑TATTGCCCCTTTTTTtaataa                                                                                                                                             >  1:136029/1‑148 (MQ=255)           ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTg                                                                                                                                       <  1:465423/151‑1 (MQ=255)           ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑AATGGCCCCCTTTTTTAATAAAGtt                                                                                                                                         >  1:342404/1‑149 (MQ=255)            taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑T‑AATGCCCCCTTTTTTCAAAAATAttt                                                                                                                                        >  1:345049/1‑149 (MQ=255)             aatTAAAGCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                                      <  2:238397/150‑1 (MQ=255)              atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAGCTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑GATGCCCCCCTTTTTCCATAATGttttc                                                                                                                                      >  2:441087/1‑147 (MQ=255)              atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTgcgc                                                                                                                                    <  1:468000/151‑1 (MQ=255)              atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑AAGCCCCCCTTTTTGCataatttt                                                                                                                                          >  2:28794/1‑139 (MQ=255)              atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑GATGCCCCCCTTTTTCCATAATttttt                                                                                                                                       >  1:295932/1‑147 (MQ=255)                 aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCt                                                                                                                                  <  1:404605/150‑1 (MQ=255)                 aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑AAAAACCCCCTTTTTTAATAAAGTTTTCGCCTc                                                                                                                                 >  1:444673/1‑151 (MQ=255)                      tgtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGT‑G‑GATGCCCCCCTTTTTTAATAATGTTTTCGcc                                                                                                                                   >  1:647938/1‑144 (MQ=255)                           tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑AATGCCCCCCTTTTTTAATAATGTTTGCGCCTCCCGCTGTAg                                                                                                                        >  1:351321/1‑150 (MQ=255)                               aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑AATGCCCCCCTTTTTTAATAAAGTTTTTGCCCCTCGCCGTAGCCCgc                                                                                                                   >  2:833427/1‑151 (MQ=255)                                  cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAATGCCCCCTTTTT‑CCATAATTTTTTCGCCCCGCCCTGTTGcc                                                                                                                      >  1:434561/1‑145 (MQ=255)                                    ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑G‑ATGGCCCCCTTTTTG‑AAAAATTTTTGTGCCTcgc                                                                                                                               >  2:172051/1‑134 (MQ=255)                                    ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑GAAGCCCCCCTTTTTTAATAATTTTTTGGCCCCGCGCCGTAGCCa                                                                                                                     >  2:335441/1‑144 (MQ=255)                                        aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGATC‑G‑GATGCCCCCCTTTTTTAATAATCGTTTCGCCTCCGGCTGTAgttcgcagaagcccc                                                                                                          >  1:945395/1‑137 (MQ=255)                                         ggATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGccccc                                                                                                         <  2:561581/151‑1 (MQ=255)                                             cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                                      <  1:652839/150‑1 (MQ=255)                                             cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑GATGCCCCCCTTTTTCCAAAATGTTTTCGCCTCCCGCCGTAGCCAGCGCAGGcccccccc                                                                                                      >  1:795187/1‑148 (MQ=255)                                                     cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑C‑AAAGCCCCCCTTTTTCCATAATGTTTTTGCCTCGCGCCGTAGCCAGCCCAGGCCCCCCCCGCccccc                                                                                               >  2:275617/1‑148 (MQ=255)                                                              cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAAc                                                                                   <  2:943907/151‑1 (MQ=255)                                                                  tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑AATGCCCCCTTTTTTCAATAATTTTTTCGCCCCTCGCTGTAGCCAGCGCAAGCCCCCCTCCCCACCCCGCGAGG‑CAACGtt                                                                                >  1:23626/1‑150 (MQ=255)                                                                   ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGtt                                                                                <  2:311733/149‑1 (MQ=255)                                                                    gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGt                                                                             <  1:19280/151‑1 (MQ=255)                                                                     aCTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTG‑G‑AATGCCCCCCTTTTTCCATAATGTTTGCGCCCCCCGCCGTAGCCAGCGCAGGCCCCCCTCCCCCCCCCGCGCTG‑CAACGtt                                                                                >  2:667693/1‑147 (MQ=255)                                                                           aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCg                                                                      <  1:705447/151‑1 (MQ=255)                                                                                 ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑CATGGCCCCCTTTTTCAATAATATTTTCGCCTCGCCCCGTAGCCCGCGCCGGCCCCCCCCCCCCCCCCCCCCTG‑GAACGTTGGGATTGCCGCTGAtt                                                                >  2:901531/1‑150 (MQ=255)                                                                                   cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCg                                                               <  2:917929/150‑1 (MQ=255)                                                                                     aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCa                                                            <  2:15247/151‑1 (MQ=255)                                                                                     aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCa                                                            <  1:48310/151‑1 (MQ=255)                                                                                            cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAg                                                     <  2:414105/151‑1 (MQ=255)                                                                                                  cGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑AATGGCCCCCTTTTTCCATAATATTTTTGCCCCGCGCCGGAGCCCGCGCAGGCCCCCCCCGCCCCCCCCCGATG‑CAAAGGTGGTATCGCCGGTGATAGCCCTCACTAgg                                                    >  2:787968/1‑146 (MQ=255)                                                                                                        gCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGA‑GCAAGGTAt                                         <  1:715873/150‑1 (MQ=255)                                                                                                            gagaTAATCGAGATCGGC‑GGGGGGTT‑G‑AAAGCCCCCCTTTTTCTAAAATATTTTCGCCCCCCGCTGTAGCCAGCGCCGGGCCCCCCCCCCCCCCCGCCGTGCCAACG‑TGGGATTGCCGGTGATTGCCCACAGTAGGAGGCAAggatc                                         >  2:849976/1‑144 (MQ=255)                                                                                                                             c‑ggGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGA‑GCAAGGTATCGAGCAATGCCAGTTCAGGATc                   >  1:596689/1‑151 (MQ=255)                                                                                                                                     tC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGA‑GCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcccc            >  2:589772/1‑151 (MQ=255)                                                                                                                                              cccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGA‑GCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATcc      >  2:350922/1‑150 (MQ=255)                                                                                                                                                cccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGA‑GCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagc   <  1:820311/151‑1 (MQ=255)                                                                                                                                                 ccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGA‑GCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagca  <  2:206814/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  AAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGA‑GCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCA  >  NC_000913/646617‑646908

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