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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	272,002	0	A	G	100.0%	22.6 / NA	17	pseudogene (239/427 nt)	insN	IS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/4);  total (13/4)

CCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑ACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGC  >  NC_000913/271763‑272073                                                                                                                                                                                                                                                 |                                                                        ccGCGATCTTGATGTCAGACTCCATGCACAGTTCTCGGTAATAGGAAAACTCAGGGCGCATTATCCTGTGCTCACACTCTGCATTGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAACAACAGTCCTGAAAAACCAGCCGGCAGACGGGCTGTATTACGCAGTAAGGTACTTGAGAC‑ACATGGCATCAGCCACGGTTCGGCCGGAGCCAGACGCATCGCCACAATGGCAACCCGGAGAggcaac                                                               >  2:561369/1‑246 (MQ=12)   gcgcTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCGACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCAGCGCCACAATGGCAACCCGGAGAGGCTACCAg                                                            <  2:90080/251‑1 (MQ=18)          tGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGTCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATgggac                                                     >  2:468080/1‑249 (MQ=18)            aTGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GATACATGGAATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGATACCAGATGGGacgc                                                   >  1:560112/1‑251 (MQ=18)                 aGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCGTGTGTTCGGGGTTCATCGCAGCAGCTCCAGATACTGGAAAAACCGTCCTGAAAAACAAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCAGGGTTCGGCCGGAGCAAGAAGCATAGCCACAATGGCAAACCGGAGAGGCTACCAGATGGGACGCTGGCt                                              >  1:595105/1‑251 (MQ=17)                  gACTCCATGAACAGTTCTCACTAATCGGGAAACTCAGAGCGCATTATCCTGTGGACACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACACTGGCAACCCGGAGAGGCTCCCCGATGGGACGCTGGCtt                                             >  2:78851/1‑251 (MQ=14)                      cccTGAACCGTTCTCGATAATCGGGAGACTCAGAGCGCATTTTCCTGTGGTCACACTCTGCCAGGTGTTCGGGGTTCATCGCAGCAGCTACAGAGACTGGAAAACCCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGTCAACCCGGAGAGGCAACCAGATGGGACGCTGGCTTGCTg                                         <  2:455198/251‑1 (MQ=14)                      cccTGAACAGTTCTCGATAAACGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCAGCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTtgcttgct                                          <  2:414735/250‑1 (MQ=17)                          gAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGc                                       >  2:25465/1‑249 (MQ=17)                             cAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGACAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAggcc                                  >  2:254822/1‑248 (MQ=14)                                 tctcGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCAt                               <  2:161756/250‑1 (MQ=14)                                  atcGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATg                              >  2:350575/2‑250 (MQ=14)                                     gATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGaa                            >  1:470428/1‑249 (MQ=14)                                                  cTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTGGt             >  1:467952/1‑251 (MQ=11)                                                   tCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTGGTc            >  2:193104/1‑251 (MQ=11)                                                             ttATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTGGTCAGCTGTcagc  >  2:72973/1‑251 (MQ=11)                                                              tATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTGGTCAGCTGTcagc  >  1:573826/1‑250 (MQ=11)                                                                                                                                                                                                                                                 |                                                                        CCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑ACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGC  >  NC_000913/271763‑272073

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