Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,259,465 C→T 9.7% V32M (GTG→ATG)  CJ019_03060 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,259,4650CT 9.7% 118.4 / 4.1 41V32M (GTG→ATG) CJ019_03060hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base T (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

ACGGCGCCGGCCGTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATACCGCTCTCCACAAGCCTCGAATGACGCCGCCAGGTG  >  KT2440/3259339‑3259595
                                                                                                                              |                                                                                                                                  
aCGGCGCCGGCCGTATCCCGGGTCGATCGCAATAACCCAGCAACAGCTCGTCGGCTTCCTTGAACGACAACCCGGCCCCCGCCCGGTCCAGCACCGCGAACCACTTGTTTTGTGGCACGGCAAACACTCGGATgc                                                                                                                            >  2:193901/1‑135 (MQ=255)
 cGGCGCCGGCCGTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGCGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACCCTCGGATgcc                                                                                                                           >  8:186891/1‑135 (MQ=255)
  ggCGCCGGCCGTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATgccg                                                                                                                          <  7:119279/135‑1 (MQ=255)
  ggCGCCGGCCGTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATgccg                                                                                                                          <  7:233217/135‑1 (MQ=255)
     gccggccgTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCTCCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGccc                                                                                                                       <  4:149204/135‑1 (MQ=255)
       cggccgTATCCCGGGTCGATCGCAATAACCCAGAAACCGCTCGTCGGCGACCGTGCACGACACCCCGGCCCCCGCCCGGTCGCGCCCCGCGCACCTNTTGGGGTGTGCCACCGCACACCCGCGGCTGCCGccccg                                                                                                                     >  5:48285/1‑134 (MQ=255)
       cggccgTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTCCCTTGAACGACCCCCCGGCACCCGCCCGGTCGCGCACCGCGCCCATCGTGGGTTCTGGCAGCGCAACCTCTCGGAGGCCGCCCCa                                                                                                                     >  5:43138/1‑135 (MQ=255)
            gTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCAttt                                                                                                                  >  8:366937/1‑133 (MQ=255)
            gTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCAttt                                                                                                                  >  3:330226/1‑133 (MQ=255)
                   gggTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTc                                                                                                         <  5:12366/135‑1 (MQ=255)
                       cGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCcgc                                                                                                     <  3:154992/135‑1 (MQ=255)
                            gCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCAt                                                                                                                    >  3:288115/1‑115 (MQ=255)
                            gCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCAt                                                                                                                    <  4:288115/115‑1 (MQ=255)
                            gCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCAGTTGTAGTCTTCCCGCGCAcc                                                                                                >  1:162639/1‑135 (MQ=255)
                                   cccAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAgg                                                                                         <  1:109182/135‑1 (MQ=255)
                                     cAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCCCCCGGCAGGCt                                                                                       >  3:309164/1‑135 (MQ=255)
                                                cgtcgGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCCCCGCAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATg                                                                            >  6:335457/1‑135 (MQ=255)
                                                 gtcgGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATgc                                                                           <  1:96312/135‑1 (MQ=255)
                                                      cTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCa                                                                      <  4:330226/135‑1 (MQ=255)
                                                            tcaaCGACAACCCGGCACCCGCCAGGTCGAGCACCGCTAACATCTTTTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTggg                                                                <  1:193901/133‑1 (MQ=255)
                                                               aCGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGcc                                                              <  6:228350/134‑1 (MQ=255)
                                                                      cccGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCAt                                                      <  2:262299/135‑1 (MQ=255)
                                                                      cccGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCAt                                                      >  5:54768/1‑135 (MQ=255)
                                                                              ccGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGtt                                              <  8:120269/135‑1 (MQ=255)
                                                                                   aGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGAc                                         <  5:132162/135‑1 (MQ=255)
                                                                                        gAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCCCCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACCTAcc                                    >  2:152518/1‑135 (MQ=255)
                                                                                        gAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATAcc                                    <  6:311113/135‑1 (MQ=255)
                                                                                                cgAACATCTTGTTTTCTGCCACCGAATACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATACCGCTCTCca                            <  8:281327/135‑1 (MQ=255)
                                                                                                 gAACATCTTGTTTTCTGCCACCGAAAACACTCGGCTGCCGCCCCATTTGTAGTCTTCCCGCGCACACGGCAGGCTCCGGCAAAATGCCGCCACCTGGGCATCGCTCATCTTCCGTTCAGACATACCGCTATCca                            >  1:341711/1‑134 (MQ=255)
                                                                                                 gAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCTCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATACCGCtctcaa                            >  5:34221/1‑132 (MQ=255)
                                                                                                   aCATCTTGTTTTCTGCCACCGAAAACATTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATa                                      >  8:31928/1‑122 (MQ=255)
                                                                                                   aCATCTTGTTTTCTGCCACCGAAAACATTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATa                                      <  3:264109/122‑1 (MQ=255)
                                                                                                   aCATCTTGTTTTCTGCCACCGAAAACATTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATa                                      >  4:264109/1‑122 (MQ=255)
                                                                                                   aCATCTTGTTTTCTGCCACCGAAAACATTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATa                                      <  7:31928/122‑1 (MQ=255)
                                                                                                   aCATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATACCGCTCTcaacaa                         >  1:111606/1‑135 (MQ=255)
                                                                                                   aCATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATACCGCTCTCCACaa                         >  4:121491/1‑135 (MQ=255)
                                                                                                   aCATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATACCGCTCTCCACaa                         >  1:237180/1‑135 (MQ=255)
                                                                                                                   cACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATACCGCTCTCCACAAGCCTCCAAAGAcgccg         >  6:240265/1‑135 (MQ=255)
                                                                                                                   cACCGAAAACACTCGCATGCGGCCCCTTTTGTAGTCTTCCCGCGCACCCGGACGGCTCAGGCAAAATGCCGCCACATGGGCCTCGCTCATCTTCCTTTCCGACATACCGCTCTCCACAAGCCTCGAATGAcgcc          <  2:76751/134‑1 (MQ=255)
                                                                                                                    aCCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCGTCCGTTCAGACATACCGCTCTCCACAAGCCTCGAATGAcgccgc        >  4:57196/1‑135 (MQ=255)
                                                                                                                         aaaCACTCGGAGGCCGCCCCTTTTGTAGGCTTCCCGCGCACCCGGCAGGCGCATGCAAAATGCCGCCACCTGGCCCGCGCTCATCTTCCTTTCAGACATACCGCTCTCCACACGCCTCCAATGACGCCGCCAggt   <  4:354398/135‑1 (MQ=255)
                                                                                                                          accacCCGGAGCCCGCCCCATTGCTCGGCACCCCGCCCCCCCGGCAGGCGCCGCCAAACGCCCCCCCCCGGGCCCGCGCGCATCTTCCGTTCAGACAACCCGCTCCCCACACGCCTCTAATGACGCCGCCAggtg  <  4:377856/133‑1 (MQ=255)
                                                                                                                              |                                                                                                                                  
ACGGCGCCGGCCGTATCCCGGGTCGATCGCAATAACCCAGAAACAGCTCGTCGGCTACCTTGAACGACAACCCGGCACCCGCCAGGTCGAGCACCGCGAACATCTTGTTTTCTGCCACCGAAAACACTCGGATGCCGCCCCATTTGTAGTCTTCCCGCGCACCCGGCAGGCTCAGGCAAAATGCCGCCACCTGGGCCTCGCTCATCTTCCGTTCAGACATACCGCTCTCCACAAGCCTCGAATGACGCCGCCAGGTG  >  KT2440/3259339‑3259595

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: