Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 1,909,455 G→A Q113* (CAG→TAG)  proQ ← RNA chaperone, putative ProP translation regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,909,4550GA96.5% 88.4 / NA 29Q113* (CAG→TAG) proQRNA chaperone, putative ProP translation regulator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (6/22);  total (7/22)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTGCGCCGGGTTTAACACCGTAAAGATAACGCCAGCTCGAAGT  >  CP009273/1909309‑1909599
                                                                                                                                                  |                                                                                                                                                
cGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCtatg                                                                                                                                                <  1:170770/149‑1 (MQ=255)
 gggttttGCCCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGGGCGTCCTCTTTCTCACCAGCCGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGc                                                                                                                                               <  2:393178/146‑1 (MQ=255)
        gcgcACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAaca                                                                                                                                        >  1:299346/1‑149 (MQ=255)
             cGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACgcgc                                                                                                                                   >  1:641398/1‑149 (MQ=255)
                ggTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCttt                                                                                                                                <  2:668800/149‑1 (MQ=255)
                   ttccGTTCAGCGCCTTCGTTGCGGCGTGGCGTAGGCGGACGTGGTCTGCCTTCGCGGCGCGGGGCGGCTTCTTTCTCACCAGCAGTTGCCGCAGCTTCGCGTTTTTTCGCTTGCTCTTCAGCACGCTATGCCTGAACACGCGCTTTCGc                                                                                                                             <  2:18289/146‑1 (MQ=255)
                     aCGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGGTGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGctt                                                                                                                           <  1:627105/149‑1 (MQ=255)
                      cGTTCAGCGCCTTCTTGGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGcttc                                                                                                                          <  1:582002/149‑1 (MQ=255)
                       gTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGcttct                                                                                                                         <  2:682086/149‑1 (MQ=255)
                                      tGCGTCGTGGCGTAGGCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCCTTCTGCCCAACCGTATCGGCAGCTTCGCGTTTTTTCTCTTGCTGTTCAGCACCCTAGGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGa                                                                                                          <  1:161729/149‑1 (MQ=255)
                                              ggcgtAGTCGGACGTGGCTTGCGTTCCCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCTAGCATGctc                                                                                                  <  2:488007/149‑1 (MQ=255)
                                                     tCGGACGTGGCTTGCGGTCGCGGCGAGGTGCGTCTTCTTTCTCACCACCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGt                                                                                           <  1:427663/149‑1 (MQ=255)
                                                          cGTGGCTTGCGGTCCCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTc                                                                                      <  1:156085/149‑1 (MQ=255)
                                                                    gTTCGCGGCGCTGTGCGTCTTCTGTCTCACCAGCAGTTAAGGAAGGTTTTCGTTTTGGCGCTTGCTGTTCAGTACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCa                                                                            <  1:397460/149‑1 (MQ=255)
                                                                        gcggcgcggTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGATTCGCGTTTTTTCGCTTGCTGTTAAGCAAGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTATACATGTTGCTCCTCTGGCTccacgc                                                                        >  2:626122/1‑144 (MQ=255)
                                                                                   cGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCg                                                             <  2:645445/149‑1 (MQ=255)
                                                                                      cttcttTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGGCa                                                          >  1:222217/1‑149 (MQ=255)
                                                                                         ctttctCACTCGCAGTTGCGGCAGCTTCGCGTTTTTTCGCGTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGa                                                       <  2:153973/149‑1 (MQ=255)
                                                                                          tttcgCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGccaccgtatt                                                      <  2:375786/149‑9 (MQ=255)
                                                                                              tcACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCCGGTCCAGCTCACCGCATGGGTTGCCGTCAAGGTCgcc                                                  >  2:28143/1‑147 (MQ=255)
                                                                                                aCCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGacac                                                >  2:319634/1‑149 (MQ=255)
                                                                                                   agcagTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCCGCTTGTGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACAcgcg                                             <  1:566904/149‑1 (MQ=255)
                                                                                                      agTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTg                                          <  2:199723/149‑1 (MQ=255)
                                                                                                          gCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCTTCAAGATCGACACGCGTTGCGcc                                      <  1:454753/149‑1 (MQ=255)
                                                                                                                    cgcgTTTTTTCGCTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTGCGCCGGGTTTAaca                            <  2:295756/149‑1 (MQ=255)
                                                                                                                           tttCGTTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGTTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTGCGCCGGGTTTAACACCGTaaa                     <  2:24310/149‑1 (MQ=255)
                                                                                                                                cTTGCTGTTCAGCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTGCGCCGGGTTTAACACCGTAAAGATaa                <  1:108015/149‑1 (MQ=255)
                                                                                                                                           gCACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTGCGCCGGGTTTAACACCGTAAAGATAACGCAAGCTCGa     <  1:690225/149‑1 (MQ=255)
                                                                                                                                            cACGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTGCGCCGGGTTTAACACCGTAAAGATAACGCCAGCTCGaa    <  1:255909/149‑1 (MQ=255)
                                                                                                                                              cGCTATGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTGCGCCGGGTTTAACACCGTAAAGATAACGCCAGCTCGAAgt  >  1:255941/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTCGCGTTTTTTCGCTTGCTGTTCAGCACGCTGTGCCTGAACACGCGCTTTCGCTTCTTCAAGCTGCTTGCGAGCATGCTCTACATGTTGCTCGTCCAGCTCACCGCATGGGTTGCCGTCAAGATCGACACGCGTTGCGCCGGGTTTAACACCGTAAAGATAACGCCAGCTCGAAGT  >  CP009273/1909309‑1909599

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