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%	144.7 / NA	49	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 (21/28);  total (21/28)

CAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACA  >  NC_000913/4507329‑4507594                                                                                                                                       |                                                                                                                                    caacTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTc                                                                                                                                 <  1:344027/139‑1 (MQ=255)   acTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCaa                                                                                                                               <  1:131791/139‑1 (MQ=255)        gTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGt                                                                                                                          <  2:753083/139‑1 (MQ=255)                 gTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAac                                                                                                                 >  2:571278/1‑139 (MQ=255)                                cAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTAt                                                                                                  <  2:35512/139‑1 (MQ=255)                                        tctcCGATAACACCAGAGCAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATAGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTcc                                                                                          <  1:571278/139‑2 (MQ=255)                                                        aaCAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAg                                                                          >  2:917645/1‑139 (MQ=255)                                                           aaaTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAgtgt                                                                       <  1:917645/139‑1 (MQ=255)                                                             aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATGATCTCTTTCGGTGTTTTGAACTTCAGTGTCt                                                                     >  2:669410/1‑139 (MQ=255)                                                             aTCGAAATACGTAAGCCGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAAt                                                                                                                              <  2:255600/82‑1 (MQ=255)                                                             aTCGAAATACGTAAGCCGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAAt                                                                                                                              >  1:255600/1‑82 (MQ=255)                                                               cGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCttt                                                                   >  1:476184/1‑139 (MQ=255)                                                                            cTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCtgttgt                                                      >  1:596736/1‑139 (MQ=255)                                                                             tGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTgttgtt                                                     <  2:681973/139‑1 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCaa                                                                                                                       <  1:781352/71‑1 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCaa                                                                                                                       >  2:781352/1‑71 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGtt                                                                                     <  2:360744/105‑1 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGtt                                                                                     >  1:360744/1‑105 (MQ=255)                                                                               aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTa                                                   >  1:221029/1‑139 (MQ=255)                                                                                ggAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                            <  1:538249/129‑1 (MQ=255)                                                                                ggAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                            >  2:538249/1‑129 (MQ=255)                                                                                     aaGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGag                                             >  2:3814/1‑139 (MQ=255)                                                                                            aaCGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                                      <  2:469566/139‑1 (MQ=255)                                                                                               ggATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAg                                   <  1:3814/139‑1 (MQ=255)                                                                                                 aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCt                                                                                                                            >  2:234423/1‑48 (MQ=255)                                                                                                 aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCt                                                                                                                            <  1:234423/48‑1 (MQ=255)                                                                                                 aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                            >  2:775439/1‑112 (MQ=255)                                                                                                 aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                            <  1:775439/112‑1 (MQ=255)                                                                                                 aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                                 <  2:596736/139‑1 (MQ=255)                                                                                                      aaTGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGt                            <  1:669410/139‑1 (MQ=255)                                                                                                      aaTGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGt                            <  1:648299/139‑1 (MQ=255)                                                                                                          gagaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCAt                        <  1:103262/139‑1 (MQ=255)                                                                                                          gagaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCAt                        <  2:919847/139‑1 (MQ=255)                                                                                                                gAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAg                  <  2:476184/139‑1 (MQ=255)                                                                                                                     atTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTatat             >  2:757157/1‑139 (MQ=255)                                                                                                                        aaaaaGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGTACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATg                       >  1:580585/1‑126 (MQ=255)                                                                                                                        aaaaaGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCAta                       <  2:580585/126‑2 (MQ=25)                                                                                                                         aaaaGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTggg         <  2:698845/139‑1 (MQ=21)                                                                                                                           aaGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCa                                                                           >  2:870990/1‑71 (MQ=255)                                                                                                                           aaGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCa                                                                           <  1:870990/71‑1 (MQ=255)                                                                                                                              gCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                            >  1:666745/1‑83 (MQ=25)                                                                                                                              gCTACTATAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                            <  2:666745/83‑1 (MQ=255)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGtt                                                                                     >  2:93626/1‑56 (MQ=255)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGtt                                                                                     <  1:93626/56‑1 (MQ=255)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTc                          <  1:33159/115‑1 (MQ=14)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTc                          >  2:33159/1‑115 (MQ=14)                                                                                                                                tACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACa  >  1:910237/1‑139 (MQ=14)                                                                                                                                 aCTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGATCTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTatat             <  2:543641/127‑1 (MQ=11)                                                                                                                                 aCTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTatat             >  1:543641/1‑127 (MQ=11)                                                                                                                                       |                                                                                                                                    CAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACA  >  NC_000913/4507329‑4507594

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