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	88.5%	51.2 / 2.6	27	G137V (GGG→GTG)	fucU	L‑fucose mutarotase
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (10/13);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

AGGCGATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑GGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACG  >  NC_000913/2939051‑2939528                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        aGGCGATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTa                                                                                                                                                                                                                                     <  1:278472/251‑1 (MQ=255)  ggCGATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTa                                                                                                                                                                                                                                     >  1:132242/1‑250 (MQ=255)  ggCGATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑GGGTa                                                                                                                                                                                                                                     >  1:78768/1‑250 (MQ=255)      aTTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACAcc                                                                                                                                                                                                                                >  1:35121/1‑251 (MQ=255)              gTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGGCGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGACCGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGGACGGGAATATTCTTTTAAAAAAAG‑TGGTACCAACGTAAtctc                                                                                                                                                                                                                        >  2:331148/1‑251 (MQ=255)                            aCAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCAATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGTGTTTTATGAACGGGCGCAAAAAGCCTTTTCGGACGTTATTCCATGCGAACGAGCGAAGGACGGGAAAATTCTTTTAAAAAAAGTGGGTA‑CCCAGGAATCTCATACCGGTAAGccc                                                                                                                                                                                                          >  2:169190/1‑251 (MQ=255)                                   tGCACAGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACCCTCTCGATCATGAAGTAGAACGACGTTACCGTAATGCGCTTTAACTAAAAGCCCCGTGTCCTGCCATCATCCGCATCCATCGTTTTTCGTTTTATACACGGGCGCAAAAAGACTTTGCGATCGTTATCACAGGCGAACGAGCGGAGAACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCTTACCGGGACGCCCGCATGAc                                                                                                                                                                                                   >  2:123683/1‑251 (MQ=255)                                               ggTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTAt                                                                                                                                                                                       >  2:75104/1‑251 (MQ=255)                                                      atgGCGGCGGAAGAAGGTGACCATCTCGATCCTGAAAGAGAAAGACGTTACCGTAATGCGCTTTCACTACAAGCCCCATGTCCTGACATCATCCGCATAACACGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGCGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑AGGTAACACCGTAATCCCATACCGGTAGGCCCGCATGaacaaggcgggttgactatg                                                                                                                                                                                 >  2:359407/1‑231 (MQ=255)                                                          cggcggTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                                                                            <  2:108708/251‑1 (MQ=255)                                                                  gAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGaaaa                                                                                                                                                                     >  1:45821/1‑250 (MQ=255)                                                                        gACACTCTCGATCCTGAATGAGAACGACGTTACCGTAATGCGCTTTTTCTACAAGCCCCGTGTCCTGACATCATCCGTATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCATGTTCGATCGTTATCACAGGCGAACGATCGAAGTACGGTAATATTCTTTTAAAAAAAG‑TCGTAACACCGTAATCTCATCCCGGTACGCCCTCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAatagga                                                                                                                                                               <  1:224245/250‑4 (MQ=255)                                                                         acacTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaaa                                                                                                                                                             <  2:146471/251‑1 (MQ=255)                                                                                cGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGCATGATGGGGTGAAAAATATGAAAGCGGCAc                                                                                                                                                      >  1:174269/1‑251 (MQ=255)                                                                                    ccTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCa                                                                                                                                                  <  2:166709/251‑1 (MQ=255)                                                                                    ccTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑GGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGGGGAAAATATGAAAACGGCACgac                                                                                                                                                   >  1:196754/1‑248 (MQ=255)                                                                                                             aTGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctgctg                                                                                                                         <  1:110720/251‑1 (MQ=255)                                                                                                                    ttCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGTTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCAt                                                                                                                   <  2:132242/250‑1 (MQ=255)                                                                                                                        cTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAcc                                                                                                                >  1:80536/1‑249 (MQ=255)                                                                                                                                                                       ttATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGGACACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGaa                                                                >  1:83275/1‑250 (MQ=255)                                                                                                                                                                         aTGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAAcc                                                             <  1:30794/251‑1 (MQ=255)                                                                                                                                                                             aCGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtc                                                         <  1:307265/251‑1 (MQ=255)                                                                                                                                                                                       aaaGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAtct                                               <  1:88095/251‑1 (MQ=255)                                                                                                                                                                                                  gatcgtgagcaacggCGAAGGAGCGAAGTACGGGAATATGCTTTTAAAAAAAG‑TGGTAACACCGTAATTTCATACCGGGACTCCCGAAGGTCGCGGGCGGTTATCGAATGATGGGGTGAAAAATAGGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCCCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTa                                     <  2:196754/237‑1 (MQ=255)                                                                                                                                                                                                                       aGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTgc               <  1:305284/251‑1 (MQ=255)                                                                                                                                                                                                                                cGGGAATATTCTTTTAAAAAAAG‑TGGTAACAACGTCATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTAACGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGATAAAGGTCAGTAAAGAAACCATTCGTCGCGATATCATTGAATTACCAACGAAGGGTACAATTCTGAGCAATCATg      >  2:372153/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaTATTCTTTTAAAAAAAG‑TGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGAcg  <  1:76715/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        AGGCGATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAG‑GGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGCGCAATCATGGACG  >  NC_000913/2939051‑2939528

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