Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,119,525 T→G 16.0% 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,5250TG16.0% 50.7 / 6.5 25V724G (GTG→GGG) PP_3622isoquinoline 1‑oxidoreductase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (12/9);  new base G (2/2);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

TACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGC  >  NC_002947/4119386‑4119662
                                                                                                                                           |                                                                                                                                         
tACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCCGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCcggcg                                                                                                                                 >  1:48300/1‑150 (MQ=255)
                       acgcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCgcgctggc                                                                                                          >  2:165349/1‑146 (MQ=255)
                       acgcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCt                                                                                                          >  2:302828/1‑150 (MQ=255)
                         gcaATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGGGTGGCTGTATGTGCGCATCATCGAAATCGGTTCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCa                                                                                                        >  1:138666/1‑150 (MQ=255)
                                      gTACCCGATCCTGCCGCCTTCGCGGATGGCCCGTCCGTATTTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGt                                                                                           <  2:120258/150‑1 (MQ=255)
                                                    cgccTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAgcgc                                                                                             >  2:446900/1‑134 (MQ=255)
                                                    cgccTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAgcgc                                                                                             <  1:446900/134‑1 (MQ=255)
                                                     gccTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCAtt                                                                            <  2:138666/150‑1 (MQ=255)
                                                        ttCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTcgca                                                                                                               <  2:404231/112‑2 (MQ=255)
                                                        ttCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTcgca                                                                                                               >  1:404231/1‑111 (MQ=255)
                                                               aTGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCAtt                                                                            >  1:221541/1‑140 (MQ=255)
                                                               aTGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCAtt                                                                            <  2:221541/140‑1 (MQ=255)
                                                                                     tcatcGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAgg                                            <  2:48300/150‑1 (MQ=255)
                                                                                                  ggCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGAc                               <  1:290260/150‑1 (MQ=255)
                                                                                                     gAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGa                                          >  1:423477/1‑136 (MQ=255)
                                                                                                     gAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGa                                          <  2:423477/136‑1 (MQ=255)
                                                                                                     gAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCa                            >  1:236492/1‑150 (MQ=255)
                                                                                                                        gAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgct         >  2:48114/1‑150 (MQ=255)
                                                                                                                            caccGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggct     >  1:200647/1‑150 (MQ=255)
                                                                                                                             accGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggctg    >  1:162708/1‑150 (MQ=255)
                                                                                                                               cGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCgctggctggc  <  1:302828/150‑1 (MQ=255)
                                                                                                                                 gccaccgccgGGGCCACGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAAcc            <  2:178414/134‑1 (MQ=255)
                                                                                                                                 gccaccgccgGGGCCACGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAAcc            >  1:178414/5‑138 (MQ=255)
                                                                                                                                    accgccgGGGCCACGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGc                                                         <  2:187898/89‑1 (MQ=255)
                                                                                                                                    accgccgGGGCCACGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGc                                                         >  1:187898/2‑90 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TACGTGGATGGCAAGCCGCGCGCACGCAATTACGACCTGTACCCGATCCTGCCGCCTTCGCGGATGGCGCGGGTGCATGTGCGCATCATCGAAAGCGGGGCGAAGATGGGCGGCATCGGCGAACCACCGCTGCCCGCCGTGGCCCCGGCGGTGGCCAACGCGGTCGCGCAGCTCACTGGCCAGCGCGTGCGCAGTTTGCCATTGAGCCGCCACACCTTCAGCTGATCGCACAAGGGACGCCTCATGACCCATCGTTTCGCAAGAACCGCTGGCTGGC  >  NC_002947/4119386‑4119662

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