Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,472,207	G→A	G436G (GGC→GGT)	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
*	NC_000913	3,472,207	0	G	A	100.0%	198.6 / NA	61	G436G (GGC→GGT)	fusA	protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/33);  total (28/33)

TCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGGCCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATCCGGGTCACACAGGGTGTCACCA  >  NC_000913/3472071‑3472336                                                                                                                                         |                                                                                                                                  tCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGAcc                                                                                                                                 <  2:1028812/139‑1 (MQ=255)    tttCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGACGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCATTAGCCAGACGACCCag                                                                                                                              >  1:1040579/2‑139 (MQ=255)       cAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGc                                                                                                                           <  1:596066/139‑1 (MQ=255)               aaTTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGAc                                                                                                                      <  2:343783/136‑1 (MQ=255)               aaTTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGAc                                                                                                                      >  1:343783/1‑136 (MQ=255)                           aTACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGt                                                                                                       >  1:173238/1‑139 (MQ=255)                                gTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGGCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCt                                                                                                  >  1:910650/1‑139 (MQ=255)                                           tgTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCgg                                                                                       >  2:788196/1‑139 (MQ=255)                                            gTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGt                                                                                      <  2:599699/139‑1 (MQ=255)                                                aGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCaa                                                                                  <  2:1040579/139‑1 (MQ=255)                                                 ggTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAAc                                                                                 <  1:788196/139‑1 (MQ=255)                                                       gTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGAcc                                                                                                                     <  1:29196/97‑1 (MQ=255)                                                       gTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGAcc                                                                                                                     >  2:29196/1‑97 (MQ=255)                                                          cGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGa                                                                            <  1:241482/135‑1 (MQ=255)                                                          cGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGa                                                                            >  2:241482/1‑135 (MQ=255)                                                               aTACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTa                                                                   <  2:173238/139‑1 (MQ=255)                                                                tACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTAc                                                                  <  1:411568/139‑1 (MQ=255)                                                                    cgcgATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGtt                                                                                            >  2:1510642/1‑109 (MQ=255)                                                                    cgcgATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGtt                                                                                            <  1:1510642/109‑1 (MQ=255)                                                                             ggTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAAt                                                     >  2:254961/1‑139 (MQ=255)                                                                               tCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTc                                                   >  2:1475051/1‑139 (MQ=255)                                                                                  ggTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCAt                                                <  1:333696/139‑1 (MQ=255)                                                                                   gTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTgg                                                                                                        >  1:644477/1‑82 (MQ=255)                                                                                   gTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTgg                                                                                                        <  2:327527/82‑1 (MQ=255)                                                                                   gTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTgg                                                                                                        >  1:327527/1‑82 (MQ=255)                                                                                   gTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTgg                                                                                                        <  2:644477/82‑1 (MQ=255)                                                                                   gTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATa                                               <  2:910650/139‑1 (MQ=255)                                                                                         ttcttcGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGAccc                                                                                                                                >  2:277437/1‑52 (MQ=255)                                                                                         ttcttcGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGAccc                                                                                                                                <  1:277437/52‑1 (MQ=255)                                                                                             tcGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGaa                                     >  1:447569/1‑139 (MQ=255)                                                                                               gTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCTGCTTTGGtttt                                                                                          <  1:621238/84‑1 (MQ=255)                                                                                               gTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCTGCTTTGGtttt                                                                                          >  2:621238/1‑84 (MQ=255)                                                                                                    tCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTAc                                                                  <  1:1563711/103‑1 (MQ=255)                                                                                                    tCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTAc                                                                  >  2:1563711/1‑103 (MQ=255)                                                                                                    tCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAgg                                                         <  2:242675/112‑1 (MQ=255)                                                                                                    tCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAgg                                                         >  1:242675/1‑112 (MQ=255)                                                                                                    tCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATg                                   <  1:545218/134‑1 (MQ=255)                                                                                                    tCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATg                                   >  2:545218/1‑134 (MQ=255)                                                                                                         acacGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCgg                              >  2:926804/1‑134 (MQ=255)                                                                                                         acacGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCgg                              <  1:926804/134‑1 (MQ=255)                                                                                                          cacGGAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATc                        <  2:447569/139‑1 (MQ=255)                                                                                                             ggAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATCCgg                     >  1:227855/1‑139 (MQ=255)                                                                                                              gAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCt                                                                                                  >  1:397298/1‑61 (MQ=255)                                                                                                              gAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCt                                                                                                  <  2:397298/61‑1 (MQ=255)                                                                                                              gAAAGACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATCCggg                    >  1:503759/1‑139 (MQ=255)                                                                                                                  gACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGa                                                                            <  2:949544/79‑1 (MQ=255)                                                                                                                  gACGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGa                                                                            >  1:949544/1‑79 (MQ=255)                                                                                                                   aCGGGTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATCCGGGTcaca               >  1:871930/1‑139 (MQ=255)                                                                                                                       gTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTc                                                                                    >  2:917787/1‑66 (MQ=255)                                                                                                                       gTCTTCTTTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTc                                                                                    <  1:917787/66‑1 (MQ=255)                                                                                                                         cttcttTAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATCCGGGTCACACAGGgt         <  1:254961/139‑1 (MQ=255)                                                                                                                               tAGCCGGACGACCCAGACTCAGACCCATTGTTTCCTGGTCAGCTTTGTTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCAACCGGGTCACGCAGGGTGTCAcc   <  2:1371597/139‑1 (MQ=255)                                                                                                                               tAGCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATCCGGGTCACACAGGGTGTCAcc   <  2:871930/139‑1 (MQ=255)                                                                                                                                 gCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATgg                                                                         <  1:1377160/67‑1 (MQ=255)                                                                                                                                 gCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATgg                                                                         >  2:1377160/1‑67 (MQ=255)                                                                                                                                 gCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAggg                                                        >  1:1120258/1‑84 (MQ=255)                                                                                                                                 gCCAGACGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAggg                                                        <  2:1120258/84‑1 (MQ=255)                                                                                                                                      aCGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGtt                                                                                            <  2:238594/43‑1 (MQ=255)                                                                                                                                      aCGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCAt                         <  2:148970/110‑1 (MQ=255)                                                                                                                                      aCGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCAt                         >  1:148970/1‑110 (MQ=255)                                                                                                                                      aCGACCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATCCGGGTCACACAGGGTGTCACCa  <  2:1028488/133‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  TCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGATTCTTCGTCAGTCCATACACGGAAAGACGGGTCTTCTTTAGCCAGACGGCCCAGAGCCAGACCCATTTTTTCCTGGTCAGCTTTGGTTTTCGGTTCAACTGCGATGGAGATTACCGGCTCAGGGAATTCCATACGTTCCAGAATGATCGGCGCATCCGGGTCACACAGGGTGTCACCA  >  NC_000913/3472071‑3472336

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: ‑