Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,255,631 A→G 8.6% intergenic (‑73/‑23) metE ← / → CJ019_04952 5‑methyltetrahydropteroyltriglutamate‑‑homocystei ne methyltransferase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,255,6310AG 8.6% 134.5 / 4.0 47intergenic (‑73/‑23)metE/CJ019_049525‑methyltetrahydropteroyltriglutamate‑‑homocystei ne methyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/27);  new base G (2/2);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

AAATTTGAGATGAACCCAAGAAATCTTCATGATCATGCTGTCATGGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCTGGCCCGT  >  KT2440/5255499‑5255764
                                                                                                                                    |                                                                                                                                     
aaaTTTGAGATGAACCCAAGAAATCTTCATGATCATGCTGTCATGGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACca                                                                                                                                    >  4:237567/1‑136 (MQ=255)
        gaTGAACCCAAGAAATCTTCATGATCATGCTGTCATGGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGt                                                                                                                            <  3:243604/136‑1 (MQ=255)
         aTGAACCCAAGAAATCTTCATGATCATGCTGTCATGGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGt                                                                                                                            <  5:19938/135‑1 (MQ=255)
             aCCCACGAACTCTTCATGATCATGCTGTCATGGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGccc                                                                                                                        <  2:345616/135‑1 (MQ=255)
              cccAAGAAATCTTCATGATCATGCTGTCATGGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGccc                                                                                                                        <  2:334239/134‑1 (MQ=255)
                      aTCTTCATGATCATGCTGTCATTGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGAGATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATg                                                                                                              <  3:243691/136‑1 (MQ=255)
                             tGATCATTCTGTCATGGTCTTCGCGCCCCATGTCTCTTTCAGATCACCACCTTGCGCTTCCGGACATGACCCGGCTCGCGTTCTGGCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCct                                                                                                         <  6:289392/134‑1 (MQ=255)
                                 cATGCTGTCATGGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCctgctg                                                                                                     >  5:229271/1‑134 (MQ=255)
                                             gTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTg                                                                                         <  1:90123/134‑1 (MQ=255)
                                                ttCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGcgcg                                                                                     <  5:247131/135‑1 (MQ=255)
                                                     gcCCCATGTCTCATTCAGCTCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAg                                                                                >  4:92769/1‑135 (MQ=255)
                                                     gcCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGGCCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAg                                                                                >  3:150280/1‑135 (MQ=255)
                                                     gcCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGGCCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAg                                                                                >  3:156370/1‑135 (MQ=255)
                                                     gcCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACa                                                                                 >  5:323473/1‑134 (MQ=255)
                                                         cATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCGGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACtt                                                                           <  3:237169/136‑1 (MQ=255)
                                                         cATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACtt                                                                           <  4:235744/136‑1 (MQ=255)
                                                                         caccacCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGa                                                            <  6:128283/135‑1 (MQ=255)
                                                                          accaccTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAg                                                           <  7:152454/135‑1 (MQ=255)
                                                                          accaccTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCACGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAAGGGAg                                                           >  1:108336/1‑135 (MQ=255)
                                                                                    gcATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTcccc                                                  >  3:179340/1‑134 (MQ=255)
                                                                                     cATCCGGACATGACCCTGCTCGCCTTCTGCCGCACCATCCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCaa                                                <  7:135645/135‑1 (MQ=255)
                                                                                      aTCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAg                                               >  7:51188/1‑135 (MQ=255)
                                                                                            aCATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTcc                                         >  8:174400/1‑135 (MQ=255)
                                                                                              aTGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGc                                       <  1:9985/135‑1 (MQ=255)
                                                                                                   ccGGCTCGCCTTCTGCCGCACAATGCCCCCATGCCCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTATTGCAGTTCGTCCCCCATCCGTCCGCTGacc                                  >  2:143646/1‑133 (MQ=255)
                                                                                                   ccGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAAc                                  >  7:149109/1‑135 (MQ=255)
                                                                                                   ccGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAAc                                  >  2:350768/1‑135 (MQ=255)
                                                                                                    cGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAg                                                           <  3:208328/109‑1 (MQ=255)
                                                                                                    cGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAg                                                           >  4:208328/1‑109 (MQ=255)
                                                                                                        tCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCaa                                                <  1:281056/116‑1 (MQ=255)
                                                                                                        tCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCaa                                                <  1:285841/116‑1 (MQ=255)
                                                                                                        tCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCaa                                                >  2:281056/1‑116 (MQ=255)
                                                                                                        tCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCaa                                                >  2:285841/1‑116 (MQ=255)
                                                                                                             ttCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTgccgc                        <  8:149109/135‑1 (MQ=255)
                                                                                                               cTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCa                      >  8:6200/1‑135 (MQ=255)
                                                                                                                  ccGCACAATGCCACCCTGGCCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATgg                   <  4:150280/135‑1 (MQ=255)
                                                                                                                  ccGCACAATGCCACCCTGGCCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATgg                   <  4:156370/135‑1 (MQ=255)
                                                                                                                  ccGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATgg                   <  4:182820/135‑1 (MQ=255)
                                                                                                                   cGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATCTGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATgg                   <  6:334501/134‑1 (MQ=255)
                                                                                                                   cGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATgg                   <  6:323473/134‑1 (MQ=255)
                                                                                                                   cGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATgg                   <  6:130304/134‑1 (MQ=255)
                                                                                                                            ccACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTGCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCt         <  3:92769/135‑1 (MQ=255)
                                                                                                                             cACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCt         <  7:52434/134‑1 (MQ=255)
                                                                                                                             cACCCAGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGATCTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCt         <  5:258966/134‑1 (MQ=255)
                                                                                                                                 cTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCTGGccc    <  7:288415/135‑1 (MQ=255)
                                                                                                                                  tGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCTGGccc    >  5:129982/1‑134 (MQ=255)
                                                                                                                                  tGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCTGGccc    <  6:129982/134‑1 (MQ=255)
                                                                                                                                    aCCCCACCAGGTGCCCGGCCAATCTGCGCCGGCTGAAAACCGCCCTGCGCGCACAGAACTTCGTCTGCGGACTGGAGTGCGTCCCCACGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCTGGCCCgg  <  4:343689/134‑2 (MQ=255)
                                                                                                                                    |                                                                                                                                     
AAATTTGAGATGAACCCAAGAAATCTTCATGATCATGCTGTCATGGTCTTCGCGCCCCATGTCTCATTCAGATCACCACCTTGCGCATCCGGACATGACCCGGCTCGCCTTCTGCCGCACAATGCCACCCTGACCACACCAGGTGCCCGGCCAATATGAGCCTGCTGAAAACCGCCCTGCGCGAACAGAACTTCGTCTGCGTAATGGAGTTCGTCCCCAAGCCTTCCGCAGAACGCTTTGCCGCCATGGAAGCGATCCTGGCCCGT  >  KT2440/5255499‑5255764

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