Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,721,870 A→C 36.4% V98G (GTC→GGC)  CJ019_01617 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,721,8700AC36.4% 13.9 / 2.3 30V98G (GTC→GGC) CJ019_01617hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/2);  new base C (8/2);  total (26/4)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCT  >  KT2440/1721741‑1721992
                                                                                                                                 |                                                                                                                          
tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCGGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGCCCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGAccccg                                                                                                                       >  8:217002/1‑133 (MQ=255)
tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTGCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATAGTGGTCCCCGTCTTACTCGGCCCAg                                                                                                                       >  2:326616/1‑135 (MQ=255)
tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGCTGCGGCCCCCGCGGGCAGCATCTTCTTCCCCGTCTTCCTCGCCCCAg                                                                                                                       >  3:259059/1‑135 (MQ=255)
tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAg                                                                                                                       >  2:287465/1‑135 (MQ=255)
tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCCGCGCTGATGCGGACCCAGCGGGCACCATCTTAATCCCCGGATTCCTCGcccccc                                                                                                                       >  5:14133/1‑133 (MQ=255)
       tAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGGCCCAGCGGTCAGCATCTTCTTCCCCGGCTTCCTCGCCCCAGCCACCCg                                                                                                                >  7:167184/1‑135 (MQ=255)
         aaCCAATTCATAAATAAACGAGCAAGCGTGCGGTTTAATTATGTACGCCAGGTTATTTCCTCTTCCCCGTCGGCGCTGATGCGGACCCAGCGGTCAGAATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGgc                                                                                                              <  2:1279/135‑1 (MQ=255)
                        aaaCGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCCCCCGGCGCGCAGCGCCCTTcc                                                                                               >  6:152645/1‑135 (MQ=255)
                          aCGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGGCTTCCTCGACCCAGCCCCCCCGCGCGCAGCGCACTTCcac                                                                                             >  6:124973/1‑135 (MQ=255)
                                aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGCCCCAGCGGGCAGCCTCTTCGTCCCCGTCTTCCTCGCCCCAGCCACCCGGCGCGCCGCGCCCGTcccccccccg                                                                                       >  4:25644/1‑133 (MQ=255)
                                aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCGCGGGGCGCCTCTTCGTCCCCGTCGTCCTCGACCCCCCCCCCCCGCGCCCCGCGCCCATCCACCCCCac                                                                                       >  2:111285/1‑134 (MQ=255)
                                aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGCCCCAGCCACCCGGCGCGCAGCGCACTTCCCCCCCCAg                                                                                       >  4:127903/1‑135 (MQ=255)
                                aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCCGGCGCGCAGCGCACTTCCACcccccg                                                                                       >  8:251367/1‑133 (MQ=255)
                                aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTAGCCCCAGCCACCCGGCGCCCAGCGCCCTTACACACCCac                                                                                       >  1:7143/1‑134 (MQ=255)
                                aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAACTTATTCCCCGGCGTCCGCGCACCAGCCACCCGGCGCGCCCCGCACGTACACACCCac                                                                                       >  5:333682/1‑134 (MQ=255)
                                     tGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCAGATGCGGACCCCGCGGGCCGCATCGTCTTCCCCGGCTTCCGCGGCCCCGCCCCCCGGCGCGCCGCGCCCGTCCCCACCCCgcgccg                                                                                  >  1:77375/1‑133 (MQ=255)
                                          tttAATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCCCGGTCCGCATCTTCTTCCCCCTCTTCCTCGACCCCGCCACCCCGCGCGCAGCGCCCATCCACACCCACCGCGGCcaaa                                                                             >  6:36836/1‑134 (MQ=255)
                                               tgatgTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCCTCTTCTTCCCCGTCGTCCTCGCCCCAGCCCCCCGGCGCGCCGCGCAAATCCACACCCCGCCCCGCCAATgcggc                                                                        >  5:158950/1‑135 (MQ=255)
                                                               aTTTCCTCTTCACCGTCTGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTACTCGACCCAGCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAAcc                                                        >  6:80186/1‑134 (MQ=255)
                                                               aTTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGGCAGCAACTTCTTCCCCGTCTTCCGCGCCCCAGCCCCCCGGCGCGCAGCGCACTTACACACCCAGCGCGCCCAATGCGGCGCCGGCGCAGGCCACg                                                        >  6:75389/1‑135 (MQ=255)
                                                                                 gcgcTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCCCGGCGTc                                      >  1:13280/1‑135 (MQ=255)
                                                                                    cTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCCCGGCGTCTgg                                   >  8:66803/1‑135 (MQ=255)
                                                                                               cccgCGGCCAGAATCTTCTCCCCCGGCTCCCCCCCCCCACCCCCCCGGCGCCCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAg                        <  4:259059/132‑1 (MQ=255)
                                                                                               ccAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCACCGTCATCGGCCCACGGCGTCTGGGCGCgtt                            >  4:274417/1‑131 (MQ=255)
                                                                                                  gCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGAt                      >  4:322832/1‑134 (MQ=255)
                                                                                                   cGGTCTGCATGTTATCCCCCGGCTCGCGCGACCCAGCCCCCGGCCCCCCCGCCCACTTCCACAGCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATAc                    <  4:307908/135‑1 (MQ=255)
                                                                                                    ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACa                   >  4:299043/1‑135 (MQ=255)
                                                                                                    ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACa                   >  2:214777/1‑135 (MQ=255)
                                                                                                     gTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCCGGCAACGTCATCGGCCCACGGCGGCGGGTCGCGGTGCAGATACac                  >  3:6070/1‑134 (MQ=255)
                                                                                                      acctcATCTTCTCCCCGGGCTGCCCCCCCCCCGCCACCGGGGGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAgg                 <  3:144603/131‑1 (MQ=255)
                                                                                                                       gTCTTCCTCGACCCAGCCACCCGGCGCGCAGCCCACTTCCACACCCAGCGCCGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGCCGTCTGGTCGCTTTACATATACAGGCTGTTCCATTTACCt  >  6:172272/1‑133 (MQ=255)
                                                                                                                                 |                                                                                                                          
TCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCT  >  KT2440/1721741‑1721992

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