Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,646,123 T→C 9.8% S329P (TCC→CCC)  aceA → Isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,646,1230TC 9.8% 120.4 / 4.8 51S329P (TCC→CCC) aceAIsocitrate lyase
Reads supporting (aligned to +/- strand):  ref base T (28/17);  new base C (2/3);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCA  >  ACB122/4646001‑4646252
                                                                                                                          |                                                                                                                                  
ccccGTACTCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGga                                                                                                                       <  4:71185/135‑1 (MQ=255)
ccccGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTCCCTTCAAC‑Tg                                                                                                                         >  7:226892/1‑133 (MQ=255)
ccccGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTCCCTTCAAC‑Tg                                                                                                                         >  7:221272/1‑133 (MQ=255)
 cccGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGTACGAAGCCCGTCGCTTCGCAGAGGCGATCAAGAATGAGTTCCCGGACCAGAGACTGTCGTAAAACTGCTCGCCTTCCTTCAAC‑Tg                                                                                                                         >  5:323789/1‑132 (MQ=255)
  ccGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGgaa                                                                                                                      >  2:87719/1‑134 (MQ=255)
      aCGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGCCCAGATCCTGTCGTACCCCTGCGCGCCTTCCTTCCAC‑TGGTAGaaaa                                                                                                                 >  1:179238/1‑135 (MQ=255)
      aCGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGCCCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGaaaa                                                                                                                 >  1:276511/1‑135 (MQ=255)
      aCGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCCGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TCGAAGaaaa                                                                                                                 >  1:187480/1‑135 (MQ=255)
      aCGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGaaaa                                                                                                                 >  5:103145/1‑135 (MQ=255)
               tGATCTGGTGCGAAACCGCCAAGCCGGCCCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGTGTCCCCGGCCCACATCCTGTCGTACAACTGCTCGCCGTCCTTCAAC‑TCGAAGACATACCTGgacg                                                                                                        >  7:348419/1‑135 (MQ=255)
               tGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGgacg                                                                                                        >  2:129193/1‑135 (MQ=255)
                       tGCTAAACCGCCAAGCCGGACCTGGACTAAGCCCGTCGGTTGGCGGAGGCGTTCAAGAAGGAGTACCCGGACCAGATCCTGTCTTACAACTGCTCGCCTTCCTTCAAC‑TGGATGAAAAACCTGGACGACGcca                                                                                                  <  8:128979/133‑1 (MQ=255)
                          gAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGCCGCCACCATc                                                                                             >  4:1683/1‑135 (MQ=255)
                          gAAACCGCCAAGCCGCACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TCGAAGAAAAACCTGGACGACGccaccaac                                                                                             >  8:195970/1‑133 (MQ=255)
                                    aGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATAAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCGTCCTTCAAC‑TGGAAGAAAAACCTGGTCGACGCCACCATCGCCAAGtt                                                                                     >  3:133724/1‑133 (MQ=255)
                                                   aaGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAACAAAGACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCgg                                                                    <  4:65025/135‑1 (MQ=255)
                                                   aaGCCCGGCGCTGCGCCGAGGCGATCAAGACGGAGAACCCGCACCAGATCCTGTCGAACACCTGCGCGCCTTCCTTCAAC‑TGGAACAAAAACCTGCACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCgg                                                                    <  5:83627/135‑1 (MQ=255)
                                                   aaGCCCGGCCCTGCGCCGAGCCGATCAAGAAGGAGTCCCCGCACCAGATCCTGTCGTACAACTGCTCGCCTCCCTTCAAC‑TGCAAGAAAACCCTGCACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCgg                                                                    <  6:103145/135‑1 (MQ=255)
                                                   aaGCACGTCGCTTCGCCGTGGCGATCAAGACGGAGAACCCGGACCAGATCCTGTCGTCAAACTGCGCCCCTTCATTCAAC‑TGGAAGAAAAACCTGCACGACGCCACCATCGCCAAGGTCCAGCGCGAACTGTCgg                                                                    <  2:108158/135‑1 (MQ=255)
                                                                 gCCGAGGCTATCACGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCGCGCCTTCATTCAAC‑TGCAACAAAAACCTGCACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAg                                                      <  5:91117/135‑1 (MQ=255)
                                                                 gCCGAGGCGATCAACAAGGAGGACCCGGACCAGATCCTGGCGTACACCTGCTCGCCTCCCTTCAAC‑TGGAAGAAAACCCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAg                                                      <  8:226892/135‑1 (MQ=255)
                                                                  ccGAGGCGATCAATAAGGAGAACCCGGACAAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGCAACAAAACCCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAg                                                      <  5:133830/134‑1 (MQ=255)
                                                                  ccGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTCCCTTCAAC‑TGGAACAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAg                                                      <  8:221272/134‑1 (MQ=255)
                                                                   cGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCa                                                    <  6:167390/135‑1 (MQ=255)
                                                                        cGATCAATAAGGCGTACCCGGCCCATATCATGTCGTACAACTGCGCGCCTTCCTTCAAC‑TGCAACAAAAACCTGCACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGt                                               <  2:80419/135‑1 (MQ=255)
                                                                        cGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACATCTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGt                                               <  6:60559/135‑1 (MQ=255)
                                                                              agatgGTGTACCCTGACCACATCCTGTCGAACAACCGCTCGCGTGCCTTTATC‑TGCAACAAAAACCGGCACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTtcatc                                          <  8:348419/134‑1 (MQ=255)
                                                                                    agaaCCCGGCCCACAGCCTGTCTAACAACTGC‑CCCCATACGTAAACGAGGAATAAAACCCCGCACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTgg                                   <  2:187480/132‑1 (MQ=255)
                                                                                        cccGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCgg                               <  7:109536/135‑1 (MQ=255)
                                                                                         ccGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCCGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCCCCCTGGCCGGc                              >  8:102622/1‑135 (MQ=255)
                                                                                         ccGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGc                              >  5:43899/1‑135 (MQ=255)
                                                                                              ccAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCca                         >  4:114945/1‑135 (MQ=255)
                                                                                                aGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGAGGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCcaca                       >  8:153677/1‑135 (MQ=255)
                                                                                                     cTGTCGTACACCTGCTCGCCTTCCTTCAAC‑TGTAATAAAAACCTGGACGACGCCACAATCGACAAGATCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACAtg                  >  4:240480/1‑135 (MQ=255)
                                                                                                     cTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACAtg                  >  1:200117/1‑135 (MQ=255)
                                                                                                     cTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACAtg                  >  1:206187/1‑135 (MQ=255)
                                                                                                         cGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGc              <  5:172192/135‑1 (MQ=255)
                                                                                                          gTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCa             >  5:9525/1‑135 (MQ=255)
                                                                                                          gTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCa             >  3:117834/1‑135 (MQ=255)
                                                                                                          caacaacaGCTCGCCTTCCTACAAC‑TGCAAGAAAACCCTGCACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCa             <  3:259377/133‑1 (MQ=255)
                                                                                                                     ctccttccttCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGAGAACTGTCGGCCATGGGTTACAAGCATCAGTTAATCACCCTGGCCGGCATCCACAACATCTGGCACGGCATGTTCa  >  8:373383/3‑135 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGt                                                            >  7:18483/1‑77 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGt                                                            <  8:18483/77‑1 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCAc            >  4:316062/1‑125 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCAc            <  3:316062/125‑1 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCAc            >  4:302062/1‑125 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCAc            <  3:302062/125‑1 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACTGCATGTTCa  >  1:118692/1‑135 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCa  >  8:84087/1‑135 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCa  >  8:264357/1‑135 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCa  >  6:283989/1‑135 (MQ=255)
                                                                                                                     cGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCa  >  5:380597/1‑135 (MQ=255)
                                                                                                                          |                                                                                                                                  
CCCCGTACGCCGACCTGATCTGGTGCGAAACCGCCAAGCCGGACCTGGACGAAGCCCGTCGCTTCGCCGAGGCGATCAAGAAGGAGTACCCGGACCAGATCCTGTCGTACAACTGCTCGCCTTCCTTCAAC‑TGGAAGAAAAACCTGGACGACGCCACCATCGCCAAGTTCCAGCGCGAACTGTCGGCCATGGGTTACAAGCATCAGTTCATCACCCTGGCCGGCATCCACAACATGTGGCACGGCATGTTCA  >  ACB122/4646001‑4646252

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