Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	770,032	G→A	G685G (GGG→GGA)	mngB →	alpha‑mannosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	770,032	0	G	A	96.3%	75.1 / ‑6.3	27	G685G (GGG→GGA)	mngB	alpha‑mannosidase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (13/13);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGG  >  NC_000913/769789‑770266                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           ggATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTg                                                                                                                                                                                                                                     <  1:237876/251‑1 (MQ=255)     atCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAATGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAAtt                                                                                                                                                                                                                                 <  1:350369/251‑1 (MQ=255)       cAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTg                                                                                                                                                                                                                               >  2:272904/1‑251 (MQ=255)                   ccAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCACCTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGGGAAGGATTACGTGCATTTGAAGTCATTGGTg                                                                                                                                                                                                                   >  2:376019/1‑251 (MQ=255)                           aCGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTgaagagaag                                                                                                                                                                                                           <  1:110866/251‑1 (MQ=255)                             gATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGCTACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGAGAAGGATTAAGTGAATTTGACGTCATCGGTgaagagaaga                                                                                                                                                                                                          >  2:274388/1‑250 (MQ=255)                                 atcGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGATAAGTGAATTTGAAGTAATCGGGGAAGAGAAGAAAAcc                                                                                                                                                                                                     >  1:247610/1‑251 (MQ=255)                                   cGCCTTCTTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCAGTGGACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTTGAAAGAAGCGCCGGTTCCGGTATGGAATATGCGCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATTTCCGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCtt                                                                                                                                                                                                   <  1:280149/251‑1 (MQ=255)                                   cGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCtt                                                                                                                                                                                                   <  2:397327/251‑1 (MQ=255)                                                     ccTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTgcg                                                                                                                                                                                 >  1:244914/1‑251 (MQ=255)                                                          cccTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTgcgtggcg                                                                                                                                                                            <  2:133207/250‑1 (MQ=255)                                                            cTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTgcgtggcgtg                                                                                                                                                                          >  2:22104/1‑251 (MQ=255)                                                                     ccGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGCAAACCTTTGCCATTAAGTTTATGCGTGGCGGGGGCTTACTg                                                                                                                                                                 >  1:386355/1‑251 (MQ=255)                                                                             gtTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAaaga                                                                                                                                                         <  1:94838/251‑1 (MQ=255)                                                                                  ggCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAATAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGAt                                                                                                                                                     <  1:273580/250‑1 (MQ=255)                                                                                      gATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCt                                                                                                                                                <  1:272904/251‑1 (MQ=255)                                                                                           gCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTaa                                                                                                                                           <  1:22104/251‑1 (MQ=255)                                                                                                  ggTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATGTGAAGTCTTCGGTGATGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTg                                                                                                                                     <  2:399987/250‑1 (MQ=255)                                                                                                   gTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTggg                                                                                                                                   <  1:161913/251‑1 (MQ=255)                                                                                                                   cGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAAtt                                                                                                                    >  2:125434/1‑250 (MQ=255)                                                                                                                                        caacTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGGGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTTGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTcac                                                                                              >  2:97593/1‑251 (MQ=255)                                                                                                                                                                                     gAATATGCTCAACGATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGGAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTAAGTTGCTGCGTGGCGTGGGCTTAATGGGCAAAGAAGATCTGCTTTTAAGGACTGGGCGGCCTTCGGGAATTAAAAAGCCAGTCACGGACTCACAAATACGTGGTCTGCTTTATTGTCGCCTAAGTTTATTGAGGtat                                                  >  2:329710/1‑250 (MQ=255)                                                                                                                                                                                         atGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGCATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGGGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAACGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACCGGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCgg                                             >  2:126436/1‑251 (MQ=255)                                                                                                                                                                                                   atggtGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGTGAATTAAAATGCCAGTCCCGGACTCACAACTACGTTGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAAcc                                   <  2:292542/247‑1 (MQ=255)                                                                                                                                                                                                                    gggCGTAACGGCATGGCTGTCTTTAGGGAAGGATTACGTGAATTTGAAGGCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACCGGGCAAAGAAGATCTTCTTTTAAGTCCTGGGCGGCCTTCTGGAATTAAAATGCCAGTCCCGGACTCACAACTATGTGGTCTGCTTTCTTGTCGCCTCAGTTTTTTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTc                   >  2:240200/1‑250 (MQ=255)                                                                                                                                                                                                                    gggCGTAACGGCATGGCTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTc                   >  2:167629/1‑250 (MQ=255)                                                                                                                                                                                                                                    cTGTCTTTAGCGAAGGATTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCATAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCGGGTGTAGCTCAGCAGGCGCGAGCATgg  >  2:394955/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           GGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGGTCCCTACACCTTTTAACACCGACAGTGTTCTGGCAGATACGCAGTTTGGTTCGCTAACGCGCCCCGTGAACGACAGTGCAATGAACAACTGGCAGCAAGAAGGCTGGAAAGAAGCGCCGGTTCCGGTATGGAATATGCTCAACTATGTTGCCTTACAGGAAGGGCGTAACGGCATGGCTGTCTTTAGCGAAGGGTTACGTGAATTTGAAGTCATCGGTGAAGAGAAGAAAACCTTTGCCATTACGTTGCTGCGTGGCGTGGGCTTACTGGGCAAAGAAGATCTGCTTTTAAGGCCTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGG  >  NC_000913/769789‑770266

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