Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	3,466,500	G→A	T125I (ACC→ATC)	fusA ←	protein chain elongation factor EF‑G, GTP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,466,500	0	G	A	100.0%	122.2 / NA	36	T125I (ACC→ATC)	fusA	protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/16);  total (20/16)

TGCGGTCCATTTTGTTAACGAACGCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGGTTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTACCGCACCATCGAGAACA  >  CP009273/3466435‑3466562                                                                  |                                                               tGCGGTCCATTTTGTTAACGAACGCAATCCGCGGATCTTTATTTTTGTTTTCCTGACGCCATACGAtt                                                              <  1:4664684/68‑1 (MQ=255)      tCCATTTTGTTAACGAACGCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGa                                                         >  1:2847502/1‑68 (MQ=255)       ccATTTTGTTAACGAACGCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGAc                                                        >  1:1201312/1‑68 (MQ=255)            ttGTTAACGAACGCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCgg                                                   >  1:1087855/1‑68 (MQ=255)            ttGTTAACGAACGCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCgg                                                   >  1:1105324/1‑68 (MQ=255)              gTTAACGAACGCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCt                                                 >  1:1204138/1‑68 (MQ=255)              gTTAACGAACGCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGAATTCAGACTGCGGCt                                                 >  1:4744042/1‑68 (MQ=255)                    gaacgCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAacac                                           <  1:4732011/68‑1 (MQ=255)                        gCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAAcacctgt                                       >  1:3239636/1‑65 (MQ=255)                        gCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACaccacc                                       >  1:3898275/1‑68 (MQ=255)                        gCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACaccacc                                       >  1:3445223/1‑68 (MQ=255)                         cAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACAccacca                                      <  1:545739/68‑1 (MQ=255)                         cAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACAccacca                                      <  1:315270/68‑1 (MQ=255)                            tGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAAcc                                   >  1:9810/1‑67 (MQ=255)                                 ggATCTTTTTATTTGATTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACGAACTGCGCa                              <  1:3017893/68‑1 (MQ=255)                                 ggAACTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCa                              <  1:1063195/68‑1 (MQ=255)                                   aaCTTTATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGt                            <  1:4051777/68‑1 (MQ=255)                                      tttATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTaaa                         <  1:2692247/68‑1 (MQ=255)                                      tttATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTaaa                         <  1:2233650/68‑1 (MQ=255)                                       ttATATTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAAc                        <  1:2729969/68‑1 (MQ=255)                                       ttATATTTGTTAGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAAc                        <  1:1527546/68‑1 (MQ=255)                                        tataTTTGTTTGCCTGACGCCATACGATTTCAGCCTTCGGCTGAACACCACCAACTGCGCAGTAAAcc                       <  1:2669796/68‑1 (MQ=255)                                         atatTTTTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGTACACCACCAACTGCGCAGTAAACCa                      <  1:1871110/68‑1 (MQ=255)                                           atTTGTTTGCCTGACGCCATACGATTTCAGACTGCTGCTGAACACCACCAACTGCGCAGTAAACCAtt                    >  1:3176308/1‑68 (MQ=255)                                           atTTGTTTGCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCAtt                    >  1:4684657/1‑68 (MQ=255)                                              tgtttgCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTAcc                 >  1:2653296/1‑68 (MQ=255)                                              tgtttgCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTAcc                 >  1:1288982/1‑68 (MQ=255)                                              tgtttgCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTAcc                 >  1:783995/1‑68 (MQ=255)                                              tgtttgCCTGACGCAATACGATTTCAGACTGCGGCTGAACCCCACCCCCTGCGCAGTAAACCATGAcc                 >  1:4295522/1‑68 (MQ=255)                                                tttgCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTACCGc               >  1:4413006/1‑68 (MQ=255)                                                tttgCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTACCGc               <  1:1726811/68‑1 (MQ=255)                                                tttcCCTGACGCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTACCGc               <  1:1996844/68‑1 (MQ=255)                                                          gCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTACCGCACCATCgaga     >  1:3619263/1‑68 (MQ=255)                                                          gCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTACCGCACCATCgaga     >  1:998084/1‑68 (MQ=255)                                                          gCCATACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGT‑AACCATTACCGCACCATCGAGaa    >  1:1424260/1‑68 (MQ=255)                                                             aTACGATTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTACCGCACCATCGAGAaca  <  1:1116782/68‑1 (MQ=255)                                                                  |                                                               TGCGGTCCATTTTGTTAACGAACGCAATGCGCGGAACTTTATATTTGTTTGCCTGACGCCATACGGTTTCAGACTGCGGCTGAACACCACCAACTGCGCAGTAAACCATTACCGCACCATCGAGAACA  >  CP009273/3466435‑3466562

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