Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,481,837 T→C 8.0% L114P (CTG→CCG)  CJ019_01397 → GTP cyclohydrolase 1 type 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,481,8370TC 8.0% 145.3 / 4.2 51L114P (CTG→CCG) CJ019_01397GTP cyclohydrolase 1 type 2
Reads supporting (aligned to +/- strand):  ref base T (21/25);  new base C (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

GGAAGGGCGAGAACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGCAACAGGTGCTGGGGC  >  KT2440/1481705‑1481971
                                                                                                                                    |                                                                                                                                      
ggAAGGGCGAGAACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTgg                                                                                                                                      >  1:143070/1‑135 (MQ=255)
  aaGGGCGAGAACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGgcg                                                                                                                                    <  6:72069/135‑1 (MQ=255)
   aGGGCGAGAACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGgcgc                                                                                                                                   <  1:10465/135‑1 (MQ=255)
   aGGGCGAGAACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGgcgc                                                                                                                                   >  2:94167/1‑135 (MQ=255)
         agaACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGAc                                                                                                                             >  6:52283/1‑135 (MQ=255)
          gaACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACt                                                                                                                            <  2:143070/135‑1 (MQ=255)
            aCCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGa                                                                                                                          <  2:214907/135‑1 (MQ=255)
               cGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGa                                                                                                                       <  3:210914/135‑1 (MQ=255)
               cGTGCATCACCGGTATTCGCCAGCGGCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGa                                                                                                                       <  6:175915/135‑1 (MQ=255)
                   cATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGAt                                                                                                                      >  3:53562/1‑132 (MQ=255)
                   cATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGAt                                                                                                                      <  4:53562/132‑1 (MQ=255)
                   cATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCAGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGCCAACAAAGAGCAGCTGGCGCGCCGACTGAAGATc                                                                                                                     >  6:270501/1‑133 (MQ=255)
                        ccGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCg                                                                                                              <  2:126658/135‑1 (MQ=255)
                           gTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGa                                                                                                                       >  3:79006/1‑123 (MQ=255)
                           gTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGa                                                                                                                       <  4:79006/123‑1 (MQ=255)
                           gTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAg                                                                                                           >  4:141024/1‑135 (MQ=255)
                                         cGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGCTCACTGTCGAAGGGCCGCTAGCCCCgg                                                                                            >  8:232383/1‑136 (MQ=255)
                                         cGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCg                                                                                             >  8:163961/1‑135 (MQ=255)
                                                aGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAAcc                                                                                      <  3:6836/135‑1 (MQ=255)
                                                  aCCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGCCGTGCAACCGGACGTTGGCAAAAATGTTCAGCTGGCGCGGCGACAGAAGAACCCTGTCCAAGCTCCGCAAGAACCGGAAAAccc                                                                                     >  5:149167/1‑134 (MQ=255)
                                                        cTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGACGGGCCGCTAGACCCGGAAAACCCCAAAg                                                                                <  4:76752/133‑1 (MQ=255)
                                                         tGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGgcgc                                                                                                                                   >  8:268771/1‑81 (MQ=255)
                                                         tGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGgcgc                                                                                                                                   <  7:268771/81‑1 (MQ=255)
                                                         tGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTgg                                                                             >  5:239965/1‑135 (MQ=255)
                                                                        gCCTGTTGGCGTTCCTCCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGtt                                                              <  6:149167/135‑1 (MQ=255)
                                                                        gCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGtt                                                              >  5:82743/1‑135 (MQ=255)
                                                                        gCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCCGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTtc                                                             >  7:125787/1‑136 (MQ=255)
                                                                        gCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCCGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCGAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATGTGTtc                                                             >  7:122399/1‑136 (MQ=255)
                                                                             ttGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCAGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGAATGGTtcgctc                                                         <  7:163961/135‑1 (MQ=255)
                                                                              tGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTtcgctc                                                         <  4:180760/134‑1 (MQ=255)
                                                                                gCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGGAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGcc                                                      >  5:113424/1‑135 (MQ=255)
                                                                                gCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTtcgctcgc                                                       >  2:25138/1‑134 (MQ=255)
                                                                                gCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGcc                                                      >  6:219565/1‑135 (MQ=255)
                                                                                       aCCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCGGATTGGTTCGCTCGCCGCATccc                                               >  2:102199/1‑134 (MQ=255)
                                                                                                    cGTGCACCCGGACGTGGGCAACAATGTGCAGCAGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTTATTGGTTCGCTAGCCGAAACCATGACcccgcgcgc                                  >  5:16591/1‑134 (MQ=255)
                                                                                                         aCCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTg                             >  5:677/1‑135 (MQ=255)
                                                                                                         aCCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGc                            >  7:71517/1‑136 (MQ=255)
                                                                                                           ccGGACGTGGGCAACAAAGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGccc                          <  7:19507/136‑1 (MQ=255)
                                                                                                               aCGTGGGCAACAATGTGCAGCCGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCccgccg                      <  8:125787/136‑1 (MQ=255)
                                                                                                               aCGGGGGCAACAATGGGCAGCCGGCGCGGCGACTGAAGATCACTGTCGCAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCccgccg                      <  8:122399/136‑1 (MQ=255)
                                                                                                                cGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCccgccg                      >  4:172062/1‑135 (MQ=255)
                                                                                                                     gCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGc                 <  1:27222/135‑1 (MQ=255)
                                                                                                                     gCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGc                 <  1:25138/135‑1 (MQ=255)
                                                                                                                     gCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGc                 <  6:3111/135‑1 (MQ=255)
                                                                                                                     gCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGAGATTTTGCCCGCCGGGTGc                 <  6:113424/135‑1 (MQ=255)
                                                                                                                      caacaaTGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGGGGGCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGCaa               <  7:157870/136‑1 (MQ=255)
                                                                                                                           aTGGGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGGCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGCAACAGGt          <  7:73172/136‑1 (MQ=255)
                                                                                                                             gtgCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGCAACAGGTg         <  4:108186/135‑1 (MQ=255)
                                                                                                                                 aGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAg                                                                                                  >  6:34540/1‑42 (MQ=255)
                                                                                                                                 aGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAg                                                                                                  <  5:34540/42‑1 (MQ=255)
                                                                                                                                  gCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGCAACAGGTGCTggg    <  6:82743/135‑1 (MQ=255)
                                                                                                                                    tGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGCAACAGGTGCTGGGGc  >  4:95653/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                      
GGAAGGGCGAGAACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCCAAAGTGGTCGGGCTGATTGGTTCGCTCGCCGAATCCATGACCGCGCGCGATTTTGCCCGCCGGGTGCAACAGGTGCTGGGGC  >  KT2440/1481705‑1481971

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