Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,180,021 T→G 24.4% E136A (GAG→GCG)  aceK ← Isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,180,0210TG24.4% 41.5 / 5.0 41E136A (GAG→GCG) aceKIsocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (3/26);  new base G (6/5);  total (10/31)
Fisher's exact test for biased strand distribution p-value = 6.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGA  >  ACB122/5179891‑5180153
                                                                                                                                  |                                                                                                                                    
ttcttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTcgcg                                                                                                                                  >  7:304487/1‑135 (MQ=255)
       aGCCGCCCCAGGGCGCCTTCCAGGTCGCCATAGGGCATGTCGACCGGGTAAGCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAACCCGATAGGTCCGGGTTGGCGCGGCGCGCTCGCGGCCgcgc                                                                                                                           <  8:71029/135‑1 (MQ=255)
          cgcgCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCGCGCGGCCgcgccgg                                                                                                                        >  4:371392/1‑132 (MQ=255)
            cgcCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGGGCGGCGCGCGCGCGGGcgcggcgggg                                                                                                                      >  7:318001/1‑129 (MQ=255)
            cgcCAGGGCGCCTTCCAGGGCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGa                                                                                                                      <  7:97506/135‑1 (MQ=255)
                   gcgcgTTCCCGGTCGCCATAGGGCTGGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCTATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCgg                                                                                                               <  8:318001/134‑1 (MQ=255)
                      ccTGCCAGGTCGCCATAGGGAATGTCGACCGGGTAATCGGCGAACACCGCGCGAAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTa                                                                                                            <  6:339417/135‑1 (MQ=255)
                            aGGTCGCCATCGGGCATGTCGACCGGTTAAGCGGCCAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTgg                                                                                                      <  3:356108/135‑1 (MQ=255)
                            aGGGCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACGGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCGCGCGGCCGCGCATGGAAGGCCGGGTAGTGTgg                                                                                                      <  2:30743/135‑1 (MQ=255)
                             ggTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGa                                                                                                     <  2:208597/135‑1 (MQ=255)
                                      aGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGc                                                                                            <  8:304487/135‑1 (MQ=255)
                                      aGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTAAAGCCGATTGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGc                                                                                            <  8:109863/135‑1 (MQ=255)
                                           aTGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCg                                                                                       <  8:32056/135‑1 (MQ=255)
                                                gAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCGCGCGGCCCGGCCTGGGAGGCCGGGTTGTGGGGTTGTACAATCATCCGtgcca                                                                                  >  4:348384/1‑131 (MQ=255)
                                                 accGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGTGAGTGCCGTCAAGCCGATAGGTAGGGGTTGGCGCGGCGCGCGCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCt                                                                                  <  1:305189/132‑1 (MQ=255)
                                                     ggtatgcGGCCAACACGGGGCGCACCAAGTTCGGGGGAGTGCCGTCATGCCCATTGGTAGGGTTTTGCGCGGCGCGCTCGCGGCCGCGGATGGAAGGCCGGGTTGTGTGGATGAACATGCAGCCGTCGCTGATCa                                                                             <  3:15771/129‑1 (MQ=255)
                                                     ggTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGAGTGCCGTCAAGCCGATAGGTTCGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCa                                                                             <  3:274757/135‑1 (MQ=255)
                                                       tAATCGGCTAACACCGCGTGCTGCATGTTCCTGAGACTGCGGTCATGCGGTTTGTTACGGGTGTGGGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAgg                                                                           <  3:154526/135‑1 (MQ=255)
                                                       tAATCGGCTAACACCGCGCGCAGCAAGTTCGTTTGACTGCCTTCACGCCGTTAGGTACGGGTTTGCGCGGCGCGCGCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAgg                                                                           <  8:331440/135‑1 (MQ=255)
                                                       tAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAgg                                                                           <  1:72859/135‑1 (MQ=255)
                                                           cGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCGCGCGGCCGCGGACGGGAGGGCGGGGAGGGGGGTTGGAAATGCAGCCGTGGCGGCTCAggggat                                                                       >  7:240570/1‑131 (MQ=255)
                                                                 acacCGGGCGCAGCAAGTTCCTGAGACTCCCGTCAAGCCGATAGGTACGGGTTGGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGa                                                                 <  7:89728/135‑1 (MQ=255)
                                                                 acacCGCGCGCTCCAAGTTCCTTGGCAGGCCGTCAACCCGATCGGAACGGGTTGGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGa                                                                 <  6:970/135‑1 (MQ=255)
                                                                   acCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCGCGCGGCCCCGCATGGGAGGGCGGGGAGTGTGGATTGACATGCAGCCCTCGCTGATCAGGTCATGGCTgc                                                                 >  8:337331/1‑132 (MQ=255)
                                                                              caAGTTCCTGAGAATGCCGGCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGGAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaaaa                                                    <  6:94651/135‑1 (MQ=255)
                                                                               aaGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAAc                                                   <  3:348384/135‑1 (MQ=255)
                                                                                  ttCCTGAGACTGCGGGCTAGCCGTGAGGTACGGGTTTGCGCGGGGCGCTCGCGGCCGCGAATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGc                                                <  2:87936/135‑1 (MQ=255)
                                                                                  gtCCTGAGAGTGACTTCAAGCCGATTGGTAGGGGTTTGCGCGGCGCTCTCGCGGCGGCGCATGGAAGGCCGGGTAGTGTGGATGAAGATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTTCAAAACAGc                                                <  7:381083/134‑1 (MQ=255)
                                                                                   tCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCg                                               <  3:181934/135‑1 (MQ=255)
                                                                                         gaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGt                                         <  6:136474/135‑1 (MQ=255)
                                                                                         gaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCACGCGGCCGCCCATGGAAGGCCGCGGTGTGGGGATTCTCATCACGCCCTCGCAGATACTGTCAAGGCTGCAGAGGGTGGAAAACACCGAGgtgt                                         >  2:21829/1‑135 (MQ=255)
                                                                                             gccgccAAGCCGATTGGTCCGGGTTGGCGCGGCGCGCGCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCa                                     <  1:155490/135‑1 (MQ=255)
                                                                                               ggTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAgg                                   <  3:340169/134‑1 (MQ=255)
                                                                                                 tCAAGCCGATAGGTACGGGTTTGCCCGGCGCGCGCCCGGCCGCGCCTGGGAGGCCGGGGTGTGTGGGTGTACAAACAGCCCGGGCGGGTCGCGTGATGGGGGAAGAGGAGGCAAAAAACCGAGGTGTACCACGtc                                 >  2:353385/1‑135 (MQ=255)
                                                                                                    aGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACAAGGTctct                              <  3:172411/135‑1 (MQ=255)
                                                                                                          tcggTCCGGGTTGGCGCGGCGCGCGCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAAc                        <  7:145789/133‑1 (MQ=255)
                                                                                                               aCGGGGTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtcgtc                   <  8:240570/135‑1 (MQ=255)
                                                                                                                              gcgcTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTc    <  4:33354/135‑1 (MQ=255)
                                                                                                                              gcgcTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTc    <  4:169677/135‑1 (MQ=255)
                                                                                                                                gcTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCGGATCCGGTCATGGCTGAAGAGGGTGCAAAAACGCGAGTTGTTCCAGGTCTCTGCCAACTCGGCGTCCAGGCGGGGGTCgc  >  7:200722/1‑134 (MQ=255)
                                                                                                                                gcTAGCGGCCGCGCATGGAAGGCCGCGTATTGTGGATGAACATCCAGCCGTCGCTGATCAGGTCATTGCTGAAGAGGCTGCAAAACAGCCAGTTGTACCACGGCTCTGCCAAATAGTCGTACACGCCGGGGTCg   >  2:54692/1‑134 (MQ=255)
                                                                                                                                  |                                                                                                                                    
TTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGA  >  ACB122/5179891‑5180153

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