Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,119,481 T→G 6.8% V724G (GTG→GGG)  PP_3622 → isoquinoline 1‑oxidoreductase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,119,4810TG 6.8% 169.1 / 2.4 56V724G (GTG→GGG) PP_3622isoquinoline 1‑oxidoreductase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (22/30);  new base G (2/2);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

AAGCGGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGCC  >  NC_002947/4119334‑4119628
                                                                                                                                                     |                                                                                                                                                   
aagCGGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCTCATCATCGAAAGCCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGt                                                                                                                                                     <  1:218571/149‑1 (MQ=255)
 agCGGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTTCCCGCCGTgg                                                                                                                                                   <  4:134671/149‑1 (MQ=255)
  gCGGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTgg                                                                                                                                                   >  8:119223/1‑148 (MQ=255)
  gCGGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTgg                                                                                                                                                   <  7:119223/148‑1 (MQ=255)
  gCGGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGc                                                                                                                                                  >  4:188110/1‑149 (MQ=255)
            tGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTg                                                                                                                                        <  3:188110/149‑1 (MQ=255)
             ggATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTgg                                                                                                                                       >  7:223195/1‑149 (MQ=255)
              gATGGCAAGCCGCGCGCTCGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGc                                                                                                                                      <  7:20513/149‑1 (MQ=255)
                   cAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAAcg                                                                                                                                 >  1:96781/1‑149 (MQ=255)
                      gCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCgg                                                                                                                              >  6:118194/1‑149 (MQ=255)
                      gCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGGGGCCAcggcgggca                                                                                                                                       <  1:166912/140‑4 (MQ=255)
                      gCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGGGGCCAcggcgggca                                                                                                                                       >  2:166912/1‑137 (MQ=255)
                       ccGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGt                                                                                                                             >  7:5936/1‑149 (MQ=255)
                        cgcgcgCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTc                                                                                                                            <  8:5936/149‑1 (MQ=255)
                                   aaTTACGACCTGTACCCGATCCTGGCGCCTTCGCGGCTGGCGCGGGTGCATGTGCGCAGCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCAc                                                                                                                 <  6:216735/149‑1 (MQ=255)
                                   aaTTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCAc                                                                                                                 >  7:9343/1‑149 (MQ=255)
                                     ttACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTg                                                                                                               <  3:286759/149‑1 (MQ=255)
                                     ttACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTg                                                                                                               <  3:203473/149‑1 (MQ=255)
                                         gACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGGGGCCAcggcggg                                                                                                                                         <  6:33992/119‑2 (MQ=255)
                                         gACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGGGGCCAcggcggg                                                                                                                                         >  5:33992/1‑118 (MQ=255)
                                           ccTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTgg                                                                                                                                       >  8:278605/1‑119 (MQ=255)
                                           ccTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTgg                                                                                                                                       <  7:278605/119‑1 (MQ=255)
                                               tACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCa                                                                                                                  >  2:67416/1‑136 (MQ=255)
                                               tACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCACTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCa                                                                                                                  <  1:67416/136‑1 (MQ=255)
                                                 cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTcgc                                                                                                                          >  4:6662/1‑126 (MQ=255)
                                                 cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTcgc                                                                                                                          <  3:6662/126‑1 (MQ=255)
                                                 cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgc                                                                                                   <  6:275929/149‑1 (MQ=255)
                                                 cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTg                                                                                                    <  8:175322/148‑1 (MQ=255)
                                                  ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                      >  7:239349/1‑129 (MQ=255)
                                                  ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                      >  7:238454/1‑129 (MQ=255)
                                                  ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                      <  8:238454/129‑1 (MQ=255)
                                                  ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                      <  8:239349/129‑1 (MQ=255)
                                                  ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgcg                                                                                                  <  1:9590/149‑1 (MQ=255)
                                                    gATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgcg                                                                                                  >  5:145619/1‑147 (MQ=255)
                                                         tGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTg                                                                                           >  3:121012/1‑149 (MQ=255)
                                                                 tCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAg                                                                                   >  7:15901/1‑149 (MQ=255)
                                                                  cgcgGATGGCGCGGGTGCATGTGCGCATGCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACt                                                                                                                >  1:175217/1‑120 (MQ=255)
                                                                     ggATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAgccgc                                                                               <  6:69986/149‑1 (MQ=255)
                                                                          gcgcGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCacac                                                                          <  3:178726/149‑1 (MQ=255)
                                                                            gcgGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGcc                                                                                         >  4:25013/1‑132 (MQ=255)
                                                                            gcgGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGcc                                                                                         <  3:25013/132‑1 (MQ=255)
                                                                             cgGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCtt                                                                       <  6:145619/149‑1 (MQ=255)
                                                                              gggTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTc                                                                      <  4:121012/149‑1 (MQ=255)
                                                                                            atcatcGAAAGCGGGGCGAAGCTGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCACCGCGTGCGCAGGTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACaa                                                        <  2:96781/149‑1 (MQ=255)
                                                                                                 tcGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCtt                                                                       <  1:37363/130‑1 (MQ=255)
                                                                                                 tcGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCtt                                                                       >  2:37363/1‑130 (MQ=255)
                                                                                                  cGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGAc                                                   >  7:224753/1‑149 (MQ=255)
                                                                                                            ggCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGa                                          <  8:15901/149‑1 (MQ=255)
                                                                                                               gAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGAcc                                        >  1:68543/1‑148 (MQ=255)
                                                                                                                 aGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCa                                      >  4:10089/1‑148 (MQ=255)
                                                                                                                       gcggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTc                               <  8:224753/149‑1 (MQ=255)
                                                                                                                       gcggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTc                               <  7:270652/149‑1 (MQ=255)
                                                                                                                       gcggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTc                               <  8:9343/149‑1 (MQ=255)
                                                                                                                         ggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGc                             <  8:223195/149‑1 (MQ=255)
                                                                                                                                  gAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgc                    >  2:211801/1‑149 (MQ=255)
                                                                                                                                    accaccGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctg                  <  7:158097/149‑1 (MQ=255)
                                                                                                                                                    gTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGcc  <  5:118194/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAGCGGTGTACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGCC  >  NC_002947/4119334‑4119628

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: