Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,344,501	G→A	G378D (GGC→GAC)	melB →	melibiose:H(+)/Na(+)/Li(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,344,501	0	G	A	95.7%	60.1 / NA	23	G378D (GGC→GAC)	melB	melibiose:H(+)/Na(+)/Li(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/9);  minor base T (0/1);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TTATGGGCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGGCGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCC  >  NC_000913/4344260‑4344735                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           ttATGGGCCGGAGCTGCTTTTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTGCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAATCTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAg                                                                                                                                                                                                                                    <  2:72163/250‑1 (MQ=255)      ggCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTg                                                                                                                                                                                                                              <  1:249058/251‑1 (MQ=255)               aTCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGCGAAGGACGGTTCAGCCTTTGCGGCttttt                                                                                                                                                                                                                     >  1:64731/1‑251 (MQ=255)                          cGGTGTTTAGCTGTGGTGTTCTCCTGTTTATGGTATTAATGATCTATCACAACGTCGTCTTCATTGTGATTGCGGGTATTTTGCTGAATGTGGTAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACTGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGTCTTGTTTCATTGCGGtt                                                                                                                                                                                                          <  2:323498/251‑1 (MQ=255)                                                                aTGAGCTATCACAACGTCGTCCTCATTGAGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAAc                                                                                                                                                                     >  1:275129/1‑250 (MQ=255)                                                                    gCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCt                                                                                                                                                                >  1:19408/1‑251 (MQ=255)                                                                          acaACGTCGGCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCaa                                                                                                                                                          <  1:21900/251‑1 (MQ=255)                                                                                            tgATTGCGGGTATTTTGCTGAATGTAGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGGTATCGTTGATTACGGTGAATATAACCTGCCCGTACGATGTGAAAGTATCG‑CTTACTCCGTACAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTACGGTTGTGTTAGGGATGCTTGGCTATGTACCGAATGTTGTACAGTCTACTCAAGCCCTAtatggtctgc                                                                                                                                          >  2:355129/1‑240 (MQ=255)                                                                                                                                                    atcatGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACtt                                                                                >  1:65808/1‑251 (MQ=255)                                                                                                                                                             tcagttatCGTTGATTTCGTTGAATATAAACTGGCCGGACGTTGTGTAAG‑ATTGTCTTGCTCCGTGCAGATTTTGGTGGTGAAGGTCTGTTGAGCTTTTGCGGCTTTTTTCATTGCTGTTGTGTTAGGGATGATTGGCTATGTACCTAATGTTGAACAGCCTACTCAAGGCCTATTAGGTGTGCATTTTATTATGATTGCGCCATCAACTCTGTTTTTCATGGTAACGCTGATTCTCTTCTTCCGTTTCt                                                                        <  2:335869/245‑1 (MQ=255)                                                                                                                                                              cAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTAt                                                                      >  1:217733/1‑251 (MQ=255)                                                                                                                                                                gATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCg                                                                    >  1:124353/1‑251 (MQ=255)                                                                                                                                                                               ggTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTAAGGAGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTAACGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAAATCTGTTTTTCATGGTAACGCTGATTCCTTACTTCCGTTTCTATAGCCTCAATGGTGacac                                                     >  1:351683/1‑251 (MQ=255)                                                                                                                                                                                        aaaCTGCACGTCCGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGAAGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCTAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCATCTCGGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCa                                             <  2:19408/250‑1 (MQ=255)                                                                                                                                                                                         aaCTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAgg                                           >  1:292971/1‑251 (MQ=255)                                                                                                                                                                                              cACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGgatcc                                       >  1:267580/1‑250 (MQ=255)                                                                                                                                                                                                    cGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGATGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTTCCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCATTGGGGACACGCTGCGCAAgaaccagatcct                                >  1:201428/1‑250 (MQ=255)                                                                                                                                                                                                        gtgAAAATATCT‑CTTACTTCGGGCAGAATATGGTTGTGAAGGACGGTTCAGCCTGTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGTTATGTACCGAATGTTGAGCAGTTTACGCATGCCTTATTAGGTATTCAGTTTATTATGATTGCTCTTTCAACTCTGTTGTTCATGGTAACGCTGATTCTCTACTTCCGCTTCTATCGCCTCAATGGTGTCACGCTGCGTAGGATCCAGATCAATc                              <  2:64731/249‑1 (MQ=255)                                                                                                                                                                                                                aTCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGATAGGGATGATTGGCTATGTAACGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATa                     >  1:374092/1‑250 (MQ=255)                                                                                                                                                                                                                            cGTGCAGACTATGGTGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAg         <  1:190954/251‑1 (MQ=255)                                                                                                                                                                                                                                cAGACTATGGGGGTGAAGGACGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCACTATTAGGTATGCAGGTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGATGATTCTCTACTTGCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCCGATCAATCTGCTGGATAAATATCGCAAAGAAcc     >  2:72812/1‑251 (MQ=255)                                                                                                                                                                                                                                   aCTATGGTGGTGAAGGACGTTTCCGCCTTTGCGGATTTTTGCATTGAGGTTGTGTTAGGGATGATTGGCGATTTACCGAATGTTGAACAGTCTACGCACGCCCTAGTAGGTATGCAGTGTATTATTATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCCCTCCTTCCTTTTCTATCGCCTCAATTGTTACACGCTGCGCAGGATCCAGATTCATCTGATGGATAAATATCGCacattacctcc  <  2:301155/251‑11 (MQ=255)                                                                                                                                                                                                                                    gtttggtggTGAATGACGGTTCAGCCTTTGCGGTTGTGTTCATTGCGGTTTGGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGGAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTTCCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAAGGGTGACACGCTGCGCAGGATCCAGATCCTTCTGCTGGATAAATATCGCAAAGTACCGcc  <  2:342238/247‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           TTATGGGCCGGAGCATCTATTCTTCCGGTGTTAAGCTGTGGTGTTCTCCTGTTAATGGCATTAATGAGCTATCACAACGTCGTCCTCATTGTGATTGCGGGTATTTTGCTGAATGTGGGAACGGCGCTTTTCTGGGTATTACAGGTCATCATGGTGGCAGATATCGTTGATTACGGTGAATATAAACTGCACGTACGCTGTGAAAGTATCG‑CTTACTCCGTGCAGACTATGGTGGTGAAGGGCGGTTCAGCCTTTGCGGCTTTTTTCATTGCGGTTGTGTTAGGGATGATTGGCTATGTACCGAATGTTGAACAGTCTACGCAAGCCCTATTAGGTATGCAGTTTATTATGATTGCTCTACCAACTCTGTTTTTCATGGTAACGCTGATTCTCTACTTCCGTTTCTATCGCCTCAATGGTGACACGCTGCGCAGGATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCC  >  NC_000913/4344260‑4344735

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