Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,119,486 C→A 7.7% P726T (CCG→ACG)  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,4860CA 7.7% 161.1 / 4.0 52P726T (CCG→ACG) PP_3622isoquinoline 1‑oxidoreductase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (20/28);  new base A (2/2);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

TGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGCC  >  NC_002947/4119346‑4119628
                                                                                                                                             |                                                                                                                                              
tGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTg                                                                                                                                        <  3:188110/149‑1 (MQ=255)
 ggATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTgg                                                                                                                                       >  7:223195/1‑149 (MQ=255)
  gATGGCAAGCCGCGCGCTCGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGc                                                                                                                                      <  7:20513/149‑1 (MQ=255)
       cAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAAcg                                                                                                                                 >  1:96781/1‑149 (MQ=255)
          gCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCgg                                                                                                                              >  6:118194/1‑149 (MQ=255)
          gCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGGGGCCAcggcgggca                                                                                                                                       <  1:166912/140‑4 (MQ=255)
          gCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGGGGCCAcggcgggca                                                                                                                                       >  2:166912/1‑137 (MQ=255)
           ccGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGt                                                                                                                             >  7:5936/1‑149 (MQ=255)
            cgcgcgCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTc                                                                                                                            <  8:5936/149‑1 (MQ=255)
                       aaTTACGACCTGTACCCGATCCTGGCGCCTTCGCGGCTGGCGCGGGTGCATGTGCGCAG‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCAc                                                                                                                 <  6:216735/149‑1 (MQ=255)
                       aaTTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCAc                                                                                                                 >  7:9343/1‑149 (MQ=255)
                         ttACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTg                                                                                                               <  3:203473/149‑1 (MQ=255)
                         ttACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTg                                                                                                               <  3:286759/149‑1 (MQ=255)
                             gACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGGGGCCAcggcggg                                                                                                                                         <  6:33992/119‑2 (MQ=255)
                             gACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGGGGCCAcggcggg                                                                                                                                         >  5:33992/1‑118 (MQ=255)
                               ccTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTgg                                                                                                                                       <  7:278605/119‑1 (MQ=255)
                               ccTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTgg                                                                                                                                       >  8:278605/1‑119 (MQ=255)
                                   tACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCa                                                                                                                  >  2:67416/1‑136 (MQ=255)
                                   tACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCACTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCa                                                                                                                  <  1:67416/136‑1 (MQ=255)
                                     cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTcgc                                                                                                                          <  3:6662/126‑1 (MQ=255)
                                     cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTcgc                                                                                                                          >  4:6662/1‑126 (MQ=255)
                                     cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgc                                                                                                   <  6:275929/149‑1 (MQ=255)
                                     cccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTg                                                                                                    <  8:175322/148‑1 (MQ=255)
                                      ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                      <  8:238454/129‑1 (MQ=255)
                                      ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                      >  7:238454/1‑129 (MQ=255)
                                      ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                      >  7:239349/1‑129 (MQ=255)
                                      ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAg                                                                                                                      <  8:239349/129‑1 (MQ=255)
                                      ccGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgcg                                                                                                  <  1:9590/149‑1 (MQ=255)
                                        gATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTgcg                                                                                                  >  5:145619/1‑147 (MQ=255)
                                             tGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTg                                                                                           >  3:121012/1‑149 (MQ=255)
                                                     tCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAg                                                                                   >  7:15901/1‑149 (MQ=255)
                                                      cgcgGATGGCGCGGGTGCATGTGCGCATGCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACt                                                                                                                >  1:175217/1‑120 (MQ=255)
                                                         ggATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAgccgc                                                                               <  6:69986/149‑1 (MQ=255)
                                                              gcgcGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCacac                                                                          <  3:178726/149‑1 (MQ=255)
                                                                gcgGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGcc                                                                                         >  4:25013/1‑132 (MQ=255)
                                                                gcgGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGcc                                                                                         <  3:25013/132‑1 (MQ=255)
                                                                 cgGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCtt                                                                       <  6:145619/149‑1 (MQ=255)
                                                                  gggTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTc                                                                      <  4:121012/149‑1 (MQ=255)
                                                                                at‑catcGAAAGCGGGGCGAAGCTGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCACCGCGTGCGCAGGTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACaa                                                        <  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:9343/149‑1 (MQ=255)
                                                                                                          gcggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTc                               <  8:224753/149‑1 (MQ=255)
                                                                                                          gcggcATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTc                               <  7:270652/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)
                                                                                                                                             |                                                                                                                                              
TGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCAT‑CATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGCTGGCGCTGCC  >  NC_002947/4119346‑4119628

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