Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,589,575	A→G	F569L (TTC→CTC)	btsT ←	pyruvate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,589,575	0	A	G	100.0%	174.4 / NA	50	F569L (TTC→CTC)	btsT	pyruvate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/26);  total (24/26)

GTGGTGCAGATAAGCAGCCATACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAACAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCAACCAGAGAATCGGTTTTTTTCAGGAAC  >  NC_000913/4589445‑4589705                                                                                                                                   |                                                                                                                                   gtggtgCAGATAAGCAGCCATACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCca                                                                                                                           <  1:148252/140‑1 (MQ=255) gtggtgCAGATAAGCAGCCATACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCca                                                                                                                           <  1:403063/140‑1 (MQ=255)   ggtgCAGATAAGCAGCCATACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCcaca                                                                                                                         <  2:590267/140‑1 (MQ=255)          aTAAGCAGCCATACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCtt                                                                                                                  >  1:950215/1‑140 (MQ=255)            aaGCAGCCATACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTaa                                                                                                                >  2:918373/1‑140 (MQ=255)                  ccATACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAAcgccgc                                                                                                          >  1:269992/1‑140 (MQ=255)                    aTACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGccc                                                                                                        <  2:4918/140‑1 (MQ=255)                    aTACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGccc                                                                                                        >  1:727275/1‑140 (MQ=255)                          ccGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGa                                                                                                  <  1:993862/140‑1 (MQ=255)                           cGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGAt                                                                                                 <  2:70229/140‑1 (MQ=255)                           cGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGAt                                                                                                 <  2:832527/140‑1 (MQ=255)                            ggAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGAt                                                                                                 <  1:1084523/139‑1 (MQ=255)                            ggAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGAt                                                                                                 >  2:1084523/1‑139 (MQ=255)                            ggAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATc                                                                                                >  2:519255/1‑140 (MQ=255)                                    aGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGc                                                                                        <  1:1007834/140‑1 (MQ=255)                                     gTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGcc                                                                                       <  2:1028428/140‑1 (MQ=255)                                        aCCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGccc                                                                                      >  2:628682/1‑138 (MQ=255)                                        aCCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGccc                                                                                      <  1:628682/138‑1 (MQ=255)                                        aCCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTg                                                                                    <  1:893283/140‑1 (MQ=255)                                         cccAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTGGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGa                                                                                   <  1:59841/140‑1 (MQ=255)                                                     gggTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAg                                                                       <  2:107808/140‑1 (MQ=255)                                                     gggTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAg                                                                       <  1:616445/140‑1 (MQ=255)                                                      ggTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGc                                                                      >  2:1080167/1‑140 (MQ=255)                                                        tGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGccc                                                                    >  2:986422/1‑140 (MQ=255)                                                        tGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGccc                                                                    >  2:884612/1‑140 (MQ=255)                                                          cgcTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCca                                                                  >  1:162810/1‑140 (MQ=255)                                                            cTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCcaca                                                                <  1:519255/140‑1 (MQ=255)                                                              tCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGa                                                              <  2:269992/140‑1 (MQ=255)                                                                      aTCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCcacaca                                                      <  1:918373/140‑1 (MQ=255)                                                                       tCAGCACAACGGTGCCCTGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAg                                                     >  2:315117/1‑140 (MQ=255)                                                                       tCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAg                                                     >  2:1015374/1‑140 (MQ=255)                                                                         aGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGcc                                                   >  1:455892/1‑140 (MQ=255)                                                                           cacaACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCg                                                 >  2:611362/1‑140 (MQ=255)                                                                           cacaACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCg                                                 >  1:1013073/1‑140 (MQ=255)                                                                           cacaACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCg                                                 >  1:593472/1‑140 (MQ=255)                                                                             caACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCa                                               <  2:1004352/140‑1 (MQ=255)                                                                             caACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCa                                               <  2:950215/140‑1 (MQ=255)                                                                                cGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTa                                            >  2:1200702/1‑140 (MQ=255)                                                                                     cccAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCACCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCgat                                       >  2:657346/1‑140 (MQ=255)                                                                                              aGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGc                              <  2:202459/140‑1 (MQ=255)                                                                                              aGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGc                              <  1:452794/140‑1 (MQ=255)                                                                                              aGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGc                              <  1:315117/140‑1 (MQ=255)                                                                                               gcgcTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCa                             >  2:495540/1‑140 (MQ=255)                                                                                                cgcTACGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCaa                            >  1:1013695/1‑140 (MQ=255)                                                                                                    aCGGCTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCAACCag                        >  2:767006/1‑140 (MQ=255)                                                                                                        cTGCCAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCAACCAGAGaa                    <  1:1200702/140‑1 (MQ=255)                                                                                                           ccAGCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCAACCAGAGAATCg                 >  1:209255/1‑140 (MQ=255)                                                                                                              gCATCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGACAACCAGAGAATCGGtt              <  2:455892/140‑1 (MQ=255)                                                                                                                 tCTGGTTGGAGATACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCAACCAGAGAATCGGttttt           <  1:986422/140‑1 (MQ=255)                                                                                                                          agaTACCGAGCAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCAACCAGAGAATCGGTTTTTTTCAGGAAc  >  2:297513/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                   GTGGTGCAGATAAGCAGCCATACAGCCGGAACAACAGTTACCCAGATGTATTGGGTGCGCTTCATCTTAATCAGCACAACGGTGCCCAGTACCAGCGCTACGGCTGCCAGCATCTGGTTGGAGATACCGAACAGCGGCCACAGGCTCTTAACGCCGCCCAGCGGATCAACCACGCCCTGATACAGCAGGTAGCCCCACAGACCCACACAGCCCGCAGTACCGATGATACCGGCAACCAGAGAATCGGTTTTTTTCAGGAAC  >  NC_000913/4589445‑4589705

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