Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,401,606 C→T 8.9% P639P (CCG→CCA)  stcD_2 ← putative N‑methylproline demethylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,401,6060CT 8.9% 132.7 / 2.8 45P639P (CCG→CCA) stcD_2putative N‑methylproline demethylase
Reads supporting (aligned to +/- strand):  ref base C (25/16);  new base T (2/2);  total (27/18)
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.87e-01

GACGCAGGGCGGCTAGAGGTCCTTGACCAGCCGCAGGGCATCGTAGATGGCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTGGTTGACCACCACCTGATCCACCAGGCGCGTTTTGGT  >  ACB122/5401474‑5401740
                                                                                                                                    |                                                                                                                                      
gACGCAGGGCGGCTAGAGGTCCTTGACCAGCCGCAGGGCATCGTAGATGGCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGc                                                                                                                                     >  2:21103/1‑136 (MQ=255)
  cGCAGGGCGGCTAGAGGTCCTTGACCAGCCGCAGGGCATCGTAGATGGCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCtt                                                                                                                                   <  3:129005/136‑1 (MQ=255)
        gcggcTAGAGGTCCTTGACCAGCCGCAGGGCATCGTAGATGGCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCTGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGAc                                                                                                                             <  6:143315/136‑1 (MQ=255)
               gagGTCCTTGACCAGCCGCAGGGCATCGTAGATGGCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCa                                                                                                                        <  8:195079/134‑1 (MQ=255)
                                gCAGGGCATCGTAGATGGCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGc                                                                                                     >  7:8340/1‑136 (MQ=255)
                                               tGGCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGc                                                                                      <  4:107902/136‑1 (MQ=255)
                                                 gCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATc                                                                                                                         <  4:111143/99‑1 (MQ=255)
                                                 gCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATc                                                                                                                         >  3:111143/1‑99 (MQ=255)
                                                      aTGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGa                                                                                                                              >  8:197580/1‑89 (MQ=255)
                                                      aTGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGa                                                                                                                              <  7:197580/89‑1 (MQ=255)
                                                      aTGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGt                                                                               >  1:29985/1‑136 (MQ=255)
                                                       tGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGtt                                                                              <  1:47044/136‑1 (MQ=255)
                                                        gggTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCCCCTTCGTTGCTGGCCAGGGGttt                                                                             >  7:189457/1‑136 (MQ=255)
                                                        gggTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGttt                                                                             >  4:150400/1‑136 (MQ=255)
                                                        gggTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGCAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGCCACCTTCGTAGCTCGCCAGGGGttt                                                                             >  4:170330/1‑136 (MQ=255)
                                                          gTATTGCGCGACGCCACTGCATCACCAATGCGGAACATCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGa                                                                           <  6:31578/136‑1 (MQ=255)
                                                                cgcgACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCGTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGGTTGAGGtgg                                                                      >  2:104213/1‑133 (MQ=255)
                                                                       ccACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGTGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAg                                                              >  4:175296/1‑136 (MQ=255)
                                                                       ccACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGTGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAg                                                              >  4:172344/1‑136 (MQ=255)
                                                                       ccACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCg                                                                                                                                       >  7:23010/1‑63 (MQ=255)
                                                                       ccACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCg                                                                                                                                       <  8:23010/63‑1 (MQ=255)
                                                                       ccACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGGTTGAGGTCGAAATACAg                                                              >  3:57542/1‑136 (MQ=255)
                                                                       ccACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGGTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAg                                                              >  5:130973/1‑136 (MQ=255)
                                                                        cACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGt                                                             >  1:142708/1‑136 (MQ=255)
                                                                             cATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCa                                                                                                                        >  1:196154/1‑72 (MQ=255)
                                                                             cATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCa                                                                                                                        <  2:196154/72‑1 (MQ=255)
                                                                                   cAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCa                                                                                                                        >  1:124529/1‑66 (MQ=255)
                                                                                   cAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCa                                                                                                                        <  2:124529/66‑1 (MQ=255)
                                                                                     aTGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGcc                                                >  4:91214/1‑136 (MQ=255)
                                                                                     aTGCGGAACAGCTGGAAGCGCCCTTC‑GGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGCAATACAGTTCGTCCAGCGGCCg                                               >  1:205149/1‑136 (MQ=255)
                                                                                        cGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGGTTGAGGTCGAAATACAGTTCGTCCAGCGGCCggg                                             >  6:110653/1‑136 (MQ=255)
                                                                                              aGCTGCATGCGCCCTTCGGGGTTGCGTACAAGGGTTTGCGGCTTGCCTAAATTCAGGTCCATGTACTCCACGGCACCTGCGTTGCGGGCCAGGGGTTTGAGGTCGAAATCCTGTTCTTCCAGCGGCCGGGTGCCGt                                       <  4:71283/136‑1 (MQ=255)
                                                                                               ggtGGACGCGCCCTGCGGGGTGGCTGACAACGCTGTGCGGCTTCCCGACAATCAGGTCCATGTACTCCACGCCCCCTTGGTTGCTGGCCAGGGGTTTGCGGTCGAAATACTGTTGGTCCAGCGGCCGGGTGCCGTg                                      <  2:87807/134‑1 (MQ=255)
                                                                                                cTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTgg                                     >  4:66675/1‑136 (MQ=255)
                                                                                                  ggAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTCCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTGCGTCCAGCGGCCGGGTGCCGTGGtt                                   <  1:21103/136‑1 (MQ=255)
                                                                                                        gcCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTGGTTGaccac                             >  6:177310/1‑136 (MQ=255)
                                                                                                            ttCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCATCGTTGCTGGCCAGGGGTTTGAGGTAGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTGGTTGaccaccacc                         >  3:71283/1‑136 (MQ=255)
                                                                                                             tCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCCGCGGCCGGGTGCCGTGGTTGACCACCACCt                        >  4:200276/1‑136 (MQ=255)
                                                                                                             tCGGGGTTGCTGACAAAGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTACAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTACGTCCAGCGGCCGGGTGACGTGGTTTACCACCCCCt                        >  2:156976/1‑136 (MQ=255)
                                                                                                                      cTGACAACGGTTTGTGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTGGTTGACCACCACCTGATccacca               <  3:175296/136‑1 (MQ=255)
                                                                                                                      cTGACAACGGTTTGTGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTCCGTCCAGCGGCCGGGTGCCGTGGTTGACCACCACCTGATccacca               <  3:172344/136‑1 (MQ=255)
                                                                                                                       tGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTGGTTGACCACCACCTGATCCACCAg              >  2:28506/1‑136 (MQ=255)
                                                                                                                                tttGCGGCTTGCCGACTATCAGGTCGATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATAGTGTTGGTCCAGCGGCCGGGTGCCGTGGTTGACCACCACCTGATCCACCAGGCGCGtttt     <  3:170330/136‑1 (MQ=255)
                                                                                                                                tttGCGGCTTGCCGACAATCAGGTCCATGTACTCCACGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTGGTTGACCACCACCTGATCCACCAGGCGCGtttt     <  3:150400/136‑1 (MQ=255)
                                                                                                                                   gCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCCCCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTGGTTGACCACCACCTGATCCACCAGGCGCGTTTTGGt  >  8:63964/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GACGCAGGGCGGCTAGAGGTCCTTGACCAGCCGCAGGGCATCGTAGATGGCCGCATGGGTATTGCGCGACGCCACTGCATCACCAATGCGGAACAGCTGGAAGCGCCCTTCGGGGTTGCTGACAACGGTTTGCGGCTTGCCGACAATCAGGTCCATGTACTCCAAGGCACCTTCGTTGCTGGCCAGGGGTTTGAGGTCGAAATACAGTTCGTCCAGCGGCCGGGTGCCGTGGTTGACCACCACCTGATCCACCAGGCGCGTTTTGGT  >  ACB122/5401474‑5401740

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