Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 546,991 G→C 10.0% Q104H (CAG→CAC fdrA → putative acyl‑CoA synthetase FdrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913546,9910GC10.0% 89.3 / 3.1 36Q104H (CAG→CACfdrAputative acyl‑CoA synthetase FdrA
Reads supporting (aligned to +/- strand):  ref base G (11/21);  new base C (2/2);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCAAACCAGGCGCTGGATCGCAACCTCAACGTGATGATGTTCTCCGAT  >  NC_000913/546849‑547126
                                                                                                                                              |                                                                                                                                       
agaaCGCCACGCCGAACGATGTTTGCGGGGCAAGTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATGATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGCCGAAg                                                                                                                                         <  2:193485/141‑1 (MQ=255)
    cgcccggccGACCGATATTTGCGTGCAAATTAGTATCGGGGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGc                                                                                                                                     <  1:203019/137‑1 (MQ=255)
    cgccacgccGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCACGTGc                                                                                                                                     >  2:102539/1‑143 (MQ=255)
        acgccGAACGAGATGTGCGTGGGAATTCGGAGCGAAGCGGCGGATGGGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCg                                                                                                                                 <  1:764445/143‑1 (MQ=255)
         ggccGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCACCTGCGTcaa                                                                                                                                <  1:604722/142‑3 (MQ=255)
            cGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTgg                                                                                                                             <  1:661282/143‑1 (MQ=255)
            cGAACGATAGTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGAGTATGCAGCAGCTTGAAGAGGCGCTAAAGCAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTgg                                                                                                                             <  2:761369/143‑1 (MQ=255)
             gAACGATATTTTCGTAGCAAAACGTAGCGAAGCCGCGGATGCGGGGATCGCCCAGCCGATTATACATCAGCTCGCAGAGGCGCTAAAACAACTGCCCCACGGGTCAAGCAGCACACCCCCCTTGACGCAGGGCCGGCGCTggc                                                                                                                            >  2:75520/1‑142 (MQ=255)
              aaCGATATTGGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCACCTGCGTcaacgcct                                                                                                                           <  2:1092247/143‑8 (MQ=255)
                      ttGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCt                                                                                                                   <  2:66990/143‑1 (MQ=255)
                      tggCGTGGCAATTGGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCt                                                                                                                   <  2:1085104/141‑1 (MQ=255)
                             gCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGaa                                                                                                            <  1:1034434/143‑1 (MQ=255)
                             gCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGaa                                                                                                            <  1:351504/143‑1 (MQ=255)
                              cAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGaaa                                                                                                           <  1:69719/143‑1 (MQ=255)
                               aaTTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCGTGACCGaaag                                                                                                          >  2:130153/1‑142 (MQ=255)
                                     gggcgcAGCGGCGGAGGCGGGGATCGCGCAGGCGACTAGGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCg                                                                                                    <  2:329353/141‑1 (MQ=255)
                                              ggcgcgtGCGGGGATCGGGCAGGCGAGTATGAAGCAGCTTGAGGAGGCGGTAAAGCAACGGCCGAAGGGGTCAGGGCGCAGGCGGGCGTTGACGCAGGTGCGGCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATc                                                                                           <  2:802628/137‑1 (MQ=255)
                                                cggATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTg                                                                                         <  1:550883/143‑1 (MQ=255)
                                                      cGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTg                                                                                   >  1:664167/1‑143 (MQ=255)
                                                       ggggATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGa                                                                                  >  1:27061/1‑143 (MQ=255)
                                                                      gATTATGCAGCAGCTTGAAGAGGCGCTAAAACACCTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTTACGTTGGTGTGTCGCTAAGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCg                                                                   <  1:819681/143‑1 (MQ=255)
                                                                       gttgtGCAGCAGCTTGAAGCGGCGCCAAAGCAAGTGGCGCGGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGa                                                                  <  2:490837/139‑1 (MQ=255)
                                                                       aTTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGa                                                                  <  1:711659/143‑1 (MQ=255)
                                                                           tGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGCCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTAt                                                              >  1:991189/1‑143 (MQ=255)
                                                                               gcagcTTGAAGAGGCGCTAAATCAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATgcgg                                                          <  2:913539/143‑1 (MQ=255)
                                                                                        ggagGGGCCAGAACTACTGGCGCGGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACATTGGCTGTCAGGAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGTCGAGTATGCGGCGGAGCTgg                                                 <  2:977074/142‑1 (MQ=255)
                                                                                          agGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCa                                               >  1:553747/1‑143 (MQ=255)
                                                                                            gcgcTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCaaa                                             >  1:453303/1‑143 (MQ=255)
                                                                                                tAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCAAACCAg                                         >  1:888960/1‑143 (MQ=255)
                                                                                                 aaaaCAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGCCTGGCGAGTATGCGGCGGAGCTGGCAAACCAgg                                        >  2:26669/1‑143 (MQ=255)
                                                                                                        ctgccccAGGGGTCAGGCAGACGCCAGGCGTTGACGCAGGTGCGTCGCGGGGACAGTGCCTGTCAGAAATTACCCCATGCAACTCTGGCTTTGATTTCAGTGGCTGGCGAGGATGCGGCGGAGCTGGCAAACCAGGCGCTGGa                                 <  1:304646/137‑1 (MQ=255)
                                                                                                          ggCGCAGGGGTCAGGCAGCAGCCAGGCGCTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCAAACCAGGCGCTGGATc                               <  1:81419/143‑1 (MQ=255)
                                                                                                              cAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTAGGCGGCGGAGCTGGCAAACCAGGCGCTGGACCGCaa                           >  2:185726/1‑143 (MQ=255)
                                                                                                                         cagcagCCAGGCGTTGACGCAGGGGCGTCGCTGGGACAGTGCCTGTCAGCAATTACCCGATGCCCATATGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGATGGCACACACGGCGAAGGAGCGCAACCTCAACGtga                >  2:9312/1‑143 (MQ=255)
                                                                                                                          agcagcCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCAAACCAGGCGCTGGATCGCAACCTCAACGtgat               <  1:662573/143‑1 (MQ=255)
                                                                                                                                      ttGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCAAACCAGGCGCTGGATCGCAACCTCAACGTGATGATGTTCTCCGa   <  1:353652/143‑1 (MQ=255)
                                                                                                                                       cgACGCACCTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATGTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCAAACCAGGCGCTGGATCGCAACCTCAACGTGATGATGTTCTCCGAt  >  2:943986/2‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                       
ATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGGGTCAGGCAGCAGCCAGGCGTTGACGCAGGTGCGTCGCTGGGACAGTGCCTGTCAGAAATTACCCGATGCCAATCTGGCGCTGATTTCAGTGGCTGGCGAGTATGCGGCGGAGCTGGCAAACCAGGCGCTGGATCGCAACCTCAACGTGATGATGTTCTCCGAT  >  NC_000913/546849‑547126

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