Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,718,735 G→T 54.7% P124Q (CCG→CAG)  kgtP ← alpha‑ketoglutarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,718,7350GT54.7% 7.9 / 10.8 19P124Q (CCG→CAG) kgtPalpha‑ketoglutarate transporter
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base T (5/5);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACC  >  CP009273/2718594‑2718868
                                                                                                                                              |                                                                                                                                     
ccgATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAACCAATAATGCCTGAGCCCa                                                                                                                                >  2:280735/1‑149 (MQ=255)
    tCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCACTAATGCCGGAGCCCACGTa                                                                                                                            >  2:224536/1‑149 (MQ=255)
      aaCGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAACACCTTAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTAcc                                                                                                                          >  2:243287/1‑149 (MQ=255)
             caTACCGAAATGATGCGTAAAAACCTATGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCCCCGCCAACAGATAATCCCTGATATAAACGAGCGAGAAGCAATAATGCCTGAGCCCACGTACCTATAGtt                                                                                                                   <  1:259222/149‑1 (MQ=255)
                           gCGTAAAAACCTTGGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACCGATAATCCCTGAAATAAAAGAGCGAGAAGCAATAATACCGGAGCCCACGGACCTACAAGTTCATAACCCGagag                                                                                                     >  1:916/1‑149 (MQ=255)
                                aaaaCCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGaggcagg                                                                                                >  2:285024/1‑149 (MQ=255)
                                   aCCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCAGAAATAAACGAGCGAGAAGCAATAATGCCTGAGCCCTCGTACCTATAGTTTCATAACCTGGGAGGCAGGCGa                                                                                             <  1:215311/149‑1 (MQ=255)
                                          cgcCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGAATATGAGAAAGCCCTGAAATAAACGAGCGAGAAGCAATAATGCCTGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                                                                      <  2:171858/149‑1 (MQ=255)
                                          cgcCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCTGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                                                                      >  1:70146/1‑149 (MQ=255)
                                                 cAACGGCAACTTCCCTCAGATAGGTGGCGCTGGTGCCATATCCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCTGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCg                                                                               <  1:316630/149‑1 (MQ=255)
                                                    cGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGaaa                                                                            <  2:160821/149‑1 (MQ=255)
                                                                        gTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCTGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGGTaa                                                        >  1:333494/1‑149 (MQ=255)
                                                                                     cATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCTGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAACCACATCATACACACCGATAACAGCATCGAtttt                                           >  1:125067/1‑149 (MQ=255)
                                                                                        atgctcCACCAACAGATAATCCCTGAAATAACCGAGCGAGAAGCAATAATGCCTGAGCCCACGTACCTATAGTTTCATAACCTGGCAGGAAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTg                                        <  2:143976/146‑1 (MQ=255)
                                                                                                      aGATAATCCCTGAAACAAACGGGCGAGAAGCAATAATGCCGGCGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTAt                           <  1:151206/149‑1 (MQ=255)
                                                                                                              ccTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATg                   >  1:43024/1‑149 (MQ=255)
                                                                                                                 gAAATAAACGAGCGAGAAGCAATAATGCCTGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATAGGCTATGCgg                >  1:358105/1‑149 (MQ=255)
                                                                                                                    aTAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCa             <  2:32265/149‑1 (MQ=255)
                                                                                                                               agaAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAAcc  <  2:290253/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                     
CCGATCAACGTCACATACTGAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACC  >  CP009273/2718594‑2718868

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