Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,911 C→A 14.8% G274C (GGT→TGT)  icd ← NADP(+)‑specific isocitrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,525,9110CA14.8% 82.8 / 2.1 27G274C (GGT→TGT) icdNADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (13/10);  new base A (2/2);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

TGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGC  >  NC_002947/4525780‑4526041
                                                                                                                                   |                                                                                                                                  
tGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGtg                                                                                                                                >  1:424823/1‑135 (MQ=255)
                      tAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATcc                                                                                                          >  2:274039/1‑136 (MQ=255)
                           cccGTGGAGGTTGAGGGGGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAggg                                                                                                     <  2:464547/135‑1 (MQ=255)
                             cGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGgcc                                                                                                   >  1:521922/1‑136 (MQ=255)
                              gttgaggttgagGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGgccg                                                                                                  >  1:488718/1‑136 (MQ=255)
                                          ggTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGt                                                                                      >  2:363947/1‑136 (MQ=255)
                                          ggTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGt                                                                                      >  2:357913/1‑136 (MQ=255)
                                                         gtACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCg                                                                       >  1:325793/1‑136 (MQ=255)
                                                                      cGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAg                                                                                                       <  2:441620/91‑1 (MQ=255)
                                                                      cGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAg                                                                                                       >  1:441620/1‑91 (MQ=255)
                                                                           ctggAGCTGCTGCAGCATGGCGTGGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTAc                                                     <  1:194837/134‑1 (MQ=255)
                                                                            aGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTAcc                                                    >  2:378660/1‑136 (MQ=255)
                                                                              gATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACGTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGACCTCGTCCCGCGCCACCTCGTAcccc                                                  >  2:486067/1‑136 (MQ=255)
                                                                                  tgctgcAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACAGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTAAGCGCCGAACTCGTCACGCGCCAcc                                                           <  2:202863/123‑1 (MQ=255)
                                                                                  tgctgcAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACAGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTAAGCGCCGAACTCGTCACGCGCCAcc                                                           <  2:348755/123‑1 (MQ=255)
                                                                                  tgctgcAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACAGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTAAGCGCCGAACTCGTCACGCGCCAcc                                                           >  1:348755/1‑123 (MQ=255)
                                                                                  tgctgcAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACAGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTAAGCGCCGAACTCGTCACGCGCCAcc                                                           >  1:202863/1‑123 (MQ=255)
                                                                                   gctgcAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTc                                             >  1:11262/1‑136 (MQ=255)
                                                                                              gCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCAc                                  <  2:398143/136‑1 (MQ=255)
                                                                                               cGTCGGCGATGGCGTCCTTGACGATGACCTCGCTACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCAcc                                 >  1:90266/1‑136 (MQ=255)
                                                                                                gTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCACAGTCCTTGAAGGCAccc                                <  1:236061/136‑1 (MQ=255)
                                                                                                   ggCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCg                             <  1:391258/136‑1 (MQ=255)
                                                                                                    gCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACATCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCgg                            <  1:307681/136‑1 (MQ=255)
                                                                                                    gCGATGGCGTCCTTGACGATGACCTCGCGACCGGTATAGGGGATCTTGAACTTCATCCAAGGGCCGCCATCAACCAGTTCAGCTCCGAACTCGTCACGCGCAACATCGTACCCACACACCTTGAACGCACCCTCgg                            >  1:375172/1‑136 (MQ=255)
                                                                                                                     gATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTAGGTGAACTTCATGATGTTg           <  2:424823/136‑1 (MQ=255)
                                                                                                                       tGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGcc         <  2:324008/136‑1 (MQ=255)
                                                                                                                              gcgACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGc  <  2:184301/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                  
TGCGGCCAGGGCATCGGACAGGTAGTCACCGTTGAGGTTGAGGGTGGCGATCACGTCGTACTCGGCCGGGCGTAGCAGGATCTGCTGCAGCATGGCGTCGGCGATGGCGTCCTTGACGATGACCTCGCGACCGGTCTTGGGGTTCTTGAACTTCATCCAAGGGCCGCCATCAAGCAGTTCAGCGCCGAACTCGTCACGCGCCACCTCGTACCCCCAGTCCTTGAAGGCACCCTCGGTGAACTTCATGATGTTGCCCTTGTGC  >  NC_002947/4525780‑4526041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: