Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,654,316 A→C 22.6% N522H (AAC→CAC)  nuoG → NADH‑quinone oxidoreductase subunit G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,654,3160AC22.6% 87.1 / 2.5 50N522H (AAC→CAC) nuoGNADH‑quinone oxidoreductase subunit G
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base C (2/11);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCGCATCGCCGACGCCCTGGTCGCCGCCAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAAC  >  ACB122/4654189‑4654447
                                                                                                                               |                                                                                                                                    
ccAGCGCATCGCCGACGCCCTGGTCGCCGCCAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTc                                                                                                                              >  5:81289/1‑135 (MQ=255)
 cAGCGCATCGCCGACGCCCTGGTCGCCGCCAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCCACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCg                                                                                                                             >  1:17587/1‑135 (MQ=255)
     gcATCGCCGACGCCCTGGTCGCCGCCAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTgc                                                                                                                         >  5:361952/1‑134 (MQ=255)
     gcATCGCCGACGCCCTGGTCGCCGCAAAGCGCCCACTGGTCGTAGCCGGTACATCGCTGGCCGACCCGGCGCTGCTCGCTCCGGCTGCCCACATCGCCAACGCCCGGAAGCGGCGCGCGAAAC‑ACGGCTCGCTGa                                                                                                                         >  5:248123/1‑135 (MQ=255)
          gCCGACGCCCTGGGCGCCGCCAAGCGCCCACTGGGCGTAGCCGGTACTGCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCACGGCCCTGACGCGGCGCGAGAAAA‑ACGGCTCGCTGAGCCtg                                                                                                                    <  6:369207/135‑1 (MQ=255)
           ccGACGCCCTGGTCGCCGCCAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCtgg                                                                                                                   >  6:205135/1‑135 (MQ=255)
             gccGCCCTGGGCGCCGCCACGCGCCCACGGGTCGTCGCCGGTGCTGCGCGGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCGGCGCGAGAAAC‑ACGGCTCGCTGAGCCtggtgg                                                                                                                <  4:179534/134‑1 (MQ=255)
                gCCGTGGGCGCCGCCAAGCCCCCACTGGGCGTAGCCGGTACTGCGCGGGCCGACCCGGCGCTTATCGACGCGGCGGCCAACACCGCCACGGCCCTGAAGCGGCGCGACAAAC‑ACGGCTCGCTGAGCCTGgtggtg                                                                                                              <  1:58450/135‑1 (MQ=255)
                ccccTGGGCGCCGCCCCCCGCCCCGTGGGCGTCGCCGGTACTTCGCGGGCCGACCCGGCGCTGAGCGACGCGGCGGCCAACATCGCCAAGGCCCTGACGCGGCGCGACAAAC‑ACGGCTCGCTGAGCCTGgtggtg                                                                                                              <  2:165284/134‑1 (MQ=255)
                        cgccgcCAAGGGCCCAGGGGGCGGAGCGGTGCCTGCGCGGGCCGACGCGCGGGTGTGCTCAGCGGCGGCCACCAGCGCAACGCCCCGTACGCGCCGCGACAAAC‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGc                                                                                                      <  2:125189/135‑1 (MQ=255)
                        cgccgcCAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAC‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGcc                                                                                                     >  7:157444/1‑136 (MQ=255)
                        cgccgcCAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGc                                                                                                      >  2:100116/1‑135 (MQ=255)
                         gccgccAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGCCCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGcc                                                                                                     >  7:45431/1‑135 (MQ=255)
                             ccACGCGCCCCCTGGGCGACGCCGGTACGTCGCTGCCCCACCGGGCGATTATCGAAGCGGCGGCCAACAACGCCCCGCCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACa                                                                                                 <  6:41842/135‑1 (MQ=255)
                                 gcgcCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTc                                                                                            >  4:201091/1‑136 (MQ=255)
                                  cgcCCACGGGGCGTAGCCGGTACTTCGCTGGCCCACCCGGCGCTGATCGAAGCGGCTGCCAACAACGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCg                                                                                           <  7:251025/136‑1 (MQ=255)
                                  cgcCCACGGGCCGAACCCGGGACTTCGCGGGCCGACCGGGCGCTGAGCGAAGCGGCTGCCAACATCCCCACGGCCCTGAAGCTGCGCGAGAAAC‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCg                                                                                           <  4:288100/136‑1 (MQ=255)
                                       aCTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGACGCGGCTGCCAACATCGCCAAGGCCCTGACGCTGCGCGAGAAAC‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTg                                                                                      <  7:372357/136‑1 (MQ=255)
                                       aCTGGGCGTCGCCGGTACGTCGCGGGCCGCCCCGGCGCTGATCGAAGCGGGGGCCAACCTCGCCAAGGCCCTGAAGCTGCCCGACAAAACACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGcc                                                                                        <  2:213561/135‑1 (MQ=255)
                                            tcgtcgCCGGCACTGCGCGGGCCGACCCGGCGCTGAGCGACGCGGCTGCCAACACCCCCACGGCCCTGACGCGGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCAt                                                                                 <  7:7682/136‑1 (MQ=255)
                                            tCGTAGCGGGTAGTACGCTGGCCGACCCGGCGCTGATCGACGCGGCGGCCAACATCGCCACGGCCCTGAAGCTGCGCGACAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCAt                                                                                 <  3:11928/136‑1 (MQ=255)
                                               tAGCCGGTACTTCGCGGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCGGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCt                                                                              <  3:142748/136‑1 (MQ=255)
                                                  ccGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGGTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCg                                                                            >  2:309999/1‑135 (MQ=255)
                                                        cTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGCAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGc                                                                       >  3:117978/1‑134 (MQ=255)
                                                        cTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGc                                                                       >  7:137635/1‑134 (MQ=255)
                                                         ttCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGa                                                                     >  2:377573/1‑135 (MQ=255)
                                                                gCCGACCCGGCGCTGATCGACGCGGCGGCCAACATCGCCACGGCCCGGAAGCTCCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAgtcggtcg                                                             <  4:82160/136‑1 (MQ=255)
                                                                 ccGACCCGGCGTTGATCGAGGGGCCTGCCAACACCCCTAAGGCCCTGAAGCTGCGCTAGAAAC‑ACGGCTCGCTGAGCGTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGa                                                            <  7:158905/136‑1 (MQ=255)
                                                                 ccGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGa                                                            <  4:56084/136‑1 (MQ=255)
                                                                  cgccccGGCGCTGATCGCCGCGGCTGCCAACAGCGCCACGGCCCTGAAGCGGCGCGATAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGa                                                            <  5:132812/132‑1 (MQ=255)
                                                                  cgccccGGCGCTGATCCCCGCGGCTGCCCCCCTCGCCCCGCCCCACCAGCTGCGCGACAAAC‑ACGGCTCGCGGAGCCGGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGa                                                            <  6:260261/132‑1 (MQ=255)
                                                                   gACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATg                                                          >  4:156824/1‑136 (MQ=255)
                                                                        ggCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCAGGCTCTGTGGCGAGTCGGTCGATgccgc                                                      >  2:19524/1‑135 (MQ=255)
                                                                            cTGATCGAAGCGGCTGCCAACATCGCCACGGCCCTGAAGCGGCGCGAGAAAC‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAg                                                 <  4:98600/136‑1 (MQ=255)
                                                                             tGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAg                                                 >  5:50067/1‑135 (MQ=255)
                                                                                  gAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCACCAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCgg                                           >  7:182127/1‑136 (MQ=255)
                                                                                  caaGCGGATGCCAACATCGCCAAGTCCATGAAGCTGCACGAGAAAA‑AAGCCTCCCTGATCGTGGTCCTGCCTTAGGCCAACAGCCTCAGCCTGGCCATGCTCGGTGGCGACTCGGAAGAAGCCGCCCTAGACCCgg                                           >  8:130554/2‑136 (MQ=255)
                                                                                           gCCAACCGCCCCACGGCCCGGACGCGGCGCGAGAAAC‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCg                                  <  4:203999/136‑1 (MQ=255)
                                                                                            ccAACATCCCCCCGGCCCTGACGCGCCGCGAGAACC‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCg                                  <  6:341970/135‑1 (MQ=255)
                                                                                                cagcGCCAAGCCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCaa                              <  6:206763/132‑1 (MQ=255)
                                                                                                cATCGCCAAGGCCCTGAAGCGGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCaa                              <  2:35306/135‑1 (MQ=255)
                                                                                                        aGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATg                      >  6:173908/1‑135 (MQ=255)
                                                                                                          gCCCTGACGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGGGCCTGAGGCCAACAGCCTCGGCCGGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCa                   <  8:59821/136‑1 (MQ=255)
                                                                                                          gCCCTGAAGCGGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCa                   <  8:245901/136‑1 (MQ=255)
                                                                                                           cccTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCa                   <  6:131144/135‑1 (MQ=255)
                                                                                                                 cgCTGCGCTAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGGGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCATCgtggt            <  7:302664/135‑1 (MQ=255)
                                                                                                                        cgAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCCTCATCGGCCAAGCCGATGCCATCGTGGTGCTGGa      >  8:151274/1‑135 (MQ=255)
                                                                                                                        cgAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGa      >  4:102931/1‑135 (MQ=255)
                                                                                                                         gagaAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCCTCCGCGGCAAAGCCGATGCCATCGTGGGGCTGGa      >  1:187317/1‑134 (MQ=255)
                                                                                                                         gagaAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGGGGTGCTGGaa     >  3:283884/1‑135 (MQ=255)
                                                                                                                          ataaac‑aCGGCTCGCTGAGCCTGGTGGGGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGaaa    <  6:12500/133‑1 (MQ=255)
                                                                                                                            aaaa‑aCGGCTCGCTGAGCCTGGTGGTGCCGGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAAc  <  2:17587/135‑1 (MQ=255)
                                                                                                                               |                                                                                                                                    
CCAGCGCATCGCCGACGCCCTGGTCGCCGCCAAGCGCCCACTGGTCGTAGCCGGTACTTCGCTGGCCGACCCGGCGCTGATCGAAGCGGCTGCCAACATCGCCAAGGCCCTGAAGCTGCGCGAGAAAA‑ACGGCTCGCTGAGCCTGGTGGTGCCTGAGGCCAACAGCCTCGGCCTGGCCATGCTCGGTGGCGAGTCGGTCGATGCCGCCCTAGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAAC  >  ACB122/4654189‑4654447

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