Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,975,617	A→G	I133T (ATT→ACT)	aas ←	acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,975,617	0	A	G	100.0%	136.8 / NA	38	I133T (ATT→ACT)	aas	acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/21);  total (17/21)

GGGAACAAGCGACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAATAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGCCTGTCGTGGTGAT  >  NC_000913/2975542‑2975693                                                                            |                                                                             gggAACAAGCGACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAAcc                                                                         <  1:552354/81‑1 (MQ=255)     aCAAGCGACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTg                                                                     >  1:552417/1‑81 (MQ=255)     aCAAGCGACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTg                                                                     <  2:552355/81‑1 (MQ=255)        aGCGACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCAc                                                                  >  2:552411/1‑81 (MQ=255)        aGCGACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCAc                                                                  >  2:552399/1‑81 (MQ=255)           gACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCAc                                                                  >  2:552412/1‑78 (MQ=255)           gACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAg                                                               <  2:552356/81‑1 (MQ=255)            aCGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGa                                                              >  2:552420/1‑81 (MQ=255)                 tAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTc                                                          <  2:640633/80‑1 (MQ=255)                   aCCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTcgcc                                                       <  2:552358/81‑1 (MQ=255)                   aCCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTcgcc                                                       <  2:552357/81‑1 (MQ=255)                      aGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTcgccgc                                                     >  1:552413/1‑80 (MQ=255)                      aGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTcgccgc                                                     >  1:552429/1‑80 (MQ=255)                       gACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGa                                                   >  2:552410/1‑81 (MQ=255)                       gACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGa                                                   >  1:552418/1‑81 (MQ=255)                         ccTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACa                                                 <  2:552359/81‑1 (MQ=255)                          cTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACaa                                                <  1:552360/81‑1 (MQ=255)                             tCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACaaa                                               <  1:552361/79‑1 (MQ=255)                              cAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACACACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAAc                                             <  2:552362/80‑1 (MQ=255)                                 gcggcTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCcgc                                         >  2:552408/1‑81 (MQ=255)                                  cggcTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCcgcg                                        >  1:552414/1‑81 (MQ=255)                                      tGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCAt                                    >  1:552404/1‑81 (MQ=255)                                        aaGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATc                                   <  2:552363/80‑1 (MQ=255)                                        aaGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATc                                   <  2:552364/80‑1 (MQ=255)                                         aGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCg                                  >  1:552415/1‑80 (MQ=255)                                         aGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGt                                 <  2:552365/81‑1 (MQ=255)                                              gTAAGTTCCGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGAtt                            >  2:552428/1‑81 (MQ=255)                                                     ccGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTtcatc                      <  2:552367/80‑1 (MQ=255)                                                     ccGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTtcatc                      <  1:552366/80‑1 (MQ=255)                                                     ccGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTtcat                       >  2:552424/1‑79 (MQ=255)                                                     ccGCCCCTTCAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTtcat                       >  1:552427/1‑79 (MQ=255)                                                              cAATACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGcc              <  1:552368/79‑1 (MQ=255)                                                               aaTACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGCCTg            <  1:552369/80‑1 (MQ=255)                                                                aTACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGCCTGTc          >  2:552432/1‑81 (MQ=255)                                                                 tACGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGCCTGTCg         <  2:552370/81‑1 (MQ=255)                                                                   cGCACAGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGCCTGTCgtg       <  1:552371/81‑1 (MQ=255)                                                                        aGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGCCTGTCGTGGTGAt  <  2:552372/81‑1 (MQ=255)                                                                        aGGAGTAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGCCTGTCGTGGTGAt  <  2:552373/81‑1 (MQ=255)                                                                            |                                                                             GGGAACAAGCGACGTTTAACCAGACCTTTCAGGCGGCTGAAGTGCGTAAGTTCCGCCCCTTCAATACGCACAGGAATAACCGTTGCACCAGACTTCGCCGCGACAAAACCCGCGCCATCGTAGATTTTCATCAGCGAGCCTGTCGTGGTGAT  >  NC_000913/2975542‑2975693

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