Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,218,207	C→T	A117A (GCG→GCA)	yehY ←	putative ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,218,207	0	C	T	100.0%	191.1 / NA	59	A117A (GCG→GCA)	yehY	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/30);  total (29/30)

TGAGAGATCGTTAAGCGTGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGCGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCACA  >  NC_000913/2218075‑2218341                                                                                                                                     |                                                                                                                                       tGAGAGATCGTTAAGCGTGCCGGAGTACAGCAACCACAGCGAAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATTCGTGGAGATTCGGCGGATGACATCGCTACAGGCCAGCAATGCCAGTGCCGca                                                                                                                                  >  2:1365015/1‑139 (MQ=255)   agagaTCGTTAAGCGTGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGcagc                                                                                                                                >  1:913304/1‑139 (MQ=255)   agagaTCGTTAAGCGTGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGcagc                                                                                                                                >  1:1443202/1‑139 (MQ=255)       aTCGTTAAGCGTGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGcagccagc                                                                                                                            <  2:1443202/139‑1 (MQ=255)                  tGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAAc                                                                                                                     <  2:1304763/135‑1 (MQ=255)                  tGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAAc                                                                                                                     >  1:1304763/1‑135 (MQ=255)                  tGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACt                                                                                                                 >  1:365438/1‑139 (MQ=255)                  tGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACt                                                                                                                 >  1:548980/1‑139 (MQ=255)                  tGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACt                                                                                                                 >  1:1411460/1‑139 (MQ=255)                   gCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTg                                                                                                                >  1:1452107/1‑139 (MQ=255)                        aGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAg                                                                                                           >  1:14974/1‑139 (MQ=255)                          tACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGc                                                                                                         <  1:1260192/139‑1 (MQ=255)                                 aCCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTAcg                                                                                                  <  2:1452107/139‑1 (MQ=255)                                    acaGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGaaa                                                                                                                      <  1:410734/116‑1 (MQ=255)                                    acaGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGaaa                                                                                                                      >  2:410734/1‑116 (MQ=255)                                         ggAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAgcg                                                                                          >  1:502425/1‑139 (MQ=255)                                         ggAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAgcg                                                                                          <  1:319689/139‑1 (MQ=255)                                         ggAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAAGCTGGTACGCGCCAgcg                                                                                          <  1:944204/139‑1 (MQ=255)                                                             atGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGc                                                                                                         <  2:329972/104‑1 (MQ=255)                                                             atGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGc                                                                                                         >  1:329972/1‑104 (MQ=255)                                                              tGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCGAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGtt                                                                     <  2:913304/139‑1 (MQ=255)                                                                 aacaacCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGACATGGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCAATTTTGCGCCCATTggg                                                                  >  2:1017048/1‑139 (MQ=255)                                                                  acaacCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGt                                                                 <  2:14974/139‑1 (MQ=255)                                                                    aacCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGcagcc                                                                                                                               >  2:1439237/1‑75 (MQ=255)                                                                    aacCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGcagcc                                                                                                                               <  1:1439237/75‑1 (MQ=255)                                                                        aGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTc                                                                <  1:1127536/134‑1 (MQ=255)                                                                        aGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTc                                                                <  1:1296466/134‑1 (MQ=255)                                                                        aGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTc                                                                >  2:1127536/1‑134 (MQ=255)                                                                        aGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTc                                                                >  2:1296466/1‑134 (MQ=255)                                                                         gcgcCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCgccgcc                                                          >  1:1234885/1‑139 (MQ=255)                                                                          cgcCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCt                                                         <  2:237041/139‑1 (MQ=255)                                                                           gcCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCG‑CAGTTGGGTCGCCGCCttt                                                       >  2:208310/1‑139 (MQ=255)                                                                              acaGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTcc                                                     <  2:548980/139‑1 (MQ=255)                                                                              acaGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTcc                                                     <  2:1234885/139‑1 (MQ=255)                                                                               caGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCa                                                    >  1:52182/1‑139 (MQ=255)                                                                                    gATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTgc                                                                            <  2:1123571/110‑1 (MQ=255)                                                                                    gATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTgc                                                                            >  1:1123571/1‑110 (MQ=255)                                                                                     aTGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCa                                                                                                                   <  2:1215745/70‑1 (MQ=255)                                                                                     aTGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCa                                                                                                                   >  1:1215745/1‑70 (MQ=255)                                                                                                   ggcggATGGCATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATAc                                <  2:86591/139‑1 (MQ=255)                                                                                                        aTGGCATCGCTACAGGCCAGCAATGCCAGTGCCGcagcc                                                                                                                               <  2:1246095/39‑1 (MQ=38)                                                                                                        aTGGCATCGCTACAGGCCAGCAATGCCAGTGCCGcagcc                                                                                                                               >  1:1246095/1‑39 (MQ=38)                                                                                                          ggCATCGCTACAGGCCAGCAATGCCAGTGCCGcagc                                                                                                                                <  2:1180823/36‑1 (MQ=38)                                                                                                          ggCATCGCTACAGGCCAGCAATGCCAGTGCCGcagc                                                                                                                                >  1:1180823/1‑36 (MQ=255)                                                                                                            cATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGaa                       <  1:1017048/139‑1 (MQ=255)                                                                                                            cATCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGaa                       <  2:502425/139‑1 (MQ=255)                                                                                                              tCGCTACAGGCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAAtt                     <  2:365438/139‑1 (MQ=255)                                                                                                                       gCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACt                                                                                                                 <  1:163301/38‑1 (MQ=255)                                                                                                                       gCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACt                                                                                                                 >  2:163301/1‑38 (MQ=255)                                                                                                                       gCCAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCG‑CAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAAt           <  1:208310/139‑1 (MQ=255)                                                                                                                         cAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGc                                                                                    <  2:729579/65‑1 (MQ=255)                                                                                                                         cAGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGc                                                                                    >  1:729579/1‑65 (MQ=255)                                                                                                                          aGCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGa         >  2:339479/1‑139 (MQ=255)                                                                                                                           gCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTccc                                              >  1:676302/1‑101 (MQ=255)                                                                                                                           gCAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTccc                                              <  2:676302/101‑1 (MQ=255)                                                                                                                            cAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAATTgc                   <  2:767241/127‑1 (MQ=255)                                                                                                                            cAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAATTgc                   >  1:767241/1‑127 (MQ=255)                                                                                                                            cAATGCCAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGt       >  2:641164/1‑139 (MQ=255)                                                                                                                                 ccAGTGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGcaca  <  1:641164/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       TGAGAGATCGTTAAGCGTGCCGGAGTACAGCAACCACAGCGGAATAATGGCAATCTGCATATGCAACAACCAGCGCCACAGCGGATGCGTGGAGATTCGGCGGATGGCATCGCTACAGGCCAGCAATGCCAGCGCCGCAGCCAGCCAGAAACCACTGCCGAGGCTGGTACGCGCCAGCGCACTGCCATTTTGCGCCAGTTGGGTCGCCGCCTTTCCAGCTCCCCACACCAGCAATACGAAGACGAATTGCGCCAGAATGAGTGCACA  >  NC_000913/2218075‑2218341

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