Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,886,771 T→C 8.7% F164L (TTC→CTC)  rarD_2 → Protein RarD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,886,7710TC 8.7% 137.8 / 4.6 46F164L (TTC→CTC) rarD_2Protein RarD
Reads supporting (aligned to +/- strand):  ref base T (17/25);  new base C (2/2);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

AGCGGCTGTCACGCCTGCAGACGGTAGCCGTGGCCTGCGCCGCACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGCCCACTGTCCAGTGCTGACCTGCAAGCACACCCC  >  ACB122/3886646‑3886902
                                                                                                                             |                                                                                                                                   
aGCGGCTGTCACGCCTGCAGACGGTAGCCGTGGCCTGCGCCGCACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTgcg                                                                                                                           >  8:6603/1‑136 (MQ=255)
 gCGGCTGTCACGCCTGCAGACGGTAGCCGTGGCCTGCGCCGCACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTgcgc                                                                                                                          <  3:180717/136‑1 (MQ=255)
   ggCTGTCACGCCTGCAGACGGTAGCCGTGGCCTGCGCCGCACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCg                                                                                                                                  <  1:140163/126‑1 (MQ=255)
   ggCTGTCACGCCTGCAGACGGTAGCCGTGGCCTGCGCCGCACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCg                                                                                                                                  >  2:140163/1‑126 (MQ=255)
                           ccGTGGCCTGCGCCGCACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTggg                                                                                                >  3:17275/1‑136 (MQ=255)
                              tGGCCTGCGCCGCACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGggcgg                                                                                             >  1:125871/1‑136 (MQ=255)
                               ggCCTGCGCCGCACTGGGCGGGGGCCACGAGATGTACCAAAATGGCAGTTTTGCGTGGGAAACCCTGGTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGggcggc                                                                                            <  5:161328/136‑1 (MQ=255)
                                gCCTGCGCCGCCGTGGGCGTGGTCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTACCGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGcc                                                                                           <  6:150050/136‑1 (MQ=255)
                                        cgcACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGgtgt                                                                                   >  5:111224/1‑136 (MQ=255)
                                             tGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACAt                                                                              >  3:86801/1‑136 (MQ=255)
                                             tGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACCTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACAt                                                                              >  2:178753/1‑136 (MQ=255)
                                             tGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACCTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACAt                                                                              >  2:181372/1‑136 (MQ=255)
                                                  gTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGc                                                                                                              >  4:43717/1‑99 (MQ=255)
                                                  gTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGc                                                                                                              <  3:43717/99‑1 (MQ=255)
                                                  gTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGcc                                                                         >  5:44595/1‑136 (MQ=255)
                                                            aGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGTCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGcc                                                               >  8:115871/1‑136 (MQ=255)
                                                            aGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCctg                                                                       >  4:96368/1‑128 (MQ=255)
                                                            aGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCctg                                                                       <  3:96368/128‑1 (MQ=255)
                                                            aGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGcc                                                               >  2:22554/1‑136 (MQ=255)
                                                            aGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGcc                                                               >  8:118304/1‑136 (MQ=255)
                                                              cTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGggc                                                                                               >  3:15081/1‑102 (MQ=255)
                                                              cTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGggc                                                                                               <  4:15081/102‑1 (MQ=255)
                                                              cTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGt                                                             <  7:91207/136‑1 (MQ=255)
                                                                          ggCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACt                                                 >  2:74381/1‑136 (MQ=255)
                                                                           gCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCctg                                                                       <  7:10771/113‑1 (MQ=255)
                                                                           gCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCctg                                                                       >  8:10771/1‑113 (MQ=255)
                                                                           gCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACt                                                 <  4:45068/135‑1 (MQ=255)
                                                                           gCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACt                                                 <  2:67527/135‑1 (MQ=255)
                                                                           gCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACt                                                 <  7:118304/135‑1 (MQ=255)
                                                                              gTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGt                                             >  8:100671/1‑136 (MQ=255)
                                                                              gTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGt                                             >  1:72093/1‑136 (MQ=255)
                                                                                    ccTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATcc                                        <  7:78168/135‑1 (MQ=255)
                                                                                        ggACACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGGGCGGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGc                                   <  6:31539/136‑1 (MQ=255)
                                                                                        ggAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGTCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCTCTGTACTTCGTGATCCAGAGc                                   <  1:159022/136‑1 (MQ=255)
                                                                                        ggAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGc                                   <  7:6603/136‑1 (MQ=255)
                                                                                        ggAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGc                                   <  8:87875/136‑1 (MQ=255)
                                                                                        ggAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGc                                   <  6:67013/136‑1 (MQ=255)
                                                                                        ggAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGc                                   <  5:3337/136‑1 (MQ=255)
                                                                                        ggAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGc                                   <  2:125871/136‑1 (MQ=255)
                                                                                         gAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGcc                                  <  1:22554/136‑1 (MQ=255)
                                                                                             cccTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGCCCACt                              <  6:111224/136‑1 (MQ=255)
                                                                                                               gCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGCCCACTGTCCAGTGCTGACCTGCa            <  3:20188/136‑1 (MQ=255)
                                                                                                                cTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGCCCACTGTCCAGTGCTGACCTGCaa           <  1:74381/136‑1 (MQ=255)
                                                                                                                    ccGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGCCCACTGTCCAGTGCTGACCTGCAAGcac       <  4:17275/136‑1 (MQ=255)
                                                                                                                         cTACCTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGCCCACTGTCCAGTGCTGACCTGCAAGCACAcccc  <  1:181372/136‑1 (MQ=255)
                                                                                                                         cTACCTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGCCCACTGTCCAGTGCTGACCTGCAAGCACAcccc  <  1:178753/136‑1 (MQ=255)
                                                                                                                             |                                                                                                                                   
AGCGGCTGTCACGCCTGCAGACGGTAGCCGTGGCCTGCGCCGCACTGGGCGTGGGCCACGAGCTGTACCAAAATGGCAGTTTTGCCTGGGAAACCCTGCTGGTAACGATCGGCTACCCGGTCTACTTCGTGCTGCGCCGACGCTGCCGCACCGACAACCTGGGCGGCCTGTGGTGTGACATGTGCCTGCTGCTGCCGTGGGCGCTGTACTTCGTGATCCAGAGCCCACTGTCCAGTGCTGACCTGCAAGCACACCCC  >  ACB122/3886646‑3886902

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