Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	316,453	G→A	*289* (TAG→TAA)	insF1 →	IS3 transposase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	316,453	0	G	A	100.0%	158.2 / NA	51	*289* (TAG→TAA)	insF1	IS3 transposase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/27);  total (24/27)

GCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTT  >  NC_000913/316349‑316582                                                                                                         |                                                                                                                                  gCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATaaa                                                                                                 <  2:989365/139‑1 (MQ=18) gCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATaaa                                                                                                 <  1:247714/139‑1 (MQ=18)     cGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTcc                                                                                             >  2:208845/1‑139 (MQ=37)        tgtTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCg                                                                                          >  1:981271/1‑139 (MQ=255)          tttAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGCGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCggg                                                                                        >  2:919461/1‑139 (MQ=255)             aaTTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                                                                     <  1:810873/139‑1 (MQ=255)                   atCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTcc                                                                               <  1:208845/139‑1 (MQ=255)                    tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGTGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                              >  2:397329/1‑139 (MQ=255)                    tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                                                                     <  2:659422/132‑1 (MQ=255)                    tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                                                                     >  1:659422/1‑132 (MQ=255)                    tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                              >  2:320350/1‑139 (MQ=255)                    tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                              >  1:274147/1‑139 (MQ=255)                    tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                              >  1:879864/1‑139 (MQ=255)                       aaTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCActc                                                                           >  2:900084/1‑139 (MQ=255)                       aaTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGGGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCActc                                                                           >  2:637263/1‑139 (MQ=255)                         tgtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCtt                                                                         >  2:381880/1‑139 (MQ=255)                         tgtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCtt                                                                         >  2:342139/1‑139 (MQ=255)                              ttACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATaaa                                                                                                 <  2:117464/110‑1 (MQ=18)                              ttACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATaaa                                                                                                 >  1:117464/1‑110 (MQ=18)                                 cAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATg                                                                 <  1:900084/139‑1 (MQ=255)                                 cAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATg                                                                 <  2:484087/139‑1 (MQ=255)                                  aaTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGt                                                                <  1:13360/139‑1 (MQ=255)                                  aaTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGt                                                                <  1:608246/139‑1 (MQ=255)                                   aTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTc                                                               <  2:552537/139‑1 (MQ=255)                                   aTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTc                                                               <  1:301956/139‑1 (MQ=255)                                     cGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg                                                             >  2:572883/1‑139 (MQ=255)                                      ggTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                              >  2:3749/1‑121 (MQ=255)                                      ggTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                              <  1:3749/121‑1 (MQ=255)                                      ggTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGc                                                            >  2:542981/1‑139 (MQ=255)                                       gTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                           <  2:425377/139‑1 (MQ=255)                                       gTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                           <  1:320350/139‑1 (MQ=255)                                        tGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCat                                                          >  1:296839/1‑139 (MQ=255)                                        tGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCat                                                          <  2:859819/139‑1 (MQ=255)                                             gcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTca                                                     <  2:274147/139‑1 (MQ=255)                                             gcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTca                                                     <  1:637263/139‑1 (MQ=255)                                                        ggTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATg                                                                 <  1:932879/116‑1 (MQ=255)                                                        ggTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATg                                                                 >  2:932879/1‑116 (MQ=255)                                                              gcggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCAct                                                                            >  1:889090/1‑99 (MQ=255)                                                              gcggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCaaa                                    <  1:572883/139‑1 (MQ=255)                                                              gcggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCaaa                                    >  2:1014426/1‑139 (MQ=255)                                                              gcggcCTCAGTCCGGAACAATGTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCAct                                                                            <  2:889090/99‑1 (MQ=255)                                                                ggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCaaa                                    >  2:605941/1‑137 (MQ=255)                                                                ggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCaaa                                    <  1:605941/137‑1 (MQ=255)                                                                                  tttGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAAt                >  1:915781/1‑139 (MQ=255)                                                                                      aaaaCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGAc            <  1:758836/139‑1 (MQ=255)                                                                                       aaaCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCat                                                          >  2:573023/1‑92 (MQ=255)                                                                                       aaaCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCat                                                          <  1:573023/92‑1 (MQ=255)                                                                                        aaCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCg          <  2:525070/139‑1 (MQ=255)                                                                                                ccTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTttgtt  <  1:1014426/139‑1 (MQ=255)                                                                                                     cTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                              <  2:405700/58‑1 (MQ=255)                                                                                                     cTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                              >  1:405700/1‑58 (MQ=255)                                                                                                         |                                                                                                                                  GCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTT  >  NC_000913/316349‑316582

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