Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,825,091 G→A 8.7% G235D (GGC→GAC)  PP_2478 → isoquinoline 1‑oxidoreductase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,825,0910GA 8.7% 125.6 / 2.3 46G235D (GGC→GAC) PP_2478isoquinoline 1‑oxidoreductase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base A (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

GAGGCGGCTGCGACGCTGCCATTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGAGCAACTTCTGGAC  >  NC_002947/2824943‑2825238
                                                                                                                                                     |                                                                                                                                                   
gAGGCGGCTGCGACGCTGCCATTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCgg                                                                                                                                                     <  2:318031/149‑1 (MQ=255)
 aGGCGGCTGCGACGCTGCCATTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGc                                                                                                                                                    <  3:160955/149‑1 (MQ=255)
   gcggcTGCGACGCTGCCATTGCCCGACGCAGTGCAGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGACAt                                                                                                                                                  >  6:173788/1‑149 (MQ=255)
         gCGACGCTGCCATTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCtggt                                                                                                                                            <  7:345929/149‑1 (MQ=255)
          cGACGCTGCCATTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCcggtg                                                                                                                                           <  7:71163/149‑1 (MQ=255)
           gACGCTGCCATTGCCCGCCGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGCAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCtggtgg                                                                                                                                          >  3:333618/1‑149 (MQ=255)
           gACGCTGCCATTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCtggtgg                                                                                                                                          >  1:236901/1‑149 (MQ=255)
                   cATTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAg                                                                                                                                   >  4:161079/1‑148 (MQ=255)
                    aTTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTa                                                                                                                                 >  7:311558/1‑149 (MQ=255)
                     ttGCCCGACGCAGTGCCGCTGAAAACCACCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCTACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTAt                                                                                                                                <  3:19820/149‑1 (MQ=255)
                       gCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCtggt                                                                                                                                            >  8:209062/1‑135 (MQ=255)
                       gCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCtggt                                                                                                                                            <  7:209062/135‑1 (MQ=255)
                       gCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGACATGCTGGTGGCCTCCAGTATc                                                                                                                               <  5:173788/148‑1 (MQ=255)
                        cccGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCAc                                                                                                                             >  4:395312/1‑149 (MQ=255)
                             cGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTc                                                                                                                        <  5:334778/149‑1 (MQ=255)
                                    ccGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCTCACTTTACCATCGACCTGCAAATCCCCgg                                                                                                                                                     <  1:334903/113‑1 (MQ=255)
                                    ccGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCgg                                                                                                                                                     >  2:340517/1‑113 (MQ=255)
                                    ccGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCgg                                                                                                                                                     <  1:340517/113‑1 (MQ=255)
                                    ccGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGAACTGCAAAACCCCgg                                                                                                                                                     >  2:334903/1‑113 (MQ=255)
                                     cGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTAt                                                                                                                <  3:161079/149‑1 (MQ=255)
                                      gCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATg                                                                                                               >  6:345582/1‑149 (MQ=255)
                                        tGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGTACTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATgg                                                                                                              >  1:96024/1‑149 (MQ=255)
                                            aaCCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATggcgggc                                                                                                         <  1:409460/149‑1 (MQ=255)
                                            aaCCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATggcgggc                                                                                                         <  6:223303/149‑1 (MQ=255)
                                                      cAGTTCAAGTTGCTGGGCACCCCCGCGCAGCCT‑CTGGATACGCCAGCGAAGGTCAACGGCAAGGCGCACTTTACCCTCGACCTGCAAATCCCCGGCATGAGGGAGGCCTCCAGTAACACCAGACCGGAGTATGGCGCGCGGCTCCGCaa                                                                                               >  3:115322/1‑149 (MQ=255)
                                                        gTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACg                                                                                             >  7:193883/1‑149 (MQ=255)
                                                           cAAGTTGCTGGGCACAGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGt                                                                                            <  1:165277/147‑1 (MQ=255)
                                                            aaGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTtg                                                                                          <  1:104605/148‑1 (MQ=255)
                                                            aaGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTCCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATAGCGGGCGGCTGCTCAACGTtga                                                                                         <  8:4270/149‑1 (MQ=255)
                                                             aGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTtgat                                                                                        <  2:96024/149‑1 (MQ=255)
                                                               ttGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGCCCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTtgatga                                                                                      <  2:236901/149‑1 (MQ=255)
                                                                 gCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACg                                                                                             >  5:183782/1‑140 (MQ=255)
                                                                 gCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACg                                                                                             <  6:183782/140‑1 (MQ=255)
                                                                   tGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGt                                                                                  >  8:352549/1‑149 (MQ=255)
                                                                   tGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGACATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGt                                                                                  >  1:211450/1‑149 (MQ=255)
                                                                     ggCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTGTACCATCGCCCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGAGGGCGGCTGCGCCACGGTGATGACAGt                                                                                  >  3:319172/1‑147 (MQ=255)
                                                                      gCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAg                                                                               >  7:107824/1‑149 (MQ=255)
                                                                          cgccgcGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGgcgc                                                                           >  4:142478/1‑149 (MQ=255)
                                                                               cgccgcGC‑CTGGCTACGCCTGCGCATCGCCACGGCACTGCGCACTTTCCCATCGACCTGCACCTCCCCGGCATGCGGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCGCAACGTTGATGACAGTCAGGCGcgacga                                                                      <  4:319172/149‑1 (MQ=255)
                                                                                                gCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTgcgc                                                       >  7:201043/1‑148 (MQ=255)
                                                                                                        aGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGACATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTg                                              <  2:211450/149‑1 (MQ=255)
                                                                                                          gTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTgcg                                            >  3:363772/1‑149 (MQ=255)
                                                                                                                     ggCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATg                                 >  3:428506/1‑149 (MQ=255)
                                                                                                                                       ccTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGa               <  3:395312/149‑1 (MQ=255)
                                                                                                                                       ccTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGa               <  8:311558/149‑1 (MQ=255)
                                                                                                                                        cTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATc                                                 <  2:429328/114‑1 (MQ=255)
                                                                                                                                        cTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATc                                                 >  1:429328/1‑114 (MQ=255)
                                                                                                                                        cTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGAg              <  3:25207/149‑1 (MQ=255)
                                                                                                                                            aaaTCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGAGCAAc          >  5:190806/1‑149 (MQ=255)
                                                                                                                                                 cccGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGAGCAACTTCTg     >  7:243084/1‑149 (MQ=255)
                                                                                                                                                  ccGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGAGCAACTTCTg     >  5:434096/1‑148 (MQ=255)
                                                                                                                                                    ggCATGCTGGTGGCCTCCAGTATCACCTGTCCGGGGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGAGCAACTTCTGGAc  <  3:298714/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAGGCGGCTGCGACGCTGCCATTGCCCGACGCAGTGCCGCTGAAAACCCCCGAGCAGTTCAAGTTGCTGGGCACCGCCGCGCAGCGT‑CTGGATACGCCTGCGAAGGTCAACGGCAAGGCGCACTTTACCATCGACCTGCAAATCCCCGGCATGCTGGTGGCCTCCAGTATCACCTGTCCGGTGTATGGCGGGCGGCTGCTCAACGTTGATGACAGTCAGGCGCGACGAGTGCCTGGGGTGCGCGACATCGTGCGCCTGGACAATGCCGTGGCCGTCACCGCGAGCAACTTCTGGAC  >  NC_002947/2824943‑2825238

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