Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	271,967	T→G	pseudogene (204/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,967	0	T	G	82.8%	38.4 / 3.9	29	pseudogene (204/427 nt)	insN	IS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (11/13);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-02

AATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTAT  >  NC_000913/271724‑272189                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                               caacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAg                                                                                                                                                                                                                         >  1:224230‑M1/39‑251 (MQ=255)   acgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGc                                                                                                                                                                                                                        <  1:286888‑M1/214‑1 (MQ=255)               ggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGcc                                                                                                                                                                                                           <  1:147007‑M1/227‑1 (MQ=255)                  gaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGAAACTGGAAAAACCGTCCTGAAAAACCCGACGGCAGACGGGCTGAATTACGCAGTAAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATcgccgc                                                                                                                                                                                                         <  2:573284‑M1/229‑3 (MQ=255)                     tgaaatattaaaaaggctACAGAGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATg                                                                                                                                                                                                     <  2:494678‑M1/233‑1 (MQ=255)                       aaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATCATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTCTTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGc                                                                                                                                                                                                   <  2:336895‑M1/235‑1 (MQ=255)                             taaaaatgctTCCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCGCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGTTCACACTCTGCCATTTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTGCGGCCGGAGCAATAAGCATCGCCACAATGGCAAccc                                                                                                                                                                                              >  2:268440‑M1/11‑250 (MQ=255)                               aaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCaacccgg                                                                                                                                                                                            <  1:540371‑M1/242‑3 (MQ=255)                                  aggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCggagag                                                                                                                                                                                        >  1:336895‑M1/6‑251 (MQ=255)                                   ggctaCCGCGCTCTTGATGTCCGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCATAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTTGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGAACTTGCGCTACATAGCATCAGACACGGTTCGGCCGGAGCAAGACGCATCGCCACAATGGCCACCCggaga                                                                                                                                                                                         >  2:643164/5‑249 (MQ=18)                                   ggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTTTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCggaga                                                                                                                                                                                         >  1:643164‑M1/5‑249 (MQ=255)                                       tccGCGCTCTTGATGTCAGTCTCCCTGAACAGTTCACGATAATCGGGAAACGCAGATCGCATTTACCTTTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAggcta                                                                                                                                                                                    <  2:552735/249‑4 (MQ=17)                                                 tGATGTCAGACTCCCGGAACAGTTCTCGAGAAGCGGGAAACTCAGAGCGCATTATCCTGTGGTCCCACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTTAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCTGGAGAGGCTACCAGATggga                                                                                                                                                                          <  2:568070/250‑2 (MQ=17)                                                  gctGTCAGACTCCCTGAACAGTGCTCGATAAACGGGAAACTCTGAGCGCCTTCTACTGTGGTCACTCTCTGGCAAGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTTGAAAAACCGTCCTGAAAAACCAGACGGCTGACGGGCTGAATTACGCAGTCAGGTACTTGAGCTACATCGCATCAGCACCGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATgggat                                                                                                                                                                         <  2:182460/248‑3 (MQ=9)                                                              acgctctgccgaactaTAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCt                                                                                                                                                              <  2:255338/234‑1 (MQ=14)                                                                     aGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTc                                                                                                                                                     >  1:370654/1‑251 (MQ=14)                                                                      gTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTc                                                                                                                                                     >  1:167215/1‑250 (MQ=14)                                                                           cGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGaa                                                                                                                                                >  1:442682/1‑250 (MQ=14)                                                                             aTAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCt                                                                                                                                                              >  1:255338/1‑234 (MQ=14)                                                                               aaTCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAagag                                                                                                                                            <  2:411672/250‑1 (MQ=14)                                                                                  cGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTgg                                                                                                                                        >  1:649254/1‑251 (MQ=14)                                                                                    ggAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTg                                                                                                                                                                     <  1:271801/220‑1 (MQ=14)                                                                                    ggAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTgggg                                                                                                                                      >  1:389421/1‑251 (MQ=11)                                                                                       aaCTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGGTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTGCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACAAGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTg                                                                                                                                   <  1:4290/251‑1 (MQ=11)                                                                                                          tgtgGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCg                                                                                                                 >  1:323451/1‑250 (MQ=255)                                                                                                                                   ggggTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCgggc                                                                                        >  1:8023/1‑248 (MQ=255)                                                                                                                                                                                      gacggcagacggGCTTTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGATCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGa                                    <  2:8023/251‑1 (MQ=255)                                                                                                                                                                                                                 cAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCTTCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGAtgt         >  1:318322/1‑251 (MQ=255)                                                                                                                                                                                                                         tGAAATGCATGGCACAAGCCACGGCTATGCCGGAGCAAGAAGCATCGCCACAAGGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGTTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACGCCCCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCTGTGATGTGACCtat  <  2:323451/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                               AATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTAT  >  NC_000913/271724‑272189

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