Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	568,052	T→C	pseudogene (18/213 nt)	renD →	protein RenD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	568,052	0	T	C	96.5%	100.1 / ‑6.7	29	pseudogene (18/213 nt)	renD	protein RenD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (12/16);  minor base G (0/1);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCT  >  NC_000913/567806‑568165                                                                                                                                                                                                                                                       |                                                                                                                  cgagaattaaccgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCgg                                                                                                               <  1:333367‑M2/21‑1 (MQ=255)   agaattaaccgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcctgggcggtagccctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCt                                                                                                             >  2:462025‑M2/229‑251 (MQ=255)   agaattaaccgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcgcaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagtaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgtcggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCgcggat                                                                                                             <  2:556800‑M2/23‑3 (MQ=255)   agaattaaccgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCt                                                                                                             <  2:314239‑M2/23‑1 (MQ=255)     aattaaccgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCAGAAATCAGGCGGTGGCTaa                                                                                                           <  2:274002‑M2/25‑1 (MQ=255)      attaaccgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTaa                                                                                                           >  1:332450‑M2/226‑250 (MQ=255)        taaccgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTaaaa                                                                                                         <  1:99370‑M2/27‑1 (MQ=255)           ccgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGgc                                                                                                     >  1:480489‑M2/221‑251 (MQ=255)            cgtggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGgc                                                                                                     >  2:332819‑M2/220‑250 (MQ=255)               ggtgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAatggcacg                                                                                                  <  2:524069‑M2/34‑4 (MQ=255)                 tgaggcgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGc                                                                                               >  1:606793‑M2/215‑251 (MQ=255)                      cgatccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGcaggtc                                                                                          >  2:29129‑M2/210‑249 (MQ=255)                        atccctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGggcaggtctt                                                                                        <  2:344145‑M2/44‑5 (MQ=255)                           cctgaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattaccgggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCacgggcaggtcttctg                                                                                     <  1:496107‑M2/47‑8 (MQ=255)                              gaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATggcaagggcaggtcttctggtt                                                                                  >  1:216927‑M2/202‑241 (MQ=255)                              gaaccagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGcacgggcaggtcttctgg                                                                                    <  1:2386‑M2/48‑9 (MQ=255)                                  cagtaaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAaaatggcacgggcaggtcttctggttatcg                                                                              >  1:331429‑M2/198‑237 (MQ=255)                                      aaaacaacttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGgctaaaatggcacgggcaggtcttctggttatcgaagg                                                                          <  2:369015‑M2/58‑19 (MQ=255)                                           aacttcctgccatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaagg                                                                     <  2:216927‑M2/63‑24 (MQ=255)                                            acttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcatagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGgccgcggctaaaatggcacgggcagggcttctggttatcgaaggtcaggt                                                                    >  1:467128‑M2/188‑227 (MQ=255)                                            acttcctgtcatgggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGgccgcggctaaaatggcacgggcaggtcttctggttatcgaagttaaggt                                                                    <  1:180954‑M2/64‑25 (MQ=255)                                                         ggcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGgtg                                                       <  1:251491‑M2/77‑1 (MQ=255)                                                          gcggtagacctctaaatcgtgcacaggctctggcgaagatcgcagaaatcaaagcgaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgacggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGTTCTGGCGAACGgtgt                                                      <  2:638803‑M2/78‑1 (MQ=255)                                                                       aaatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACcggtttgc                                         >  2:227532‑M2/161‑248 (MQ=255)                                                                        aatcgtgcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCt                                        >  1:729090‑M2/160‑251 (MQ=255)                                                                              gcacaggctctggcgaagatcgcagaaatcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCag                                   <  2:608027‑M2/97‑1 (MQ=255)                                                                                                         atcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg       <  1:283532‑M2/125‑1 (MQ=255)                                                                                                          tcaaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg       >  1:473267‑M2/126‑250 (MQ=255)                                                                                                           caaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaa     <  1:419521‑M2/127‑1 (MQ=255)                                                                                                           caaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGa      >  1:416149‑M2/125‑250 (MQ=255)                                                                                                            aaagctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACgaac    <  2:199011‑M2/128‑2 (MQ=255)                                                                                                              agctaagttcggactgaaaggagcaagtgtatgacgggcaaagaggcaattattcattacctggggacgcataatagcttctgtgcgccggacgttgccgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCacgaacct  <  1:227532‑M2/130‑4 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                  CTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCT  >  NC_000913/567806‑568165

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