Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,023,365 T→C 72.2% S40P (TCT→CCT)  proP_6 → Proline/betaine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,023,3650TC72.2% 54.9 / 22.0 36S40P (TCT→CCT) proP_6Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base T (3/7);  new base C (12/14);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

AGCTTGAACATGAGCATGCCTCGCAAAGCCGTGGTCGCAGCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCATGACT  >  KT2440/5023239‑5023493
                                                                                                                              |                                                                                                                                
aGCTTGAACATGAGCATGCCTCGCAAAGCCGTGGTCGCAGCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCt                                                                                                                          <  7:15889/135‑1 (MQ=255)
        cATGAGCATGCCTCGCAAAGCCGTGGTCGCAGCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCTCTCAATCTGAAGCCCt                                                                                                                  >  3:267975/1‑135 (MQ=255)
               aTGCCTCGCAAAGCCGTGGTCGCAGCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTc                                                                                                           >  4:185192/1‑135 (MQ=255)
               aTGCCTCGCAAAGCCGTGGTCGCAGCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTc                                                                                                           <  4:219028/135‑1 (MQ=255)
                                  tCGCAGCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTa                                                                                        <  3:185192/135‑1 (MQ=255)
                                  tCGCAGCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTa                                                                                        >  5:175218/1‑135 (MQ=255)
                                       gCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCg                                                                                   <  5:56878/135‑1 (MQ=255)
                                       gCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCg                                                                                   >  8:69594/1‑135 (MQ=255)
                                         tGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGtt                                                                                 >  7:258650/1‑135 (MQ=255)
                                         tGCAGTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGtt                                                                                 >  7:61055/1‑135 (MQ=255)
                                              tAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTa                                                                                        >  8:51174/1‑123 (MQ=255)
                                              tAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTa                                                                                        <  7:51174/123‑1 (MQ=255)
                                                 gCAACGCGACGGAATGGGTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCAtt                                                                          <  3:91340/134‑1 (MQ=255)
                                                        gACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACt                                                                  >  3:138351/1‑135 (MQ=255)
                                                            gAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGc                                                              <  3:290275/135‑1 (MQ=255)
                                                            gAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGc                                                              <  6:97933/135‑1 (MQ=255)
                                                                      aTTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGttt                                                     <  3:180382/134‑1 (MQ=255)
                                                                         aCGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGGATGTTTTggg                                                 >  8:191157/1‑135 (MQ=255)
                                                                          cGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTgggg                                                >  6:148485/1‑135 (MQ=255)
                                                                                    gCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAgg                                       <  4:192359/134‑1 (MQ=255)
                                                                                            gCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAAt                              >  2:111400/1‑135 (MQ=255)
                                                                                              tCACCTTCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTg                            >  2:300827/1‑135 (MQ=255)
                                                                                                   ttCATCTCAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCg                        >  3:102829/1‑134 (MQ=255)
                                                                                                         tcAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCgt                 >  4:266741/1‑135 (MQ=255)
                                                                                                         tcAGCCGCCTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCgt                 <  8:40693/135‑1 (MQ=255)
                                                                                                            gccgccTTGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCgt                 <  6:112603/132‑1 (MQ=255)
                                                                                                            gccgccTTGTTTCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCgt                 <  1:300827/132‑1 (MQ=255)
                                                                                                              cgccTTGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTACGTGACCGAATTGGACGCAAGCGTGTGCt            >  4:30527/1‑135 (MQ=255)
                                                                                                                   tGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCa       <  4:232524/135‑1 (MQ=255)
                                                                                                                     tttCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCa       <  5:282716/133‑1 (MQ=255)
                                                                                                                     tttCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCa       <  1:80749/133‑1 (MQ=255)
                                                                                                                     tttCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCa       <  6:253446/133‑1 (MQ=255)
                                                                                                                     tttCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCa       <  6:175218/133‑1 (MQ=255)
                                                                                                                     tttCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCa       <  1:106748/133‑1 (MQ=255)
                                                                                                                     tttCCTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCAt      <  3:268319/134‑1 (MQ=255)
                                                                                                                         cTGGCCCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCATGACt  <  6:26697/134‑1 (MQ=255)
                                                                                                                              |                                                                                                                                
AGCTTGAACATGAGCATGCCTCGCAAAGCCGTGGTCGCAGCTGCATTAGGCAACGCGACGGAATGGTTCGATTACGGAATCTATGCTTACGGGCTCACCTTCATCTCAGCCGCCTTGTTTCCTGGCTCTCAATCTGAAGCCCTTTACTTCGGCCTGGCAACATTCGCTATCTCGTTCCTCATTCGACCACTTGGCGGTATGTTTTGGGGGCCGCTAGGTGACCGAATTGGCCGCAAGCGTGTGCTGGCCATGACT  >  KT2440/5023239‑5023493

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: ‑