Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	271,946	G→A	pseudogene (183/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,946	0	G	A	100.0%	65.1 / NA	30	pseudogene (183/427 nt)	insN	IS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/10);  total (20/10)

CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑ACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGC  >  NC_000913/271706‑272073                                                                                                                                                                                                                                                  |                                                                                                                                aagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCa                                                                                                                       >  2:162770‑M1/57‑251 (MQ=255)           aaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCcacggttc                                                                                                             >  1:95799‑M1/47‑248 (MQ=255)             aagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAgccacggttcgg                                                                                                           >  1:118805‑M1/45‑246 (MQ=255)                      cgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCa                                                                                                   <  1:458651‑M1/215‑1 (MQ=255)                       gcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCa                                                                                                   <  1:231200‑M1/215‑1 (MQ=255)                               atggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCg                                                                                         >  2:223867‑M1/27‑251 (MQ=255)                                  gagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGcc                                                                                       <  1:373173‑M1/227‑1 (MQ=255)                                    gaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCcaca                                                                                    >  1:255084‑M1/22‑251 (MQ=255)                                      atgaaatattaaaaaggctACCGCGCGCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGTCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAAt                                                                                  >  2:588519‑M1/20‑251 (MQ=255)                                        gaaatattaaaaaggctACCGCGCTCTTGATGTCAGCCTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGGGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCGGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAAt                                                                                  <  1:455198‑M1/232‑1 (MQ=255)                                               taaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAacccg                                                                         <  1:575213‑M1/241‑2 (MQ=255)                                               taaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCACGATAATCGGGAAACTCAGAGCGCATTATCCGGTGGTCACACTCTGCCATGTGTTCGGGGTTCAGCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTGGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAAccc                                                                          >  2:422112‑M1/11‑250 (MQ=255)                                                        tACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑GCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCAcaatggcaacccggcgaggct                                                                 <  1:334191‑M1/249‑10 (MQ=255)                                                          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)                                                                                                                                                                                                                                                  |                                                                                                                                CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGA‑ACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGC  >  NC_000913/271706‑272073

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