Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,389,315 T→A 16.1% L926Q (CTG→CAG)  dld2 → D‑lactate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,389,3150TA16.1% 69.6 / 8.8 31L926Q (CTG→CAG) dld2D‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (15/11);  new base A (2/3);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01

TTTACCACCCCCGAGCTGAACGCCCACTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGCCACCACCGCATTCGCGG  >  NC_002947/5389170‑5389461
                                                                                                                                                 |                                                                                                                                                  
tttACCACCCCCGAGCTGAACGCCCACTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTa                                                                                                                                                <  2:462960/150‑1 (MQ=255)
      accCCCGAGCTGAACGCCCACTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACcgc                                                                                                                                          <  1:526420/150‑1 (MQ=255)
        ccccGAGCTGAACGCCCACTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCAGGTAGAccaggc                                                                                                                                        <  2:595088/150‑5 (MQ=255)
                       ccACTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCATAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTGGACCGCGTCACCCGGCCACGCat                                                                                                                         >  1:380165/1‑149 (MQ=255)
                        cACTCGTTGCGCAGCCTGAAGGTCGCCGTGCAGTATTGCAGCGAAGTTAGTGCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGCTATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGc                                                                                                                        <  1:78252/150‑1 (MQ=255)
                        cACTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCAGGTAGAccagg                                                                                                                                         <  1:258087/133‑4 (MQ=255)
                        cACTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCAGGTAGAccagg                                                                                                                                         >  2:258087/1‑130 (MQ=255)
                        cACTCGTTGCGCAGCCTCAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACAACGGCATGGTCTACCTGTGAGACCGCGTCACGCGGCCACGCAGc                                                                                                                        >  2:314796/1‑150 (MQ=255)
                          cTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCct                                                                                                                      >  1:197146/1‑150 (MQ=255)
                                  gcAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTgg                                                                                                                                                  <  2:97008/114‑1 (MQ=255)
                                  gcAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTgg                                                                                                                                                  >  1:97008/1‑114 (MQ=255)
                                             gACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCAGGTAGAcca                                                                                                                                           <  2:597762/110‑2 (MQ=255)
                                             gACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCAGGTAGAcca                                                                                                                                           >  1:597762/1‑109 (MQ=255)
                                              aCGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACCCAGCCTCTAACGCGTTGATATAGGTg                                                                                                  >  2:98538/1‑150 (MQ=255)
                                                gCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTggg                                                                                                <  2:125571/150‑1 (MQ=255)
                                                 ccGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTgggg                                                                                               >  2:536475/1‑150 (MQ=255)
                                                      cAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGt                                                                                                             >  1:371990/1‑131 (MQ=255)
                                                      cAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGt                                                                                                             <  2:371990/131‑1 (MQ=255)
                                                                         tttCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTAc                                                                       >  1:397435/1‑150 (MQ=255)
                                                                                                  aGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGc                                              >  2:383234/1‑150 (MQ=255)
                                                                                                                    cGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTc                            <  1:124988/150‑1 (MQ=255)
                                                                                                                           gACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCa                     >  1:385789/1‑150 (MQ=255)
                                                                                                                                 cACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGccacc               >  1:10300/1‑150 (MQ=255)
                                                                                                                                 cACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGccacc               >  1:316943/1‑150 (MQ=255)
                                                                                                                                  aCGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGccacca              <  2:10300/150‑1 (MQ=255)
                                                                                                                                      ccTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGCCACCACCGc          >  2:532848/1‑150 (MQ=255)
                                                                                                                                           tCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGCCACCACCGCATTcg     >  1:582340/1‑150 (MQ=255)
                                                                                                                                             tACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGCCACCACCGCATTcgcg   <  1:430126/150‑1 (MQ=255)
                                                                                                                                             tACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGCCACCACCGCATTcgcg   <  1:98538/150‑1 (MQ=255)
                                                                                                                                              aCCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGCCACCACCGCATTCGCgg  >  1:449850/1‑150 (MQ=255)
                                                                                                                                              aCCTGGTAAACCGCGTCACCCGGCCGCGCAGCCTCTAACGCATTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGCCACCACCGCATTCGCgg  <  2:397435/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTTACCACCCCCGAGCTGAACGCCCACTCGTTGCGCAGCCTGAAGGACGCCGTGCAGTATTGCAGCGAAGGTATTTCCACCAGCCGCACCTGCGAGATAGGCCTGTCGAGCCACAGCGGCATCGACTACCACGGCCTGGTCTACCTGGTAGACCGCGTCACCCGGCCACGCAGCCTCTAACGCGTTGATATAGGTGGGGAGCGCTGCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCCTGCTGGCCAATCCCCACCCACTGTCCACCGCAGCCACCACCGCATTCGCGG  >  NC_002947/5389170‑5389461

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