Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	383,878	T→C	T20A (ACG→GCG)	yaiP ←	putative glucosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	383,878	0	T	C	100.0%	104.0 / NA	31	T20A (ACG→GCG)	yaiP	putative glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/22);  total (9/22)

GCGCACCACCTTTACCGGTATTTTTTTGCGTGACGGCAACAAAGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTT  >  NC_000913/383635‑384111                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          gcgcACCACCTTTACCGGTATTTTTTTGCGTGACGGCAACAAAGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATAACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTaa                                                                                                                                                                                                                                     >  2:185566/1‑250 (MQ=255)  cgccccaccTTTACCGGTTTTTTGTGGCGTGACGGCAACAAAGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTaaaa                                                                                                                                                                                                                                   <  2:237023/247‑1 (MQ=255)         ccTTTACCGGTATTTTTTTGCGTGACGGCAACAAAGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAg                                                                                                                                                                                                                            <  1:286990/251‑1 (MQ=255)                                gACGGCAACAAAGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                                                                                     <  1:279261/251‑1 (MQ=255)                                 aCGGAAACAAAGCGGGCGCCCCATTTGCGTGTGACATCCGCCATGACCGCTTCGGGATTGTCCGTGGAGCCGTCGTTAAAGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGGGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTGCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACACACAGAt                                                                                                                                                                                                    <  2:290884/251‑1 (MQ=255)                                  cGGCAACAAAGCGGTCGCCCCATTTGCTTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                                                                                                                                   <  2:305419/251‑1 (MQ=255)                                   ggCAACAAAGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                                                                                                                                   >  1:96002/1‑250 (MQ=255)                                         aaaGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                                                                                                                                            <  1:111848/251‑1 (MQ=255)                                          aaGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATGATCCAGTGACTGCGCCAGACACGGGCCGTCGTTATACGCAGGTATAAGGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCGTAGCGCACGTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCa                                                                                                                                                                                           <  1:334643/251‑1 (MQ=255)                                               gTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                                                                                                                                      <  1:287572/251‑1 (MQ=255)                                                                 gACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATGACCCGGCAAAAATAAGGGGTACGCAGTAGATTATCCAGTGACGGCGCCAGACACGGGCCTTCGGTATACGCAGGTATAATGGCGTCAATACAGCCTGTCTGTGGACTTGGTGTACGTCTTAGCGCACTGAAAGACCATAGCAACATTGCTATGGACATACGGTTAAAGATCCAGGGTTGCATGGAAGTTGCTCaaa                                                                                                                                                                    <  2:18232/251‑1 (MQ=255)                                                                     tCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGaa                                                                                                                                                                 <  1:114910/250‑1 (MQ=255)                                                                               ccGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGAGTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTGTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGa                                                                                                                                                      <  2:267268/251‑1 (MQ=255)                                                                                      ggTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACa                                                                                                                                                >  1:324740/1‑250 (MQ=255)                                                                                      ggTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACa                                                                                                                                                >  2:328108/1‑250 (MQ=255)                                                                                         aTTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAg                                                                                                                                             <  2:327471/250‑1 (MQ=255)                                                                                              ccGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCa                                                                                                                                       <  2:336935/251‑1 (MQ=255)                                                                                                 tGGAGCCGTCGATAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                                                                                                                                    >  2:235754/1‑251 (MQ=255)                                                                                                                                     aaGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGc                                                                                                 <  2:112634/250‑1 (MQ=255)                                                                                                                                                aGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGc                                                                                       >  1:209904/1‑249 (MQ=255)                                                                                                                                                    gATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGa                                                                                 <  1:328108/251‑1 (MQ=255)                                                                                                                                                                       cAGACACGGGCCGTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATGTTATGAATTACAAAGCCTTCACCCAGATCGCGATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATcgc                                                              <  2:38554/251‑1 (MQ=255)                                                                                                                                                                                gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCt                                                     >  2:95387/1‑251 (MQ=255)                                                                                                                                                                                 ccTTCGTTATACTCAGGTATAATTGCGTCTATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATACAGGTTTTCATGCACTTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTTTGCAACTGCACGCAGGCGATCATGACGCATATAATCTCGTCGTTCCTa                                                    <  1:287887/251‑1 (MQ=255)                                                                                                                                                                                         ataCTCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTTTCCGCTAAATTGTTCAACTGCACGCAGGCGATCATTACGCATATAATCGCGTCGTTCCTAGCTTTTa                                             <  2:58799/250‑1 (MQ=255)                                                                                                                                                                                          taCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATg                                           >  1:100634/1‑251 (MQ=255)                                                                                                                                                                                                   ataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTaa                                  <  1:148683/251‑1 (MQ=255)                                                                                                                                                                                                        ggCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTGTTTAATgca                              <  1:296499/250‑1 (MQ=255)                                                                                                                                                                                                              aaTACATCCTGTCTTTTGACTGGGCTTACGTCTTAGCGCACTTAAAAACCATAGCAACATAGCTATGGACATACCGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAACTTGTGCAACTGCACGCAGGCGATCATGAAGCTTATAATCGCGTCGTTCATAGCTTTTATGTTTGTTTAATGCAGCAAGtt                       >  2:109363/1‑251 (MQ=255)                                                                                                                                                                                                                               gACTGGGTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATAttct      <  1:258197/251‑1 (MQ=255)                                                                                                                                                                                                                                   gggTTTACGTCTTAGCGCACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTcttctt  <  1:312157/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          GCGCACCACCTTTACCGGTATTTTTTTGCGTGACGGCAACAAAGCGGTCGCCCCATTTGCGTTTGACTTCCGCCATGACCGCTTCGGTATTGTCCGTGGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATAATCGCGTCGTTCCTAGCTTTTATGTTTTTTTAATGCAGCAAGTTTGACTGTGAAATATTCTTCTT  >  NC_000913/383635‑384111

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