Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,111,259	Δ1 bp	coding (718/900 nt)	rfbD ←	dTDP‑4‑dehydrorhamnose reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,111,259	0	G	.	89.5%	58.5 / 2.4	19	coding (718/900 nt)	rfbD	dTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (12/5);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACC‑TGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTC  >  NC_000913/2111052‑2111501                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   caTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCGCGACAACC‑TGACAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGCGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAAcc                                                                                                                                                                                                                     >  1:188832/1‑238 (MQ=255)                                                            gctcttcatatcccAACCATGCCAGTCAGGCAAGACAAGCGCAAAGGTCTTCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGCCGACGAGATGGTTTAGGATAGGCTGTTGTTGGTACTGCGTGGAGCTGGTTAATTGCAAGGA‑AATGCCGGTTTTGCGCGCCTCGTCAAAAACCAGCGCATCATAATCGTACCACGTGGTGGTACGACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGtt                                                                                                                                               <  1:323531/238‑1 (MQ=255)                                                                         naaCC‑TGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCGGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGACTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGaca                                                                                                                                 >  1:275404/2‑250 (MQ=255)                                                                                   ngTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGAGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAAcc                                                                                                                                                                                                                     >  1:116566/2‑157 (MQ=255)                                                                                         gCAAGAAAAGCGCAAAGTTCTGCTGAAATTTTTCTGAATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCACATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTg                                                                                                                 >  1:379390/1‑251 (MQ=255)                                                                                                     cAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTc                                                                                                                                                                                                                            >  1:96675/1‑132 (MQ=255)                                                                                                     cAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGNGGACGACGANNTGNTGTAGNANNNGNNNTTNTTGGTACTGCGTTGAGCTTNTTNNGTGCAANGG‑AATGCCTGCTTTGCGCGCCTCTTc                                                                                                                                                                                                                            <  2:96675/132‑1 (MQ=255)                                                                                                                          cTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTgcgcgcn                                                                                                                                                                                                                                 <  1:79501/106‑2 (MQ=255)                                                                                                                                             ntgtgGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGAACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCTCAAACCAGCGCAGCATAATCGTAACAGGTTGTGGTACCACTGGCTACCAAATGGTAc                                                                                                                                                                  >  1:50849/2‑150 (MQ=255)                                                                                                                                                    cgacgaGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCA‑CGGGAATGACTGCTTTGCGAGCCTCTTCCCACACCAGCGCAGCCTAATCGTACCAGGTTGTGGTAcc                                                                                                                                                                                     >  1:216282/1‑124 (MQ=255)                                                                                                                                                             tgtgtANGNTAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCNTGCTTTGCGCNCCTNNTCAANANCNAGCGCAGCATAATCGTACCAGGTTGTGGTNCCACTGGCTACCAAATGGTACAANNCTNNNACNNCCGGt                                                                                                                                                <  2:335783/151‑1 (MQ=255)                                                                                                                                                             ggTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTCCAAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGCCATCCGGt                                                                                                                                                >  1:335783/1‑152 (MQ=255)                                                                                                                                                             ggTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGc                                                                                                                                                                               >  1:146512/1‑121 (MQ=255)                                                                                                                                                                gtAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAACAAACAGCGCAGCATACTCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGGGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCACACTGATCGTTAATAACCGCTAATTCTTCaagc                                          >  1:310851/1‑248 (MQ=255)                                                                                                                                                                         gCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGcc                                 <  1:345375/251‑1 (MQ=255)                                                                                                                                                                                  ggTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGAAATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGCTCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTaa                         >  1:333018/1‑250 (MQ=255)                                                                                                                                                                                   gTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAAc                        >  1:189353/1‑250 (MQ=255)                                                                                                                                                                                               gCTTGGTGAGTGCA‑GGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTtgtg                                                                                                                                                                                          <  1:231690/76‑1 (MQ=255)                                                                                                                                                                                                        gTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGAACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGAATGTGCTGTACACTCAGACAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttc  >  1:446220/1‑251 (MQ=255)                                                                                                                                                                                                            naagggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGaca                                                                                                                                 >  1:265802/4‑121 (MQ=255)                                                                                                                                                                                                            caagggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGaca                                                                                                                                 <  2:265802/118‑1 (MQ=255)                                                                                                                                                                                                            caagggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACaa                                                                                                            >  1:420020/4‑142 (MQ=255)                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   CATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACC‑TGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTC  >  NC_000913/2111052‑2111501

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