Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	271,876	A→G	pseudogene (113/427 nt)	insN →	IS911A regulator fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	271,876	0	A	G	100.0%	55.8 / NA	18	pseudogene (113/427 nt)	insN	IS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/8);  total (10/8)

CAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTC  >  NC_000913/271635‑272042                                                                                                                                                                                                                                                  |                                                                                                                                                                       aactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaaacgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGaa                                                                                                                                                               >  2:208541‑M1/128‑251 (MQ=255)          atgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTcc                                                                                                                                                      >  1:260631‑M1/119‑251 (MQ=255)                gtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgacatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGaaa                                                                                                                                                 <  1:113516‑M1/138‑1 (MQ=255)                gtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGaaa                                                                                                                                                 <  1:41181‑M1/138‑1 (MQ=255)                    gggcaaaacaccaaaagcctctccgataaaaccagaacaaatcgaaatacgtaagctgaggaaaaagatacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATACTCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGAAGCTACAGATACTGGAAAAACCGtcctgaagaacc                                                                                                                                             <  2:109730‑M1/142‑6 (MQ=255)                                       tctccgataacaccagaacaaatcgaaatacgtaagctgaggcaaaagctacagcgcattgaagtagagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCACTGAACAGTTCTCGATAATCGGGCAACTCAGGGCGCATTATCCTGTGGTAACACTCTGCCGTGTGATCGGAGTGCATCGCAGCAGCTACAGATACTGGAAAAGCCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTAt                                                                                                                         <  1:64543‑M1/162‑1 (MQ=255)                                            gataacaccagaacaaatcgaaatacgtaagctgaggaaaacgctaccacgcattgaaatggagaatgacatattaacaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGc                                                                                                                    <  1:239261‑M1/167‑1 (MQ=255)                                                     agaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTa                                                                                                           >  1:87414‑M1/76‑251 (MQ=255)                                                     agaacaaatcgaaatacggaagcagaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGTAAAAACCGCCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTa                                                                                                           >  2:122177‑M1/76‑251 (MQ=255)                                                                 aatacgtaagcggaggaaaaagcgacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCTCTCTTGATGTCAGACTCCCTGAACAGTTCTCTATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTTGAAAAACCGTCCTGAAAAGCCATACGGCAGACGGGCCGTAGTACGCAGTCTGttacttgagctacat                                                                                               >  2:171125‑M1/64‑244 (MQ=255)                                                                                   aaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTTGGCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATGCTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGccacggttcg                                                                             >  1:89523‑M1/46‑247 (MQ=255)                                                                                   aaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGccacggttcg                                                                             >  2:2174‑M1/46‑247 (MQ=255)                                                                                   aaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGccacggttcg                                                                             <  1:2174‑M1/206‑5 (MQ=255)                                                                                            acgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGc                                                                     <  1:59693‑M1/214‑1 (MQ=255)                                                                                                 ttgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCaaga                                                                 <  1:41345‑M1/218‑1 (MQ=255)                                                                                                                       aaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAacccg                                          >  1:109730‑M1/10‑249 (MQ=255)                                                                                                                                 ccGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCCCACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTACTGAAAACCCAGCCGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCAACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTAcc                               >  1:49575/1‑251 (MQ=18)                                                                                                                                                               gTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACAAGACGGCAGACGGGCTGTATTACGCAGTCAGGGACTTGAGCTACATGACATCAGCAACGGATCGGCCGAAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGCTGGGACGCTGGCTTGCTGGCAgggtc  >  1:225738/1‑246 (MQ=17)                                                                                                                                                                                                                                                  |                                                                                                                                                                       CAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTC  >  NC_000913/271635‑272042

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