Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,051,344 T→C D325G (GAC→GGC)  ubiI ← 2‑octaprenylphenol hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,051,3440TC95.0% 131.7 / ‑0.1 40coding (974/1203 nt)ubiI2‑octaprenylphenol hydroxylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (18/20);  minor base . (1/1);  total (19/21)
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.36e-01

CGCAGCAGTTTTTTCGCCGGATTGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGTCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAGCGCCAGACGGTG  >  NC_000913/3051215‑3051474
                                                                                                                                 |                                                                                                                                  
cGCAGCAGTTTTTTCGCCGGATTGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTg                                                                                                                          >  1:632342/1‑140 (MQ=255)
       gTTTTTTCGCCGGATTGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGc                                                                                                                   <  1:148643/140‑1 (MQ=255)
             tCGCCGGATTGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGt                                                                                                             >  2:173644/1‑140 (MQ=255)
             tCGCCGGATTGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGt                                                                                                             >  1:420610/1‑140 (MQ=255)
               gCCAGATTGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGttt                                                                                                           >  2:664213/1‑140 (MQ=255)
                 cGGATTGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCa                                                                                                         <  1:137177/140‑1 (MQ=255)
                      tGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCg                                                                                                    <  1:392188/140‑1 (MQ=255)
                              gAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAg                                                                                            <  1:644730/140‑1 (MQ=255)
                                   cAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTg                                                                                       <  1:664213/140‑1 (MQ=255)
                                       tCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTgcagc                                                                                   <  2:632342/140‑1 (MQ=255)
                                         gcgGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCAt                                                                                 >  2:484/1‑140 (MQ=255)
                                              aTCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATa                                                                            <  1:394765/140‑1 (MQ=255)
                                                  cTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGc                                                                        >  2:734367/1‑140 (MQ=255)
                                                             ccagcATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTAcc                                                             >  2:12384/1‑140 (MQ=255)
                                                               agcATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCAGCACCACCCATAAAGCCGAGATTTAcccc                                                           >  1:312686/1‑140 (MQ=255)
                                                                 cATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCCCGGCAGCATCCATAACGCCGAGATTTACCCCCt                                                         >  2:323416/1‑140 (MQ=255)
                                                                 cATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTg                                                        >  1:489402/1‑140 (MQ=255)
                                                                         ccgcACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGcc                                                 <  2:30221/140‑1 (MQ=255)
                                                                           gcACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCCTAAAGCCGGGATTTACCCCCTGCCCCGCCAg                                               >  2:308545/1‑140 (MQ=255)
                                                                           gcACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAg                                               <  2:483719/140‑1 (MQ=255)
                                                                            cACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGc                                              >  1:509685/1‑140 (MQ=255)
                                                                             aCTGTGCTTGCGGCTACGCTCATAGGGACACAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCg                                             <  1:323416/140‑1 (MQ=255)
                                                                               tgtgCTTGCGGCTACGCACATATCGACGCAGCTTACTGTTCTGCCCGAGGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCggg                                           <  2:312686/140‑1 (MQ=255)
                                                                                gtgCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGt                                          <  2:99226/140‑1 (MQ=255)
                                                                                    ttGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAAt                                      <  1:173644/140‑1 (MQ=255)
                                                                                                    aGCGACGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCg                      >  1:416935/1‑140 (MQ=255)
                                                                                                        aCGCAGATAAATGTACTGCCCGATGCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCccac                 <  2:489402/140‑1 (MQ=255)
                                                                                                         cGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCccac                 >  2:148754/1‑140 (MQ=255)
                                                                                                         cGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCccac                 <  1:108414/140‑1 (MQ=255)
                                                                                                         cGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCccac                 <  2:416935/140‑1 (MQ=255)
                                                                                                         cGCAGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCccac                 <  1:547590/140‑1 (MQ=255)
                                                                                                         cGCAGATAAATGTACTGCCCGACGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCccac                 <  1:236189/140‑1 (MQ=255)
                                                                                                            aGATAAATGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAg              >  1:382022/1‑140 (MQ=255)
                                                                                                                 aaTGTACTGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAGCGCCa         <  1:12384/140‑1 (MQ=255)
                                                                                                                        tGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGTCAGCGGGTTAATGGTAGGTGCGGCGCCGCCTACCAGCGCCAGACGGTg  >  2:241939/1‑140 (MQ=255)
                                                                                                                        tGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAGCGCCAGACGGTg  >  2:32046/1‑140 (MQ=255)
                                                                                                                        tGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAGCGCCAGACGGTg  <  2:420610/140‑1 (MQ=255)
                                                                                                                        tGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAGCGCCAGACGGTg  >  1:17658/1‑140 (MQ=255)
                                                                                                                        tGCCCGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAGCGCCAGACGGTg  <  2:700804/140‑1 (MQ=255)
                                                                                                                        ggCCCGACGCCTTCCCCCTGACGTTCCAGCCGTTTCAGTTCGAGAATCAGCCCTGCAGCATCCATAAAGCGGAGATTTCCCCCCTCCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAGCGCCAGACGGTg  <  2:382022/139‑1 (MQ=255)
                                                                                                                          cccGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCccac                 >  2:446538/1‑123 (MQ=255)
                                                                                                                          cccGATGCCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCccac                 <  1:446538/123‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
CGCAGCAGTTTTTTCGCCGGATTGGTACCGGAAAACAGATCGCGGAATCCCTGCATACCAGCCAGCATCAACGCCGCACTGTGCTTGCGGCTACGCTCATAGCGACGCAGATAAATGTACTGCCCGATGTCTTTCCCCTGACGATGCAACCGTTTCAGTTCGGCAATCAGCTCTGCAGCATCCATAAAGCCGAGATTTACCCCCTGCCCCGCCAGCGGGTGAATGGTATGTGCGGCGTCGCCCACCAGCGCCAGACGGTG  >  NC_000913/3051215‑3051474

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