Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,670,581	G→A	G295R (GGA→AGA)	ynfM →	putative arabinose efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,670,581	0	G	A	100.0%	130.9 / NA	41	G295R (GGA→AGA)	ynfM	putative arabinose efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/22);  total (19/22)

GTTCGTCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCGGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTACTCTTCT  >  NC_000913/1670457‑1670707                                                                                                                             |                                                                                                                               gttcgtCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGaccacc                                                                                                                  <  2:354356/139‑1 (MQ=255)     gtCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCt                                                                                                              <  2:84797/139‑1 (MQ=255)     gtCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCt                                                                                                              <  2:300584/139‑1 (MQ=255)                    tACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAg                                                                                               <  1:442627/139‑1 (MQ=255)                        tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  1:154529/1‑127 (MQ=255)                        tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  2:154529/127‑1 (MQ=255)                        tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGAt                                                                                           >  2:377084/1‑139 (MQ=255)                        tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGAt                                                                                           >  2:263148/1‑139 (MQ=255)                        tCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGAt                                                                                           >  1:764854/1‑139 (MQ=255)                              aTCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATgttgtt                                                                                     <  2:205044/139‑1 (MQ=255)                                    tGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               <  2:190239/139‑1 (MQ=255)                                     gATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  1:158622/1‑114 (MQ=255)                                     gATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  2:158622/114‑1 (MQ=255)                                            tcCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTa                                                                       <  2:602847/139‑1 (MQ=255)                                            tcCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTa                                                                       >  2:281960/1‑139 (MQ=255)                                                   agCATGTCAGTCAGGCCGTGGTTGGCTTTCTATCGCTTGCTTCTTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCCCCCGCTATGCGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTCTGCTGtt                                                                >  2:566739/2‑139 (MQ=255)                                                        gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGTACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGttt                                                                                    <  1:389170/114‑1 (MQ=255)                                                        gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGTACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGttt                                                                                    >  2:389170/1‑114 (MQ=255)                                                        gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  1:642814/1‑95 (MQ=255)                                                        gtcagtcaGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  2:642814/95‑1 (MQ=255)                                                              cagGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACtga                                                     <  1:602050/139‑1 (MQ=255)                                                                     tGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               >  2:484826/1‑106 (MQ=255)                                                                     tGGATGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               <  1:484826/106‑1 (MQ=255)                                                                           gCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAg                                        <  1:281960/139‑1 (MQ=255)                                                                           gCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAg                                        >  1:23900/1‑139 (MQ=255)                                                                                         gCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               >  1:738024/1‑86 (MQ=255)                                                                                         gCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               <  2:738024/86‑1 (MQ=255)                                                                                              tttGACCGGTACATGGAGCTCACCCAAAGCCAGAACTATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACggggg                                                                          >  1:757493/1‑86 (MQ=255)                                                                                              tttGACCGGTACATGGAGCTCACCCAAAGCCAGAACTATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACggggg                                                                          <  2:757493/86‑1 (MQ=255)                                                                                                 gACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       <  1:625799/54‑1 (MQ=255)                                                                                                 gACCGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCg                                                                                                       >  2:625799/1‑54 (MQ=255)                                                                                                   ccGGTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCgg                >  2:299116/1‑139 (MQ=255)                                                                                                     ggTACATGGAGCTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTTATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGaa              <  2:467314/139‑1 (MQ=255)                                                                                                                cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               >  1:338276/1‑63 (MQ=255)                                                                                                                cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGAc                                                                               <  2:338276/63‑1 (MQ=255)                                                                                                                cTCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTACtcttc   >  1:148923/1‑139 (MQ=255)                                                                                                                 tCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGAc                                                 >  2:663518/1‑92 (MQ=255)                                                                                                                 tCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGAc                                                 <  1:663518/92‑1 (MQ=255)                                                                                                                 tCACCCAAAGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTACtcttct  <  2:764854/139‑1 (MQ=255)                                                                                                                        aaGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACtgatga                                                  >  2:639536/1‑84 (MQ=255)                                                                                                                        aaGCCAGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACtgatga                                                  <  1:639536/84‑1 (MQ=255)                                                                                                                             |                                                                                                                               GTTCGTCACGCTGTTTAATTACATCGGCTATCGGTTGATGCTCTCCCCCTGGCATGTCAGTCAGGCCGTGGTTGGCTTATTATCGCTGGCTTATTTGACCGGTACATGGAGCTCACCCAAAGCCGGAACCATGACCACCCGCTATGGGCGTGGTCCAGTGATGTTGTTTTCGACGGGGGTTATGCTGTTTGGTTTACTGATGACCTTATTCAGCTCGCTGTGGCTGATCTTTGCCGGAATGTTACTCTTCT  >  NC_000913/1670457‑1670707

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