Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,507,463	G→A	intergenic (+12/+3)	insO → / ← insI1	pseudogene, IS911 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related/IS30 transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,507,463	0	G	A	100.0%	180.3 / NA	59	intergenic (+12/+3)	insO/insI1	pseudogene, IS911 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related/IS30 transposase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/33);  total (26/33)

CTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGG  >  NC_000913/4507332‑4507587                                                                                                                                    |                                                                                                                             cTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAAt                                                                                                                       <  2:1526676/139‑1 (MQ=255)     gTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGt                                                                                                                   <  1:717361/139‑1 (MQ=255)      tgatgaGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACGACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTc                                                                                                                  >  1:815638/1‑139 (MQ=255)                   ggCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACAcccc                                                                                                     <  2:411447/139‑1 (MQ=255)                    gCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                    <  1:1074903/139‑1 (MQ=255)                     cAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCtt                                                                                                   >  1:1130627/1‑139 (MQ=255)                        aaCACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCa                                                                                                >  1:242279/1‑139 (MQ=255)                                   ccTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCttt                                                                                     <  2:1242837/139‑1 (MQ=255)                                      ctcCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCgg                                                                                  >  1:491384/1‑139 (MQ=255)                                       tcCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGgt                                                                                 >  1:1385509/1‑139 (MQ=255)                                                acCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAAc                                                                        >  1:1398108/1‑139 (MQ=255)                                                   aGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTc                                                                     <  2:1144757/139‑1 (MQ=255)                                                      aCAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAgt                                                                  <  2:815638/139‑1 (MQ=255)                                                          aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAac                                                                                                          <  1:1031275/95‑1 (MQ=255)                                                          aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAac                                                                                                          >  2:1031275/1‑95 (MQ=255)                                                          aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                    >  2:1259749/1‑101 (MQ=255)                                                          aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                    <  1:1259749/101‑1 (MQ=255)                                                          aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCt                                                              >  2:42807/1‑139 (MQ=255)                                                          aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCt                                                              >  2:863155/1‑139 (MQ=255)                                                               aaTACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTtctc                                                         >  2:1202265/1‑139 (MQ=255)                                                                 tACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCgg                                                       >  2:917995/1‑139 (MQ=255)                                                                           gAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCg                                                                                   <  2:426870/101‑1 (MQ=255)                                                                           gAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCg                                                                                   >  1:426870/1‑101 (MQ=255)                                                                            aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTa                                            >  1:1448930/1‑139 (MQ=255)                                                                            aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTa                                            >  2:223315/1‑139 (MQ=255)                                                                             ggAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAg                                           >  2:1000791/1‑139 (MQ=255)                                                                              gAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGc                                          >  1:336001/1‑139 (MQ=255)                                                                                   aGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACAcccc                                                                                                     >  1:444251/1‑75 (MQ=255)                                                                                   aGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACAcccc                                                                                                     <  2:444251/75‑1 (MQ=255)                                                                                    gCTACAACGGACTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagca                                    <  1:1000791/139‑1 (MQ=255)                                                                                         aaCGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                               <  1:1232495/139‑1 (MQ=255)                                                                                          aCGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAcc                              <  1:42807/139‑1 (MQ=255)                                                                                              aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTTTGTTGTTTAGCTGAGCAGCAACCAGAt                          <  2:1398108/139‑1 (MQ=255)                                                                                              aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          <  1:223315/139‑1 (MQ=255)                                                                                              aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          <  2:1130627/139‑1 (MQ=255)                                                                                              aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          >  2:886807/1‑139 (MQ=255)                                                                                              aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          <  2:336001/139‑1 (MQ=255)                                                                                              aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          <  2:1448930/139‑1 (MQ=255)                                                                                                   aaTGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGt                     <  2:366498/139‑1 (MQ=255)                                                                                                   aaTGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGt                     <  2:242279/139‑1 (MQ=255)                                                                                                         gaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGt               <  2:1385509/139‑1 (MQ=255)                                                                                                            tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGa            <  2:491384/139‑1 (MQ=255)                                                                                                             gAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAg           <  1:68360/139‑1 (MQ=255)                                                                                                              aaaTATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCAt                 >  1:891489/1‑132 (MQ=255)                                                                                                              aaaTATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCACATCTAgcttat                 <  2:891489/132‑6 (MQ=255)                                                                                                                 tatTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGtata       <  1:863155/139‑1 (MQ=38)                                                                                                                    tAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                               >  2:936261/1‑112 (MQ=255)                                                                                                                    tAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                               <  1:936261/112‑1 (MQ=38)                                                                                                                      aaaaGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTggg  <  1:886807/139‑1 (MQ=37)                                                                                                                           gCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGtata       >  2:433944/1‑129 (MQ=18)                                                                                                                           gCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGtata       >  2:276649/1‑129 (MQ=17)                                                                                                                           gCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGtata       <  1:276649/129‑1 (MQ=18)                                                                                                                           gCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGtata       <  1:433944/129‑1 (MQ=18)                                                                                                                             tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                     >  2:85710/1‑81 (MQ=17)                                                                                                                             tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                     <  1:85710/81‑1 (MQ=17)                                                                                                                             tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAg                                           <  2:724448/91‑1 (MQ=17)                                                                                                                             tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAg                                           >  1:724448/1‑91 (MQ=17)                                                                                                                              aCTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagcag                                   <  1:1172272/98‑1 (MQ=14)                                                                                                                              aCTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagcag                                   >  2:1172272/1‑98 (MQ=14)                                                                                                                                    |                                                                                                                             CTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGG  >  NC_000913/4507332‑4507587

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