Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,313,950 T→G 42.8% R612R (AGG→CGG)  CJ019_01235 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,313,9500TG42.8% 8.3 / 5.9 29R612R (AGG→CGG) CJ019_01235hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (12/3);  new base G (2/11);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCACCTGC  >  KT2440/1313819‑1314085
                                                                                                                                   |                                                                                                                                       
ggTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTgcgc                                                                                                                                     >  3:219418/1‑136 (MQ=255)
       ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGGCCGGCGCCGGGCg                                                                                                                               >  5:252700/1‑135 (MQ=255)
       ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCGGCGCAGGGCgg                                                                                                                              >  3:594/1‑135 (MQ=255)
                      tCCTGTGGATGCCGCCAGCGCACCAACGCCTCGACCGCGACGATACGCAGGTCGTCGAGGCGCTAGCGTGGCTTGTTATGCAACTAGAACTCCTCATTGAGCAGGGACCAGCGCAGGTCGTTCGACCGCTggcgg                                                                                                                >  1:224214/1‑135 (MQ=255)
                                   gccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATct                                                                                                   >  8:295173/1‑135 (MQ=255)
                                               aGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCAttt                                                                                       >  8:127831/1‑135 (MQ=255)
                                                        ccGCGCCGATTCGCAGGTCGTCCCGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCGGGCGGCGGTCCTGGATCTGCGGGTTCATGTCGGCAAc                                                                                >  8:96638/1‑133 (MQ=255)
                                                        ccGCGCCGATGCGCCGGGCGGCCCGGGGGTAGGGTGGCTGGTAAGCCAACTCCACCCCCACATTGCCCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGa                                                                              <  8:302216/135‑1 (MQ=255)
                                                        ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGTTCTGCTGGTTCATTTCGGCCAc                                                                                >  7:109411/1‑133 (MQ=255)
                                                        ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGCTCTGCTGGTTCATTTCGGCCAc                                                                                >  7:133654/1‑133 (MQ=255)
                                                              cgaTTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGGCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGGAACGAAGTAgc                                                                        >  1:226194/1‑135 (MQ=255)
                                                                atgcGCTGGGCGGCCCGGCGGTAGCGTGGCGGGAAATCAAACTCCAACTCCTCATGGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAgcgc                                                                      <  8:295629/132‑1 (MQ=255)
                                                                 tgcgtgggTCGTGCGGGGTGTAGGGTGGCTGGTAATCCAACTCCACCTCCTCTTGGCGGGGGGCGCGCCGAAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAgcgc                                                                      <  8:239040/128‑1 (MQ=255)
                                                                             ccggggggTAGCGTGGCTGGTAATCCAACTCGAACTCCTCATTCCGCCGTGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTc                                                         <  8:133654/129‑1 (MQ=255)
                                                                                  cGGTTGCGTGGCGGGTAATGCAACTCCACCGCCTCATTGCGCAGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                                   <  4:86254/136‑1 (MQ=255)
                                                                                  cGGTTGCGGGGCTGGAAATGCACCTCCAACGCCCCTTGGCGCCGGGCCGGGCGCAGGTCGTTCTCCAGCTGGCGGCGGGACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                                   <  3:232193/136‑1 (MQ=255)
                                                                                   ggTTGCGTGGCTTGTACTCCAACTCTACCACCTCATTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                                   <  7:245483/135‑1 (MQ=255)
                                                                                   ggTCGCGTGGGTGTTAATCAAACTCCAACTCCTCATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                                   <  5:64545/135‑1 (MQ=255)
                                                                                       ggggggCTGGAAATGCAACGCGAACACCTCATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                              <  3:291855/132‑1 (MQ=255)
                                                                                       gCGGGGCTGGAAATCCACCTCCACCTCCTAATTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCtt                                               <  4:594/135‑1 (MQ=255)
                                                                                        ggTGGCTGGTACCCAAAGTCGAACTCGTCTTGGCTCAGGCCCCGGCGCAGGTCGTTCTCGAGCTGGCGGCGGTACTGGATCTGCTGTTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                              <  3:211576/134‑1 (MQ=255)
                                                                                            ggtttaaaTGCCATTCCAACTCCGCTTTGCGCTGGGCCGGCCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTg                                         <  3:330228/130‑1 (MQ=255)
                                                                                                gTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGGGCGGGTGGTTCTCCAGCTGGCGGCGGTTCTGGTTCTTCTTGGTCAGTTAGGGCAAGAAGTTGCGCCAAGTGATCTTTCCAGCCGCCTTGGCCAGGtt                                      >  2:185026/1‑134 (MQ=255)
                                                                                                               ctcctcATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGgcg                      <  3:309384/136‑1 (MQ=255)
                                                                                                                          gccgggCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCg                 <  4:72405/127‑1 (MQ=255)
                                                                                                                          gcAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCg                 >  3:72405/1‑130 (MQ=255)
                                                                                                                          gcAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCa            >  7:349278/1‑135 (MQ=255)
                                                                                                                               gCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCa       >  1:229372/1‑135 (MQ=255)
                                                                                                                               gCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCAc      >  4:301298/1‑136 (MQ=255)
                                                                                                                                   tGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCACCTGc  <  3:301298/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
GGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCACCTGC  >  KT2440/1313819‑1314085

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