Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,021,386 C→A 7.9% intergenic (+203/‑46) hisC_6 → / → CJ019_02839 Histidinol‑phosphate aminotransferase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,021,3860CA 7.9% 136.9 / 4.0 51intergenic (+203/‑46)hisC_6/CJ019_02839Histidinol‑phosphate aminotransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/18);  new base A (2/2);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GCAAAAGGGTTATGCCCTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGACCCTGCTGTTTTCGACCGCCAGTGCCTG  >  KT2440/3021258‑3021520
                                                                                                                                |                                                                                                                                      
gCAAAAGGGTTATGCCCTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCTGCACACCCACCCTCgta                                                                                                                                   <  4:261537/134‑3 (MQ=255)
gCAAAAGGGTTATGCCCTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAat                                                                                                                                  <  8:235261/135‑1 (MQ=255)
 cAAAAGGGTTATGCCCTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAatc                                                                                                                                 >  8:210193/1‑135 (MQ=255)
    aaGGGTTATGCCCTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAatcatc                                                                                                                              <  7:210193/135‑1 (MQ=255)
            tGCCCTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGaa                                                                                                                      >  7:131194/1‑135 (MQ=255)
               ccTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACCCCCACCCTCGAATCATCGATCAGAAgtg                                                                                                                   >  2:31256/1‑135 (MQ=255)
               ccTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAgtg                                                                                                                   >  1:16076/1‑135 (MQ=255)
               ccTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAgtg                                                                                                                   >  7:167054/1‑135 (MQ=255)
                    aGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAgtggtg                                                                                                                <  8:3913/133‑1 (MQ=255)
                    aGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGcc                                                                                                              <  2:207296/135‑1 (MQ=255)
                    aGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGc                                                                                                               >  2:166474/1‑134 (MQ=255)
                     gggTTTCGAGGCCTCTGGGTCTCAACGCAATCCGTTCGGTTACCGCACATAACCGCATCGCCGCGCCTCAAGCCACCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCCATCAGAAGtgctgcac                                                                                                             >  6:36132/1‑133 (MQ=255)
                           cGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCcaacaa                                                                                                         >  2:138912/1‑133 (MQ=255)
                           cGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAg                                                                                                        <  8:167054/134‑1 (MQ=255)
                            gAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCCCCCTCGAATCATCGATCAGAAGTGGTGCCAACAAgca                                                                                                      >  1:130624/1‑133 (MQ=255)
                            gAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACAGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGa                                                                                                      >  4:156762/1‑135 (MQ=255)
                                    tGGGTCTAAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAACCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTAGAATCATCGAACAGAAGTGGTGCCAACAACGACTACTGcc                                                                                              >  2:117829/1‑135 (MQ=255)
                                     gggTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGc                                                                                               <  2:16076/133‑1 (MQ=255)
                                     gggTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTACGATAGCAGCACACCCACCCTCGAATCATCGATCACAAGTGGTGCCTACAATGACTACTGccc                                                                                             >  5:123051/1‑135 (MQ=255)
                                             aCGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCtt                                                                                     >  1:15605/1‑135 (MQ=255)
                                              cGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCTCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCAAACAAGGAGTACTGCCCGCATGCtt                                                                                     >  8:146403/1‑134 (MQ=255)
                                              cGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTg                                                                                    <  5:262251/135‑1 (MQ=255)
                                               gCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTg                                                                                    >  6:148146/1‑134 (MQ=255)
                                               gCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTg                                                                                    >  5:74867/1‑134 (MQ=255)
                                               gCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTg                                                                                    >  8:231943/1‑134 (MQ=255)
                                                   tCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAAc                                                                               >  7:60889/1‑135 (MQ=255)
                                                             tACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTc                                                                     <  5:225157/135‑1 (MQ=255)
                                                             tACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTc                                                                     <  3:113283/135‑1 (MQ=255)
                                                                  cacaTAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCATCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCAcc                                                                >  2:176210/1‑135 (MQ=255)
                                                                  cacaTAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCATCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCAcc                                                                <  1:176210/135‑1 (MQ=255)
                                                                    caTAACCGCATCGCGGCGCGTCAAGCCTGCGCCAGTACCTCCGATAGGAGCACACCGACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCgg                                                              <  1:31256/135‑1 (MQ=255)
                                                                               cgcgGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTggg                                                   >  2:293958/1‑135 (MQ=255)
                                                                               cgcgGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTggg                                                   >  4:166106/1‑135 (MQ=255)
                                                                               cgcgGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTggg                                                   >  8:188954/1‑135 (MQ=255)
                                                                               cgcgGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCATCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGtt                                                        <  4:13840/130‑1 (MQ=255)
                                                                               cgcgGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCATCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGtt                                                        >  3:13840/1‑130 (MQ=255)
                                                                                  ggcgCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGgtgt                                                >  3:43759/1‑135 (MQ=255)
                                                                                    cgcCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACCCCCACCCTCGAATCATCGATCAGAAGTGGTGCCAGCAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGtgttg                                              >  1:18954/1‑135 (MQ=255)
                                                                                    cgcCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCAcc                                                                <  5:121412/117‑1 (MQ=255)
                                                                                    cgcCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCAcc                                                                >  6:121412/1‑117 (MQ=255)
                                                                                    cgcCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGtgttg                                              <  6:55500/135‑1 (MQ=255)
                                                                                           aGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCt                                         <  1:293958/133‑1 (MQ=255)
                                                                                              cTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCa                                    <  1:166474/135‑1 (MQ=255)
                                                                                                    cAGTACCTCCGATAGCAGCACCCCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGa                              >  7:259983/1‑135 (MQ=255)
                                                                                                    cAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCcaac                                                                                                           <  3:180380/58‑1 (MQ=255)
                                                                                                    cAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCcaac                                                                                                           >  4:180380/1‑58 (MQ=255)
                                                                                                    cAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGa                              >  5:92806/1‑135 (MQ=255)
                                                                                                    cAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGa                              >  7:137932/1‑135 (MQ=255)
                                                                                                    cAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGa                              >  5:112246/1‑135 (MQ=255)
                                                                                                                     gcaCACCCACCCTCGACTCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTGCTGGGTGTTGTCGCTGTTCAGCCTGACCCTGCTGGTTTCGAcc             <  2:163541/135‑1 (MQ=255)
                                                                                                                      cacaCCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGACCCTGCTGTTTTCGACCg            <  5:148146/135‑1 (MQ=255)
                                                                                                                         acccacccTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGACCCTGCTGTTGTCGACCGCCa         >  6:207350/1‑135 (MQ=255)
                                                                                                                                cTCGAATCATCGATCAGAAGTGGTGTCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGACCCTGCTGTTTTCGACCGCCAGTGCCtg  >  2:132300/1‑135 (MQ=255)
                                                                                                                                |                                                                                                                                      
GCAAAAGGGTTATGCCCTAAAGGGTTTCGAGGCCTCTGGGTCTCAACGCATTCCGTTCGGTTACCGCACATAACCGCATCGCGGCGCCTCAAGCCTGCGCCAGTACCTCCGATAGCAGCACACCCACCCTCGAATCATCGATCAGAAGTGGTGCCAACAAGGAGTACTGCCCGCATGCTTGCAAACCTGAAAATTCGCACCGGAATGTTCTGGGTGTTGTCGCTGTTCAGCCTGACCCTGCTGTTTTCGACCGCCAGTGCCTG  >  KT2440/3021258‑3021520

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