Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,451,238	G→A	A61T (GCC→ACC)	ytfR →	galactofuranose ABC transporter putative ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,451,238	0	G	A	100.0%	88.1 / NA	27	A61T (GCC→ACC)	ytfR	galactofuranose ABC transporter putative ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/17);  total (10/17)

AAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGATAATC  >  NC_000913/4451102‑4451370                                                                                                                                         |                                                                                                                                     aaGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCg                                                                                                                                   >  1:333401/1‑140 (MQ=255)    tAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATc                                                                                                                                <  1:125977/140‑1 (MQ=255)       aTTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATcgcg                                                                                                                             >  1:106876/1‑140 (MQ=255)                cGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATct                                                                                                                    >  2:201846/1‑140 (MQ=255)                   cGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATctg                                                                                                                   <  2:45307/138‑1 (MQ=255)                   cGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATctg                                                                                                                   >  1:45307/1‑138 (MQ=255)                      cAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATctggctgg                                                                                                              <  1:322909/140‑1 (MQ=255)                       aaaGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGa                                                                                                             <  2:139432/140‑1 (MQ=255)                           cGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGc                                                                                                         <  1:201846/140‑1 (MQ=255)                               aGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAgg                                                                                                     >  2:261270/1‑140 (MQ=255)                                          aTTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCg                                                                                          >  1:193070/1‑140 (MQ=255)                                                        cGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCa                                                                            <  1:165014/140‑1 (MQ=255)                                                             cGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACgcgc                                                                       <  1:261270/140‑1 (MQ=255)                                                                   cATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCaacaac                                                                 <  2:333401/140‑1 (MQ=255)                                                                             ctCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTgg                                                       <  2:227010/140‑1 (MQ=255)                                                                                ggTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCAc                                                    <  2:106876/140‑1 (MQ=255)                                                                                     aaaCGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCt                                               >  2:63617/1‑140 (MQ=255)                                                                                             cGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGaa                                       <  2:193070/140‑1 (MQ=255)                                                                                                           cTAATCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTAcc                         <  1:65577/140‑1 (MQ=255)                                                                                                               tCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAat                     <  1:63617/140‑1 (MQ=255)                                                                                                               tCAAAGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAat                     >  2:36184/1‑140 (MQ=255)                                                                                                                   aGCATTAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATgt                 >  2:275080/1‑140 (MQ=255)                                                                                                                       ttAACTGGTGTATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATgtcggt             <  1:107417/140‑1 (MQ=255)                                                                                                                              gtgtATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGAt      <  1:356996/140‑1 (MQ=255)                                                                                                                              gtgtATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGAt      <  1:270059/140‑1 (MQ=255)                                                                                                                               tgtATACCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGATa     >  1:39459/1‑140 (MQ=255)                                                                                                                                  ataCCACACCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGATAATc  <  1:275080/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                     AAGTAAATTTTTCCCCGGCGTCAAAGCGTTAGACAACGTTGATTTCAGCCTGCGCCGTGGCGAAATCATGGCGCTGCTCGGTGAAAACGGGGCGGGAAAATCAACGCTAATCAAAGCATTAACTGGTGTATACCACGCCGATCGCGGCACCATCTGGCTGGAAGGCCAGGCTATCTCACCGAAAAATACCGCCCACGCGCAACAACTCGGCATTGGCACCGTCTATCAGGAAGTCAACCTGCTACCCAATATGTCGGTCGCTGATAATC  >  NC_000913/4451102‑4451370

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