Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,072,630 T→G 28.3% intergenic (‑308/+126) CJ019_02884 ← / ← yabJ_2 putative protein/2‑iminobutanoate/2‑iminopropanoate deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,072,6300TG28.3% 55.5 / 2.4 41intergenic (‑308/+126)CJ019_02884/yabJ_2putative protein/2‑iminobutanoate/2‑iminopropanoate deaminase
Reads supporting (aligned to +/- strand):  ref base T (10/17);  new base G (11/2);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 6.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCTGGGGATGCAGGGAATGTGAGGGGCACACCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAG  >  KT2440/3072498‑3072759
                                                                                                                                     |                                                                                                                                 
gTTCTGGGGATGCAGGGAATGTGAGGGGCACACCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAAAGGGGCGCAGGCCGGTGCCAgg                                                                                                                                 >  8:97363/1‑135 (MQ=255)
     ggggATGCAGGGAATGTGAGGGGCACACCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGGGTGAGGCGCAAAGGGGCGCACGCCGGTGGCGggggggc                                                                                                                            >  2:293001/1‑134 (MQ=255)
             agtgAATGTGAGGGGCACACCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAACAGGGAGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGaaaaa                                                                                                                      <  3:160253/130‑1 (MQ=255)
                 aaTGTGAGGGGCACACCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGt                                                                                                                <  8:296127/135‑1 (MQ=255)
                   ggggAGGGGCACACCCGAGCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAACGGGCGCTTTTGACCCGGGCGCTGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTAt                                                                                                              <  6:163913/132‑1 (MQ=255)
                              cacCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGTTTGAACCGAAATGGGGCGCAGGCCGGTGCCGGGCGGGAGAAAAAGGGAGTATTCCGGGTGCAt                                                                                                   >  4:191194/1‑135 (MQ=255)
                              cacCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGTTTGAACCGAAATGGGGCGCAGGCCCGGGGCGGGGGGGAGCAAAAGGGAGTATTACGGGTGCAt                                                                                                   >  4:170967/1‑135 (MQ=255)
                               acCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATg                                                                                                  <  8:221367/135‑1 (MQ=255)
                                        gtgCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGTTTGTAACCAAAATGGGCGCCGGCCGGTGCCGGGCGGGTGAAAAAGGGAGTATTGAGGGTGCATTACGTggggg                                                                                         >  4:185151/1‑134 (MQ=255)
                                          gctgcCTGGCCGCAAGGTCCGGCGCATCAATAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGa                                                                                      <  4:90392/135‑1 (MQ=255)
                                          gctgcCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTgggg                                                                                          <  4:21947/132‑1 (MQ=255)
                                          gctgcCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCg                                                                                        <  6:208093/134‑1 (MQ=255)
                                              ccTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGTTTGCCCCGCAATGGGGCGCCGGCCGGTGCGGGGCGGGTGAAAAAGGGAGTATTCCGGGTGCTTGAcggggggcgggggg                                                                                   >  1:163858/1‑135 (MQ=255)
                                                tGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAggggg                                                                                 <  7:113952/135‑1 (MQ=255)
                                                 ggCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTGTTGGACCGAAATGGGGCGCCGGCCGGTGCCGGGGGGGAGAAAAAGGGAGTATTGCGGGTGCATGCCGTGGGGGGGAGGGGgc                                                                                >  1:172694/1‑135 (MQ=255)
                                                 ggCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGTTTGACCCGATAGGGGGCACAGGCCGGTGCCGGGCGGGGGAACAACGGAGTATTCCGGGTGGATGCCGGGGGGCgggggggg                                                                                 >  1:172689/1‑134 (MQ=255)
                                                 ggCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGGTTGAGCCCCAATGGGGCGCCGGCCGGGGCCGGGCGGGAGCAAAAGGGAGGTTTCCGGGGGCCAGACGGGGGGCGGCgggg                                                                                  >  3:23612/1‑133 (MQ=255)
                                                      cAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGTTTGTACCGCAATGGGGCGCCGGCCGGTGCCGGGGGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGGGGGGCGGAGGGGgcatgca                                                                           >  5:132588/1‑135 (MQ=255)
                                                           gCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGGTTTGGACCGGAATGGGGCGCCGGGCGGGGCCGGGCGGGGGCAAAAGGGAGTTTTCCGGGTGCATTCCGGGGGGGGGGGGGGggctgcatagcg                                                                      >  8:132628/1‑131 (MQ=255)
                                                               gcgcATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGGCGTGGGGCGGAGGGGGCATGCATCGAGcgcc                                                                  >  4:188524/1‑135 (MQ=255)
                                                                 gcATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGGCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCCTGACGTGGGGCGGAGGGGGCATGCATCGAGcgccgc                                                                >  6:252772/1‑135 (MQ=255)
                                                                   aTCAAAACCGGGCGCTTTGGCCCCGGGCGCGGGTTTTGACCCGAAAGGGGGCGCCGGCCGGGGCCGGGCGGGAGAAAACGGGAGAATTCCGGGGGCATGACGCGGGGCGGAGGGGGCATGCAGCGAGCGCcgcgc                                                              <  2:299980/135‑1 (MQ=255)
                                                                   aTCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCcgcgc                                                              <  7:97363/135‑1 (MQ=255)
                                                                   aTCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCGGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCcgcgc                                                              <  3:191194/135‑1 (MQ=255)
                                                                        aaaaGGGCGCTTTTGACCCGGGCGCGGGTTTTGCACCGAAATGGGGCGCCGGCCGGGGGCTGGGGGGAGGAAAACGGAGGATTCCGGGTGGATGACGGGGGGCGCAGGGGGCCTGCATACAGCGCCGCccgggcg                                                         >  5:85109/1‑135 (MQ=255)
                                                                        aaaaGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGCCGTGGGGCGGAGGGGGCATGCCTCGAGCGCCGCgcgggcg                                                         >  6:117919/1‑135 (MQ=255)
                                                                        aaaaGGGCGCTTTTGACCCGGGCGCGGGGTTTTAAACGAAATGGGGCGCAGGCCGGTGGCGGGCGGGAGAAAAACGGAGTATTACGGGTGCATGACGTGGGGCGGCGGGGGCATGCATCGAGCGCCGCgcgggcg                                                         >  5:254221/1‑135 (MQ=255)
                                                                                gcTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGgaga                                                                                                                          <  7:234053/62‑1 (MQ=255)
                                                                                gcTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGgaga                                                                                                                          >  8:234053/1‑62 (MQ=255)
                                                                                  ttttGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTgg                                                                                                                                 <  7:245867/53‑1 (MQ=255)
                                                                                  ttttGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTgg                                                                                                                                 >  8:245867/1‑53 (MQ=255)
                                                                                     tGACCCGGGCGCGGGTTTTGAACCGAAAGGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTg                                              <  5:17239/133‑1 (MQ=255)
                                                                                     tGACCCGGGCGCGGGGTTTGAACCGAAATTGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATTCATCGAGCGCCGCGCGGGCGGCACTCGATTgcg                                            >  3:101975/1‑135 (MQ=255)
                                                                                      gccccGGGCGCGGGTGTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAACAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTgcg                                            <  3:185151/132‑1 (MQ=255)
                                                                                      gACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTgcg                                            <  6:132588/134‑1 (MQ=255)
                                                                                                           aCCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTgaga                       >  7:136818/1‑135 (MQ=255)
                                                                                                                 aTGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGgcg                 >  3:142486/1‑135 (MQ=255)
                                                                                                                            gCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGAc                      <  6:293854/119‑1 (MQ=255)
                                                                                                                            gCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGAc                      >  5:293854/1‑119 (MQ=255)
                                                                                                                                gTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCag  <  6:115965/135‑1 (MQ=255)
                                                                                                                                gTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGGGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCag  <  5:206251/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                 
GTTCTGGGGATGCAGGGAATGTGAGGGGCACACCCGATCTGTGCTGCCTGGCCGCAAGGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAG  >  KT2440/3072498‑3072759

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: