Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,361,249	G→A	Q56* (CAG→TAG)	rhmD ←	L‑rhamnonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,361,249	0	G	A	96.0%	71.6 / ‑6.3	25	Q56* (CAG→TAG)	rhmD	L‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (14/10);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACACAAGAAATCGTATTCATCACCAGGCCACCAGAGCCGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCG  >  NC_000913/2361006‑2361489                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  cGACACAAGAAATCGTATTCATCACCAGGCCACCAGAGCCGGAGTAATTCAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTTACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCGAAGAACGTTAATGCCAAATGACTAGCGTGAc                                                                                                                                                                                                                                            <  1:221897/251‑1 (MQ=255)                    catcaCCAGGCCACCAGAGCCGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGGTCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAAAGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTAtt                                                                                                                                                                                                                         >  2:75962/1‑251 (MQ=255)                       ccccAGGCCACCAGAGCCGGAGTAATACAGGGGGGCACTGAGCATGTGTTCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATCTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACCTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGGGACTGCTCGTAATCGCGGTAttt                                                                                                                                                                                                                        <  2:84085/247‑1 (MQ=255)                                 ccaGAGCCGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑t                                                                                                                                                                                                           <  2:65357/251‑1 (MQ=255)                                                   aGGGGGGCACTGAGCATTTTAGCGTGGATCAGGTTGATATCACTGACACATTTACCCTCAAAGCAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGATCTGGACGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAATTGACTAGCGTGACTGCTCGTAATCGCTGTATTTACTCATCGGGG‑TGGCAAGATGATCGTCTAt                                                                                                                                                                                         <  2:79083/251‑1 (MQ=255)                                                    gggTGGCACTGAGCATTTGATCGTGGAGCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTATCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCTATATGACCGGCTATc                                                                                                                                                                                        <  2:136938/251‑1 (MQ=255)                                                                gCATTTGAACGTGTATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTACGATGTTTTTCGAAAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTTCTCGTAATCGCGGTATTTACTCATCGGGG‑TGTCAATATGATCGTCTATCCAGTGATTCGcc                                                                                                                                                                            <  2:246736/251‑1 (MQ=255)                                                                         cGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCg                                                                                                                                                                   >  1:25059/1‑251 (MQ=255)                                                                                    ttGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCAc                                                                                                                                                         <  1:141865/250‑1 (MQ=255)                                                                                     tGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGTG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCAcc                                                                                                                                                        <  1:15033/250‑1 (MQ=255)                                                                                                acacATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGcc                                                                                                                                            <  1:111135/251‑1 (MQ=255)                                                                                                 cacaTTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCt                                                                                                                                           >  1:1749/1‑251 (MQ=255)                                                                                                                       tttAAGATGTTTGTCGACAATAAAACAGCCCATTTCACCGGCTGACGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCAACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAAACGCGGTATTTACTCATCGGGG‑TTGCCATATGATCGGCTATCCAGTGGTTCGCCCCCGAGTCGTGCTACTAACCGCAGCCAACGCATTTTTCTGCTGTTGCACCGccgg                                                                                                                     >  2:146214/2‑249 (MQ=255)                                                                                                                       gTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAAAAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑AGGCAATATGATCGTCTATCAAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGcc                                                                                                                       >  1:229625/1‑249 (MQ=255)                                                                                                                            gATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAc                                                                                                                >  1:161435/1‑251 (MQ=255)                                                                                                                              tGTTTTTCGACACTAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAAAGTTAATGCCAAATGACTAGCGTGACTGCCCGTAATCGCGGTATTTACTCATCGGGG‑CGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAcc                                                                                                               >  2:54175/1‑250 (MQ=255)                                                                                                                                                                                tCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCCCCGCCAGTAAACCAGGCCCGAACCTGTTTAATTTTTGGTAGGGACATGATGTTCTCCATTGtt                                                             >  2:80672/1‑250 (MQ=255)                                                                                                                                                                                      cTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTACTTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGc                                                      >  2:19220/1‑251 (MQ=255)                                                                                                                                                                                           cGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGt                                                  >  2:13226/1‑250 (MQ=255)                                                                                                                                                                                                            tCAAAAACCAACGTGCCTCGAACGTTAATGCCAACTGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCCGTAACCCCGGCGCGCACCTGTTTAATTTTTGGTCGGGTCATGATGTTCTCCATTGTTATTAGGCTTTTAAGTCAAAGGGACTTTTc                                >  2:156669/1‑251 (MQ=255)                                                                                                                                                                                                                caccaacctggCCAAGAACTTTAATCCCAAATGACTGGCGTGACTGCTCGAAACCGCGGTATTTCCTCACCGGGG‑TGGCAATATAACCGTCTATCCAGTGATCCCCCCCCTGGTCGTGATAACCCCCGCGGCCGGCGCCTTTGTCTGTTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATc                            <  1:59752/241‑1 (MQ=255)                                                                                                                                                                                                                                    ttAATGCCAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATcc        >  2:73268/1‑251 (MQ=255)                                                                                                                                                                                                                                    tccaatttaaaggACTAGCGTGACTTCTCGAAACCGCGTTTTTTACTC‑TCGGGGGTGCAAATATGTCCGTCTATCCAGTGATTCCCCCCCGGGTCGATAAAACCCCCGCCGCCGCCTCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATc         <  2:207151/238‑1 (MQ=255)                                                                                                                                                                                                                                          ccAAATGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGgcg  >  1:149495/1‑251 (MQ=255)                                                                                                                                                                                                                                            aaaTGACTAGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATc         >  1:207151/1‑242 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  CGACACAAGAAATCGTATTCATCACCAGGCCACCAGAGCCGGAGTAATACAGGGTGGCACTGAGCATTTGATCGTGGATCAGTTTGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGG‑TGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCG  >  NC_000913/2361006‑2361489

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