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%	143.8 / NA	46	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/25);  total (21/25)

GTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGG  >  NC_000913/4507336‑4507587                                                                                                                                |                                                                                                                             gTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGt                                                                                                                   <  1:884137/139‑1 (MQ=255)          gTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAac                                                                                                          >  2:394099/1‑139 (MQ=255)                                      gATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGtt                                                                              <  1:1049575/139‑1 (MQ=255)                                       aTAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGttt                                                                             >  1:691143/1‑139 (MQ=255)                                            acCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                    >  1:519923/1‑111 (MQ=255)                                            acCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCt                                                                                                    <  2:519923/111‑1 (MQ=255)                                               aGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTc                                                                     <  2:1413215/139‑1 (MQ=255)                                                      aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCt                                                              >  1:1151629/1‑139 (MQ=255)                                                               cGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTc                                                     <  2:1151629/139‑1 (MQ=255)                                                                gTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCt                                                    >  2:307142/1‑139 (MQ=255)                                                                  aaGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCtgt                                                  >  2:321900/1‑139 (MQ=255)                                                                     cTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGtt                                                                              <  2:1438135/108‑1 (MQ=255)                                                                     cTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGtt                                                                              >  1:1438135/1‑108 (MQ=255)                                                                        aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGt                                                                                                                   <  1:435156/68‑1 (MQ=255)                                                                        aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGt                                                                                                                   >  2:435156/1‑68 (MQ=255)                                                                        aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAgt                                                                  >  1:723187/1‑117 (MQ=255)                                                                        aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAgt                                                                  <  2:723187/117‑1 (MQ=255)                                                                        aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTgttgttca                                            >  1:1373883/1‑137 (MQ=255)                                                                        aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTa                                            >  1:677385/1‑139 (MQ=255)                                                                        aGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTa                                            >  2:174220/1‑139 (MQ=255)                                                                         ggAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAg                                           >  1:914339/1‑139 (MQ=255)                                                                          gAAAAAGCTACAACGGATTGAAATGGAGACTGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGc                                          <  1:1020853/139‑1 (MQ=255)                                                                             aaaGCTACAACCGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAGTTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGa                                       >  2:873484/1‑139 (MQ=255)                                                                                gCTACAACGGATTGAAATGGAGTATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagca                                    <  2:691143/139‑1 (MQ=255)                                                                                 cTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagcag                                   <  1:307142/139‑1 (MQ=255)                                                                                 cTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagcag                                   <  1:174220/139‑1 (MQ=255)                                                                                      aCGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAcc                              <  1:321900/139‑1 (MQ=255)                                                                                          aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          <  2:1082575/139‑1 (MQ=255)                                                                                          aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          <  2:914339/139‑1 (MQ=255)                                                                                          aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          >  2:311316/1‑139 (MQ=255)                                                                                          aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGAt                          >  1:116053/1‑139 (MQ=255)                                                                                                 tGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAg                                           >  2:332389/1‑115 (MQ=255)                                                                                                 tGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAg                                           <  1:332389/115‑1 (MQ=255)                                                                                                    agaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCGATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTGTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagc                                     >  1:927323/1‑118 (MQ=255)                                                                                                    agaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATg                <  1:873484/139‑1 (MQ=255)                                                                                                    agaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATg                <  2:565852/139‑1 (MQ=255)                                                                                                    agaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTCTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagc                                     <  2:927323/118‑1 (MQ=255)                                                                                                        tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCttt                                                                                     <  1:937614/66‑1 (MQ=255)                                                                                                        tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCttt                                                                                     >  2:937614/1‑66 (MQ=255)                                                                                                        tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGaa                                                                         >  2:1427707/1‑78 (MQ=255)                                                                                                        tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGaa                                                                         <  1:1427707/78‑1 (MQ=255)                                                                                                        tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGa            <  1:394099/139‑1 (MQ=255)                                                                                                             tatTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagca                                    >  1:693419/1‑110 (MQ=255)                                                                                                             tatTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagca                                    <  2:693419/110‑1 (MQ=255)                                                                                                             tatTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTCAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGtata       <  2:1373883/139‑1 (MQ=255)                                                                                                                  aaaaGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTggg  <  2:116053/139‑1 (MQ=21)                                                                                                                                |                                                                                                                             GTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGG  >  NC_000913/4507336‑4507587

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