Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,724,990 C→A 8.5% L748L (CTG→CTT)  sucA ← 2‑oxoglutarate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,724,9900CA 8.5% 180.1 / 4.0 56L748L (CTG→CTT) sucA2‑oxoglutarate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (16/35);  new base A (3/2);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GTGGCGCAGCAGCGACTTCGGCGTCAGGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGAT  >  ACB143/4724858‑4725121
                                                                                                                                    |                                                                                                                                   
gtgGCGCAGCAGCGACTTCGGCGTCAGGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGcaca                                                                                                                                    <  4:292180/134‑1 (MQ=255)
    cgcAGCAGCGACTTCGGCGTCAGGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCa                                                                                                                              <  4:78115/136‑1 (MQ=255)
          agcGACTTCGGCGTCAGGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAAc                                                                                                                        >  1:42800/1‑136 (MQ=255)
          agcGACTTCGGCGTCAGGACTACCAGCGGCGTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAAc                                                                                                                        >  2:66071/1‑136 (MQ=255)
                   ggCGTCAGGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGAc                                                                                                               <  6:29858/136‑1 (MQ=255)
                     cGTCAGGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGAcgc                                                                                                             <  2:42800/136‑1 (MQ=255)
                           taCTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCggag                                                                                                       <  7:11962/135‑1 (MQ=255)
                                 cAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGc                                                                                                 <  1:191091/136‑1 (MQ=255)
                                      gCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGTAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCgg                                                                                            <  7:157009/136‑1 (MQ=255)
                                         tGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGccc                                                                                         >  2:32862/1‑136 (MQ=255)
                                          gcgcAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGccc                                                                                         <  7:201796/135‑1 (MQ=255)
                                          gcgcAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGccc                                                                                         <  7:236552/135‑1 (MQ=255)
                                          gcgcAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGccc                                                                                         <  3:138127/135‑1 (MQ=255)
                                          gcgcAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGccc                                                                                         <  2:148474/135‑1 (MQ=255)
                                                   gCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGCCCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCa                                                                               >  6:160260/1‑136 (MQ=255)
                                                    cGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCAt                                                                              <  1:170791/136‑1 (MQ=255)
                                                       aTGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGc                                                                           <  4:177550/136‑1 (MQ=255)
                                                        gcgccTGACGACGCAGCATGTGGTAAACCTGTGCCGGGGTAGACGGTACGCAAACCTGCATGTTTTGCTCAGCACAAAGTTGCAGATAACGTTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCGTAAcc                                                                          >  4:131126/4‑136 (MQ=255)
                                                               gCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTgg                                                                    >  1:215003/1‑135 (MQ=255)
                                                                 gACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcag                                                                 <  5:277574/136‑1 (MQ=255)
                                                                 gACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcag                                                                 <  7:221186/136‑1 (MQ=255)
                                                                 gACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcag                                                                 <  3:18996/136‑1 (MQ=255)
                                                                  aCGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcag                                                                 >  1:212059/1‑135 (MQ=255)
                                                                   cGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCa                                                                                                                  <  3:44157/85‑1 (MQ=255)
                                                                   cGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCa                                                                                                                  >  4:44157/1‑85 (MQ=255)
                                                                      agcagGTGGTAAATCTTTGACTGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcagcagca                                                            >  1:203109/1‑136 (MQ=255)
                                                                      agcagGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGgcagcagca                                                            >  7:268341/1‑136 (MQ=255)
                                                                         agGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATgg                                                         <  2:295576/136‑1 (MQ=255)
                                                                         agGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATgg                                                         <  6:139598/136‑1 (MQ=255)
                                                                           gtggtAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTc                                                       <  4:1286/136‑1 (MQ=255)
                                                                                aaaTCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGAc                                                                                                               >  2:145780/1‑75 (MQ=255)
                                                                                aaaTCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGAc                                                                                                               <  1:145780/75‑1 (MQ=255)
                                                                                 aaccTGTGCCGGGGTAGACGGTACGCAAACCTGCATGTTTTGCTCAGCACAAAGTTGCAGATAACGTTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCGTAACCGTGCGGCAGCAACATCACCAGAcca                                                 >  8:225118/4‑135 (MQ=255)
                                                                                        gCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggc                                          <  8:189464/136‑1 (MQ=255)
                                                                                         ccGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcg                                         >  2:81654/1‑136 (MQ=255)
                                                                                         ccGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggc                                          <  5:75663/135‑1 (MQ=255)
                                                                                          cGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcgg                                        >  3:317907/1‑136 (MQ=255)
                                                                                          cGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggcgg                                        >  2:236017/1‑136 (MQ=255)
                                                                                                    ggcaCGCAGACCTGGAGCTTGTGCGCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCTGCCCCTGCCCTTCATAGCCGTGTGTCATCATCTTGGTCAGACCGCACAGGCGGCCCCACTtgt                              <  4:158777/133‑1 (MQ=255)
                                                                                                    ggTACGCAGACCTGGCTGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTACAGACGCTCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCTTGTTGAATCAGTATTTTAATAACTCACAGGCGGCCCCACTtgt                              <  8:78381/136‑1 (MQ=255)
                                                                                                    ggTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTtgt                              >  7:37155/1‑136 (MQ=255)
                                                                                                    ggTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTtgt                              <  6:147029/136‑1 (MQ=255)
                                                                                                          cAGACCTGGATGTTGTGGTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTcac                        <  1:309556/136‑1 (MQ=255)
                                                                                                          cAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTcac                        <  1:150643/136‑1 (MQ=255)
                                                                                                          cAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTcac                        <  5:217922/136‑1 (MQ=255)
                                                                                                          cAGACCTGGATGTTGTGCTCGGCGCAAAGCTGCAGGTAACGCTCCAGACGCTCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTcac                        <  8:69071/136‑1 (MQ=255)
                                                                                                           aaaCCTGCATGTTTTGCTCAGCACAAAGTTGCAGATAACGATCCAGACGCGCGGAGGAGTGCTCCGGCCCATGCCCTTCGTAACCGTGCGGCAGCAACATCACCAGACCACACATCCGGCCCCATTTCTGTTcgc                        >  5:28990/3‑133 (MQ=255)
                                                                                                                   aTGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGtga               <  8:223158/136‑1 (MQ=255)
                                                                                                                      ttttgCTCAGCACAAAGTTGCAGATAACGTTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCGTAACCGTGCGGCAGCAACATCACCAGACCACACATCCGGCCCCATTTCTGTTCGCCAGAGGagatg             <  7:225118/132‑1 (MQ=255)
                                                                                                                         tgCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACt         <  1:17972/136‑1 (MQ=255)
                                                                                                                          gCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTg        >  1:16269/1‑136 (MQ=255)
                                                                                                                            tCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggc                                          <  5:81609/100‑1 (MQ=255)
                                                                                                                            tCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAggc                                          >  6:81609/1‑100 (MQ=255)
                                                                                                                              ggCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCg    <  3:142349/136‑1 (MQ=255)
                                                                                                                              ggCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCg    <  6:85711/136‑1 (MQ=255)
                                                                                                                                cgcACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGAt  <  3:163731/136‑1 (MQ=255)
                                                                                                                                cgcACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGAt  <  5:9536/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GTGGCGCAGCAGCGACTTCGGCGTCAGGACTACCAGCGGCTTGCGCAGCGGGCGAATGACCTGGCGACGCAGCAGGTGGTAAATCTGTGCCGGCGTCGTCGGTACGCAGACCTGGATGTTGTGCTCGGCGCACAGCTGCAGGTAACGCTCCAGACGCGCGGAGGAGTGCTCCGGCCCCTGCCCTTCATAGCCGTGTGGCAGCAGCATGGTCAGACCGCACAGGCGGCCCCACTTGTGCTCACCACTGGTGATGAACTGGTCGAT  >  ACB143/4724858‑4725121

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