Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,747	0	T	G	42.9%	58.7 / 15.6	44	T255P (ACC→CCC)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (5/19);  new base G (18/0);  total (25/19)
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

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAC  >  NC_000913/646599‑646895                                                                                                                                                     |                                                                                                                                                           aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑c                                                                                                                                                         <  1:31700/151‑1 (MQ=255)          aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑cccc                                                                                                                                            <  1:606640/151‑1 (MQ=255)              ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCtt                                                                                                                                         <  2:466947/150‑1 (MQ=255)                aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑T‑GA‑AT‑G‑GCCCCctttt                                                                                                                                      >  2:619753/1‑151 (MQ=255)                aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GA‑AT‑G‑CCCCCctttt                                                                                                                                      >  1:665038/1‑151 (MQ=255)                 tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GA‑AT‑G‑ACCCCctttt                                                                                                                                      >  1:756660/1‑150 (MQ=255)                  cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTg                                                                                                                                     <  1:305724/150‑1 (MQ=255)                   aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GA‑AT‑G‑CCCCCttttt                                                                                                                                      >  1:516966/1‑148 (MQ=255)                      tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GG‑AA‑G‑CCCCCcttttttaa                                                                                                                                  >  1:259793/1‑145 (MQ=255)                       gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑G‑GG‑AA‑G‑CCCCCCTTTTTTaata                                                                                                                                >  2:455203/1‑150 (MQ=255)                         ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTGC‑GAGGT‑C‑CCCCCtttttaa                                                                                                                                    >  1:710122/1‑144 (MQ=255)                           aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGG‑C‑GA‑AT‑G‑CCCCCCTTTTTCaataa                                                                                                                               >  2:283757/1‑147 (MQ=255)                            aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGttt                                                                                                                         <  2:336980/151‑1 (MQ=255)                            aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTG‑C‑GG‑GT‑G‑CCCCCTTTTTGCGAAAA‑TGtt                                                                                                                          >  1:281886/1‑150 (MQ=255)                             ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTg                                                                                                                        <  2:556799/151‑1 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTgc                                                                                                                       <  1:238093/150‑1 (MQ=255)                                 ttAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑T‑GA‑AT‑G‑CCCCCCTTTTTCAAAAA‑AGTTTgcggc                                                                                                                    >  1:144138/1‑149 (MQ=255)                                   aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GG‑AT‑G‑CCCCCCTTTTGCCATAA‑TGTTTTCGcc                                                                                                                    >  1:775688/1‑149 (MQ=255)                                     aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTTCGCCTccc                                                                                                                >  1:344797/1‑149 (MQ=255)                                       cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGACGAACTGCCGGAGATAATCGAGATCGGCGGGGGAT‑T‑GA‑ATCC‑CCCCCTTTTTGCAA‑AA‑TGCGTTCGCCTctc                                                                                                                >  2:836576/1‑147 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAg                                                                                                         <  2:98560/149‑1 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GA‑AT‑G‑CCCCCCTTTTTCCATAA‑AGTTTTCGCCCCGCGCCGGAGcc                                                                                                       >  2:489133/1‑151 (MQ=255)                                                 aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GC‑AT‑G‑CCCCCCTTTTTTCATAA‑TTTTTTCGCCTccccccggagcc                                                                                                       >  2:476771/1‑137 (MQ=255)                                                   ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAgcgc                                                                                                  <  1:73275/151‑1 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                       <  2:257934/150‑1 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                       <  1:294930/150‑1 (MQ=255)                                                                       cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GG‑AT‑G‑CCCCCCTTTTTCAATAA‑TTTTTTCGCCCCGCGCCGGTGCCCGCGCAGGCCCCCCCCCCCAcc                                                                                 >  1:949119/1‑148 (MQ=255)                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑T‑GG‑AT‑G‑CCCCCCTTTTTCAATAA‑TGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCcccgc                                                                             >  1:230382/1‑151 (MQ=255)                                                                              gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGGGATAATCGAGATCGGCGGGGGTT‑T‑GA‑AA‑G‑CGCCCCCTTTGCACAAA‑TCTTTTTCGCGCCCGCCGGGGCCCGGGCAGGCCCCCCCCCCCACcccgcgc                                                                           >  2:928578/1‑146 (MQ=255)                                                                              gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGCG‑T‑TG‑AT‑G‑TCCCCCTTTTTTTATTA‑TAGTTGCGCCCCCCGCCGTTGCCACCGCAGGGCCCCCTCGCCCCCCCGCGAt                                                                          >  1:396/1‑148 (MQ=255)                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GA‑AT‑GACCCCCTTTTTTC‑ATAA‑TGTTTGCGCCCCCCGCTGTCGCCAGCGCCGGGCCCCCCCCCCACCCCGCGGAGCCACGTTgggtt                                                             >  1:723100/1‑146 (MQ=255)                                                                                                 cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                                    <  2:221360/151‑1 (MQ=255)                                                                                                 cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GG‑AT‑G‑CCCCCCTTTTTTAATAATTTTTTTCGCC‑CCCCCCGTCGCCCGCGCCGGCCCCCCCCCCCCCCCCCCCCAGCAACGGTGGTATCGcc                                                         >  2:844407/1‑146 (MQ=255)                                                                                                    tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                  <  1:691343/150‑1 (MQ=255)                                                                                                         acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTTCGCCTCCCGCTGTAGCCAGCGCAGGCCCCCCCCGCCACCCCGCGGTGCAACGTTGGGATCGCCGTTGATCGcc                                               >  1:314276/1‑148 (MQ=255)                                                                                                          cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑CGGA‑AT‑C‑CCCCCTTTTT‑CAATAA‑TGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGGCCCCCCCGCCCCCCCCCGAGGGCACGTTGGTATCGCCGTTTAttcccccc                                            >  1:162910/1‑143 (MQ=255)                                                                                                                    cGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCaa                                 <  1:367483/151‑1 (MQ=255)                                                                                                                     gAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑C‑GG‑AT‑G‑CCCCCCTTTTTCAATAA‑TGTTTTTGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGAGCAACGTTGGTATCGCCGGTGATCGCCCTCAGTAgggg                                    >  1:79834/1‑145 (MQ=255)                                                                                                                           ccGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                          <  1:870224/151‑1 (MQ=255)                                                                                                                                 aTAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAt                    <  1:685489/151‑1 (MQ=255)                                                                                                                                       gagaTCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGt              <  2:432123/151‑1 (MQ=255)                                                                                                                                          aTCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa           <  1:594003/151‑1 (MQ=255)                                                                                                                                                gggggTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt     <  1:967105/151‑1 (MQ=255)                                                                                                                                                   ggTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAc  >  1:567505/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                           AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTT‑C‑GA‑AT‑G‑CCCCCTTTTTGCAATAA‑TGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAC  >  NC_000913/646599‑646895

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