Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	316,398	A→G	H271R (CAC→CGC)	insF‑1 →	IS3 element protein InsF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	316,398	0	A	G	100.0%	158.2 / NA	46	H271R (CAC→CGC)	insF‑1	IS3 element protein InsF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/24);  total (22/24)

GGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCA  >  NC_000913/316347‑316531                                                    |                                                                                                                                      gggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGTGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaa                                               >  2:1135973/1‑140 (MQ=17)  ggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaaa                                              >  2:38198/1‑140 (MQ=18)   gCAACGGTGTTTAATTATATCGAATGTGATTACAGTCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTGCATATTACGTGGGTAGGATCATAAAg                                             >  2:1470295/1‑140 (MQ=255)   gCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAg                                             >  2:678730/1‑140 (MQ=18)         gtgtTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCg                                       <  1:1603537/140‑1 (MQ=255)              tAATTATATCGAATGTGATTACAATCGGTGGGGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                  <  2:1628995/140‑1 (MQ=255)              tAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                  <  1:951514/140‑1 (MQ=255)              tAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                  <  1:179449/140‑1 (MQ=255)                  tataTCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGt                              >  2:1618094/1‑140 (MQ=255)                    tatCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTcc                            >  1:1110223/1‑140 (MQ=255)                    tatCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTcc                            >  2:1474310/1‑140 (MQ=255)                     atCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                           <  1:677925/140‑1 (MQ=255)                     atCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                           <  2:594260/140‑1 (MQ=255)                      tCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTcc                            >  2:1183318/1‑138 (MQ=255)                      tCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTcc                            <  1:1183318/138‑1 (MQ=255)                      tCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCAc                          >  2:1182561/1‑140 (MQ=255)                      tCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCAc                          >  1:1014800/1‑140 (MQ=255)                      tCGAATGTGATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCAc                          >  1:1606775/1‑140 (MQ=255)                           tgtgATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCttt                     >  1:774429/1‑140 (MQ=255)                              gATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTgg                  >  1:1182544/1‑140 (MQ=255)                              gATTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTgg                  >  1:1272346/1‑140 (MQ=255)                               aTTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGc                 <  1:678730/140‑1 (MQ=255)                               aTTACAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGc                 >  2:1417558/1‑140 (MQ=255)                                  aCAATCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATg              <  2:878247/140‑1 (MQ=255)                                    aaTCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTc            >  1:412459/1‑140 (MQ=255)                                    aaTCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTc            <  2:806227/140‑1 (MQ=255)                                    aaTCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTc            <  1:158751/140‑1 (MQ=255)                                      tCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                           <  1:796188/123‑1 (MQ=255)                                      tCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                           >  2:796188/1‑123 (MQ=255)                                      tCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg          <  2:1272346/140‑1 (MQ=255)                                      tCGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCAGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg          <  1:480464/140‑1 (MQ=255)                                       cGGTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGc         >  2:85881/1‑140 (MQ=255)                                        ggTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        <  1:226100/140‑1 (MQ=255)                                        ggTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        <  1:1519489/140‑1 (MQ=255)                                        ggTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        <  1:1417558/140‑1 (MQ=255)                                        ggTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        <  1:1373713/140‑1 (MQ=255)                                        ggTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        <  1:780995/140‑1 (MQ=255)                                        ggTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        <  2:32427/140‑1 (MQ=255)                                         gTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCat       >  1:1109310/1‑140 (MQ=255)                                         gTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        <  2:524748/139‑1 (MQ=255)                                         gTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        >  1:965638/1‑139 (MQ=255)                                         gTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        >  1:524748/1‑139 (MQ=255)                                         gTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa        <  2:965638/139‑1 (MQ=255)                                         gTGGCGGCGGCGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGGCGTTGTCCACTCTTTTGGCATGTCAGCat       >  2:889778/1‑140 (MQ=255)                                           ggcggcggcGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCatat     <  2:1606775/140‑1 (MQ=255)                                              ggcggcGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTca  <  2:1109310/140‑1 (MQ=255)                                              ggcggcGCAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTca  <  1:38198/140‑1 (MQ=255)                                                    |                                                                                                                                      GGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCA  >  NC_000913/316347‑316531

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