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%	190.3 / NA	61	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 (28/33);  total (28/33)

CAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGA  >  NC_000913/4507329‑4507592                                                                                                                                       |                                                                                                                                  caacTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTc                                                                                                                               >  2:160378/1‑139 (MQ=255)    cTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAAt                                                                                                                            <  2:906933/139‑1 (MQ=255)     tGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCaa                                                                                                                             <  1:107893/137‑1 (MQ=255)     tGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCaa                                                                                                                             >  2:107893/1‑137 (MQ=255)     tGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATc                                                                                                                           <  2:535543/139‑1 (MQ=255)          gatgaGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCa                                                                                                                      <  2:26146/139‑1 (MQ=255)          gatgaGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCa                                                                                                                      >  1:26146/1‑139 (MQ=255)             gaGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATg                                                                                                                   <  2:585395/139‑1 (MQ=255)                  tCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAaca                                                                                                              <  2:1308153/139‑1 (MQ=255)                    aGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACAcc                                                                                                            >  1:1610921/1‑139 (MQ=255)                       gCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                         <  1:160378/139‑1 (MQ=255)                        cAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAAt                                                                                                                            <  1:1047455/119‑1 (MQ=255)                        cAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAAt                                                                                                                            >  2:1047455/1‑119 (MQ=255)                          aaaCACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTc                                                                                                      >  1:65375/1‑139 (MQ=255)                              acCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAAtt                                                                                                  <  2:854263/139‑1 (MQ=255)                                          tcCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                         <  2:629064/120‑1 (MQ=255)                                          tcCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                         >  1:629064/1‑120 (MQ=255)                                           ccGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGgtg                                                                                     <  1:804619/139‑1 (MQ=255)                                                    ccAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                         <  1:800004/110‑1 (MQ=255)                                                    ccAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                         >  2:800004/1‑110 (MQ=255)                                                    ccAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCg                                                                                        >  2:1593242/1‑127 (MQ=255)                                                    ccAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCg                                                                                        <  1:1593242/127‑1 (MQ=255)                                                     aaGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAAGGGATTGAAATGGAGAATGAAATATTAAAAAGGCTAATGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCCTTCGGTGTTTTGAACtt                                                                           >  1:1041477/2‑139 (MQ=255)                                                        aaCAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCaa                                                                                                                             <  2:1248983/86‑1 (MQ=255)                                                        aaCAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCaa                                                                                                                             >  1:1248983/1‑86 (MQ=255)                                                             aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCt                                                                   >  1:981869/1‑139 (MQ=255)                                                                  aaTACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTtctc                                                              >  2:688081/1‑139 (MQ=255)                                                                      cGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGGCTTTCTCGCTc                                                          <  2:981869/139‑1 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTGAAAAGGCTACTGTAAAtt                                                                                                                                <  2:941728/60‑1 (MQ=17)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAAtt                                                                                                                                >  1:941728/1‑60 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                         >  2:36791/1‑83 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                         <  1:36791/83‑1 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                          <  1:495073/130‑1 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                          >  2:495073/1‑130 (MQ=255)                                                                                ggAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTAt                                                                                                <  1:776180/91‑1 (MQ=255)                                                                                ggAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTAt                                                                                                >  2:776180/1‑91 (MQ=255)                                                                                          aCAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAgcagca                                      <  1:688081/139‑1 (MQ=255)                                                                                           cAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACAcc                                                                                                            <  1:471128/68‑1 (MQ=255)                                                                                           cAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACAcc                                                                                                            >  2:471128/1‑68 (MQ=255)                                                                                                 aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                                    <  2:1110521/134‑1 (MQ=255)                                                                                                 aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                                    >  1:1110521/1‑134 (MQ=255)                                                                                                 aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                               >  1:222305/1‑139 (MQ=255)                                                                                                  ttGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCttt                                                                                                       >  1:1674635/1‑66 (MQ=255)                                                                                                  ttGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCttt                                                                                                       <  2:1674635/66‑1 (MQ=255)                                                                                                        tGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTc                        <  2:222305/139‑1 (MQ=255)                                                                                                          gagaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCAt                      >  2:1297319/1‑139 (MQ=255)                                                                                                           agaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTtc                                                                >  1:624301/1‑96 (MQ=255)                                                                                                           agaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTtc                                                                <  2:624301/96‑1 (MQ=255)                                                                                                            gaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGt                    <  2:1610921/139‑1 (MQ=255)                                                                                                             aaTGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGtt                   <  1:1297319/139‑1 (MQ=255)                                                                                                               tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCg                                                                                        <  2:1403622/68‑1 (MQ=255)                                                                                                               tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCg                                                                                        >  1:1403622/1‑68 (MQ=255)                                                                                                                    tatTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGtata            <  2:65375/139‑1 (MQ=255)                                                                                                                         aaaaGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTGTGTTGTTTACCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTggg       <  2:1041477/139‑1 (MQ=255)                                                                                                                              gCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                          <  2:69145/83‑1 (MQ=18)                                                                                                                              gCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                          >  1:69145/1‑83 (MQ=18)                                                                                                                              gCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGa  >  2:1297990/1‑139 (MQ=17)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                          >  2:1290070/1‑81 (MQ=17)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                          <  1:1290070/81‑1 (MQ=17)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                                    >  2:1302986/1‑103 (MQ=14)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                                    <  1:1302986/103‑1 (MQ=14)                                                                                                                                       |                                                                                                                                  CAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGA  >  NC_000913/4507329‑4507592

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