Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,544,403 T→C T148A (ACC→GCC)  bioH ← pimeloyl‑acyl carrier protein methyl ester esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,544,4030TC95.5% 64.4 / ‑4.4 22T148A (ACC→GCC) bioHpimeloyl‑acyl carrier protein methyl ester esterase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (11/10);  total (12/10)
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.96e-01

GACGGAGATCGACCGTTTTCAGGATTTCCAGCCCGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACGAGCACTAAAACAAGGTGACGACGCCACGGTGACCAG  >  NC_000913/3544277‑3544541
                                                                                                                              |                                                                                                                                          
gACGGAGATGGGCCGTTTTCAGGATTTCCAGCCCGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAgg                                                                                                                               <  1:44782/140‑1 (MQ=255)
     agaTCGCCCGTTTTCAGGATTTCCAGCCCGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCg                                                                                                                          <  1:642325/140‑1 (MQ=255)
     agaCCGACCGTTTTCAGGATTTCCAGCCCGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCg                                                                                                                          <  2:655196/140‑1 (MQ=255)
          gACCTTTTTCAGGATTTCCAGCCCGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCa                                                                                                                     <  1:25841/140‑1 (MQ=255)
               ttttCAGGATTTCCAGCCCGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGACCCCATGGCTTGTAACGCCAGTAACCGCTCCACTGtc                                                                                                                >  1:186064/1‑139 (MQ=255)
                          tCCAGCCCGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATTCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAatc                                                                                                     >  2:471871/1‑140 (MQ=255)
                                ccGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCGCCGCTGTACGCTGAAAATCATCACt                                                                                               >  1:234879/1‑140 (MQ=255)
                                  gCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGa                                                                                             >  1:488453/1‑140 (MQ=255)
                                           cACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTgctgct                                                                                    >  1:17849/1‑140 (MQ=255)
                                              gTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGaa                                                                                 >  1:540545/1‑140 (MQ=255)
                                               tCAACCTCCGGCATCTGTAACGCCACAACTGTTTTCTCCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGaaa                                                                                <  1:188104/140‑1 (MQ=255)
                                                   ccTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATTATTACTGAGTTGCTGCTGAAATccc                                                                            >  1:211518/1‑140 (MQ=255)
                                                        ggCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAg                                                                       >  1:18560/1‑140 (MQ=255)
                                                        ggCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAg                                                                       >  1:399455/1‑140 (MQ=255)
                                                         gCATCGGTAACGCCAGAGCGGTCTTCTTCACTGCCCGCGTGTCCTGGCACGCCGTTTCAGCCCCCATGGCTTGAAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTCCTGAAATCCCGCCAGc                                                                      <  2:186064/140‑1 (MQ=255)
                                                             cGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGt                                                                  >  1:429569/1‑140 (MQ=255)
                                                                                        cGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCGGCCAGCAAGTCCGGTTTTATCCCCGGCCACTCGTCAc                                       >  1:288193/1‑140 (MQ=255)
                                                                                         gCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACg                                      <  2:18560/140‑1 (MQ=255)
                                                                                         gCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACg                                      <  2:36871/140‑1 (MQ=255)
                                                                                             gcgcATCCTGGCGCGCCGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACGAGCa                                  <  2:488453/140‑1 (MQ=255)
                                                                                                  tCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCGCCGGCCACTCGTCACGAGCACTaaa                             >  1:415723/1‑140 (MQ=255)
                                                                                                           gcCGTTTCAGTCCCCATGGCTTTTAACGCCAGGAATCGCTCCACTGTATGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACGAGCACTAAAACAAGGTga                    <  2:234879/140‑1 (MQ=255)
                                                                                                             cGTTTCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACGAGCACTAAAACAAGGTgacg                  >  2:57050/1‑140 (MQ=255)
                                                                                                                ttCAGTCCCCATGGCTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACGAGCACTAAAACAAGGTGacgacg               <  1:471871/140‑1 (MQ=255)
                                                                                                                             gTTTGTAACGCCAGCAACCGCCCCACTGTACGCTGAAAACCACCACGGAGTGGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACGAGCACTAAAACAAGGTGACGACGCCACGGTGACCAg  <  2:415723/140‑1 (MQ=255)
                                                                                                                              |                                                                                                                                          
GACGGAGATCGACCGTTTTCAGGATTTCCAGCCCGCCATTAAGCACGTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACGAGCACTAAAACAAGGTGACGACGCCACGGTGACCAG  >  NC_000913/3544277‑3544541

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: