Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,939,297	G→T	G137V (GGG→GTG)	fucU →	L‑fucose mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,939,297	0	G	T	96.5%	82.5 / NA	29	G137V (GGG→GTG)	fucU	L‑fucose mutarotase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (13/15);  minor base A (1/0);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTAAATATATCC  >  NC_000913/2939057‑2939542                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      ttATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCg                                                                                                                                                                                                                                             >  1:312978/1‑251 (MQ=255)   aTCCCGTTATTTGAACTGGACAGTTATGCCCCGCCGCTGGTGATGATGGCGGCGGTAGAAGGGGACACGCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCTCTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTa                                                                                                                                                                                                                                           <  2:477885/251‑1 (MQ=255)           aTTTGAACTGGACAGTTATCCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCAt                                                                                                                                                                                                                                    >  1:369851/1‑250 (MQ=255)                 aCTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGATAACGGAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATGATACGCATCAATTGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCTCTCGCTATCAAAGGCGAACTAGCGAAGGACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGCAATCTCATACAGGt                                                                                                                                                                                                                              >  2:436408/1‑250 (MQ=255)                    ggACAGTTATGCACCGCCGCGGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTTAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGTTCTTTATCGCAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTAc                                                                                                                                                                                                                            <  1:434185/249‑1 (MQ=255)                                aCCGCCGCTGGTGATGATGGCGGCGGTAGAAGGAGACACTCTCGATCCTGAAGCAGAACGACGTTACCGCAATGCGCTTTCAATACAAGCACCGTGTCCTGAACTCCTCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTACAAAAAGTGGTAACACCGACATCTCATACCGGTACGCCCGCATGAcgcg                                                                                                                                                                                                              >  2:459145/1‑251 (MQ=255)                                  cgccgcTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAACCGTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcggg                                                                                                                                                                                                            <  1:30296/251‑1 (MQ=255)                                          gTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATc                                                                                                                                                                                                    <  1:463343/251‑1 (MQ=255)                                              tgatgGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                                                                                                >  1:354012/1‑251 (MQ=255)                                                atgGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACTACGTTACCGTAATGCGCTTTCACGACGAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGTGCGTAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatga                                                                                                                                                                                              <  1:514013/251‑1 (MQ=255)                                                                                              gTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAc                                                                                                                                                <  1:320812/251‑1 (MQ=255)                                                                                                           gcTTTCCCTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGAGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTCTCACAGGCGAACGAGCGCAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTTTCGAATGATGGGGGGAAAAATATGAACGCGGCCCGCCAGCAAGCGATCGTAGCCCTGCTGCTGaaac                                                                                                                                   >  2:359216/1‑249 (MQ=255)                                                                                                                          cccGTGGCCTGACATCATACGAATCAATCGTTATGCGTTTTATGAACGGGGGCAAAAAGCCTTTGAGATAGTTATAACAGGCGAACGAGCGAAGCACGGGAATATTCTTTTAAAAAAAGAGGAAACAACGGAATCTCATACCGGAACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGGGAAAAATAAGAAAGCGGCACGCCAGAAAGGGATAGACGAAGTGCTGCGGAACCATAcacgccagaca                                                                                                                     >  2:198484/1‑240 (MQ=255)                                                                                                                                    gACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctc                                                                                                           <  1:6372/250‑1 (MQ=255)                                                                                                                                              gCATCAATCGATTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAAc                                                                                                  >  2:447331/1‑249 (MQ=255)                                                                                                                                                   aaTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTaaa                                                                                           <  2:66544/251‑1 (MQ=255)                                                                                                                                                      cGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGt                                                                                        <  1:41428/251‑1 (MQ=255)                                                                                                                                                          ttGCGTTTTATGAACGGCCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGCATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGCCGCGGGCGGTTATCGACTGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATCGTCGCCCTGCTGCTGAACCATACCAGACTGACCACGGAAGCTCTCTCTGAACCGCTAAAAGTCAGt                                                                                    >  2:266729/1‑251 (MQ=255)                                                                                                                                                           tGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGACGCTCTCTCTGAACAGCTAAAGGTCAGt                                                                                    >  1:272692/1‑250 (MQ=255)                                                                                                                                                             cGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCAGACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTaaa                                                                                 <  1:404239/251‑1 (MQ=255)                                                                                                                                                                       aCGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtc                                                                       <  1:176452/251‑1 (MQ=255)                                                                                                                                                                                                                         aCGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGAACGCCCGCATGACGCGGGCGGTTCCCGCATGATGGGGTGAAAAATATGAAAACGGCACGCCAGCAAGCGCTAGTCGACCTGCTGCTGCACCCTACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCAt                     >  1:466188/1‑251 (MQ=255)                                                                                                                                                                                                                            ggAATATTCTTTGAAAAAAAGTGGTAACACCGTAATCTCATTACGGTACGCCCGCATGACGCGTGCGGTGATCGAATGAGGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATCGGCGTCCTGCTGCTGATCCATCCCAGCCTGACCACGGGAGCTCTCTCGGAAAAGCTAAAGGGCAGTAAAGAAACCATTCGTCTCGATCTCAATGAATTGCAGACGCAGGGTAAAATCCTGCGCAATCATGGa                  <  2:394383/251‑1 (MQ=255)                                                                                                                                                                                                                                 atTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGAcgcg              <  1:359216/250‑1 (MQ=255)                                                                                                                                                                                                                                   tCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTAGCGTATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAATCTCGCTCGGAACAGCTAAAGGTCAGTAAAGAAACCATTCGGCGCGATCTCAGTGAATGACATAGGCATGTTGAAATGCTGTTCCATCATGGACGCGCt            <  2:242226/250‑1 (MQ=255)                                                                                                                                                                                                                                         aaaaaaaGTGGTAACACCGTAATCGTATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCCAGCGATAGACGACCTGCTGCTGAAGCATCCCAGCCTGACCACGGTAGCTCTCTCTGAACAGCTAAAGGACAGTAAAGAAACCATTCGTCGTGATCTCAATGAAGTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTAAatata     <  2:3920/251‑1 (MQ=255)                                                                                                                                                                                                                                           aaaaaGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTAAATATATc   >  1:169000/1‑251 (MQ=255)                                                                                                                                                                                                                                            aaaaGTGGTAACACCGTAATCTCATACAGGTACACCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGCACAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGATGAAACATACCAGCCTGACCACGGAAGATCTCTCTGAACAGCTAAATGTCAGAAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGCCACAGGGAAAAATTCTGCGCATTCATGGACGCGCTAAATATATcc  >  2:308343/1‑251 (MQ=255)                                                                                                                                                                                                                                             aaaGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTAAATATATcc  >  1:15264/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      TTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACGCGCTAAATATATCC  >  NC_000913/2939057‑2939542

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: ‑