Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,135,562 G→A 10.0% A411V (GCG→GTG)  fadE ← medium‑long chain acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,135,5620GA10.0% 101.4 / 3.6 40A411V (GCG→GTG) fadEmedium‑long chain acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (21/15);  new base A (2/2);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

GGTGGTACTTGAGGATCGCCGACAGCACCGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGTTCATCAACATCATCCAGCCCTTGCCGA  >  NC_002947/2135422‑2135709
                                                                                                                                            |                                                                                                                                                   
ggtggtACTTGAGGATCGCCGACAGCACCGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCa                                                                                                                                            <  2:666523/150‑1 (MQ=255)
  tggtACTTGAGGATCGCCGACAGCACCGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCACCAATGGCACg                                                                                                                                          >  1:361768/1‑150 (MQ=255)
    gtACTTGAGGATCGCCGACAGCACCGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGCAGGCCGCCAATGGCACGtt                                                                                                                                        <  1:436268/150‑1 (MQ=255)
        ttGAGGATCGCCGACAGCACCGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAAc                                                                                                                                    >  2:326939/1‑150 (MQ=255)
            ggATCGCCGACAGCACCGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCTGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCaa                                                                                                                                                 <  2:552859/133‑1 (MQ=255)
            ggATCGCCGACAGCACCGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCaa                                                                                                                                                 >  1:552859/1‑133 (MQ=255)
                            cGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAAGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTa                                                                                                                >  1:530970/1‑150 (MQ=255)
                             gAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACCGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCAGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTAc                                                                                                               >  2:480261/1‑150 (MQ=255)
                                 gCTTTTCGCCCATGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGc                                                                                                           <  2:354683/150‑1 (MQ=255)
                                   ttttCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCg                                                                                                         >  2:637836/1‑150 (MQ=255)
                                   ttttCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCACCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCg                                                                                                         <  2:361768/150‑1 (MQ=255)
                                   ttttCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCAGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCg                                                                                                         >  2:637758/1‑150 (MQ=255)
                                         cccAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCAcc                                                                                                   >  2:141694/1‑150 (MQ=255)
                                          ccAGGTCCACGGCCTTGGCGGTCCGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCa                                                                                                  <  2:488450/150‑1 (MQ=255)
                                          ccAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCa                                                                                                  >  2:146204/1‑150 (MQ=255)
                                                 cACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCACCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGgt                                                                                           >  2:430925/1‑150 (MQ=255)
                                                           gCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTggcg                                                                                 >  1:435686/1‑150 (MQ=255)
                                                                   cagGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGt                                                                         >  1:475365/1‑150 (MQ=255)
                                                                        gggCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGa                                                                    <  1:601107/150‑1 (MQ=255)
                                                                              tGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGGCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGAccgccg                                                              >  1:646023/1‑150 (MQ=255)
                                                                              tGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGAccgccg                                                              >  1:251042/1‑150 (MQ=255)
                                                                                   aTCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGt                                                         <  1:146204/150‑1 (MQ=255)
                                                                                     cAGCCAGGCGTTGCCGCAAATGCGTTCCAGCGATTCCTGTATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTgc                                                       <  1:401931/150‑2 (MQ=255)
                                                                                            gCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAAGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAg                                                <  2:530970/150‑1 (MQ=255)
                                                                                                  ccgccAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGAGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGccc                                          >  1:12428/1‑150 (MQ=255)
                                                                                                           cGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCACCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAgg                                 <  1:430925/150‑1 (MQ=255)
                                                                                                                cAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGAccgcc                                                               >  2:264673/1‑115 (MQ=255)
                                                                                                                cAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGAccgcc                                                               <  1:264673/115‑1 (MQ=255)
                                                                                                                cAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCAGCCCGCCGGAAGTGAGATGGAGCGGCCCACCGACAGGCAGtt                            >  2:528869/1‑150 (MQ=255)
                                                                                                                      ttCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGTTCATCaa                      <  1:141694/150‑1 (MQ=255)
                                                                                                                         cTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTg                                                                                             >  1:78216/1‑76 (MQ=255)
                                                                                                                         cTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTg                                                                                             <  2:78216/76‑1 (MQ=255)
                                                                                                                          tGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGTTCATCAAcatc                  >  2:392795/1‑150 (MQ=255)
                                                                                                                                 cTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGt                             <  2:254499/132‑1 (MQ=255)
                                                                                                                                 cTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGt                             >  1:254499/1‑132 (MQ=255)
                                                                                                                                 cTTCGAAGGCCGCCAATGGCACGTTGAACAGATAGCGGATAGTTGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGCGCGGCCCACCGACAGGCAGTTCGTCAACATCATCCAGc           <  2:435686/150‑1 (MQ=255)
                                                                                                                                  ttCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAg                                                <  1:504080/112‑1 (MQ=255)
                                                                                                                                  ttCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGTTCATCAACATCATCCAGcc          >  2:701505/1‑150 (MQ=255)
                                                                                                                                  ttCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCACTGAGATGGAg                                                >  2:504080/1‑112 (MQ=255)
                                                                                                                                     gAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGTTCATCAACATCATCCAGCCCtt       >  1:306762/1‑150 (MQ=255)
                                                                                                                                          ccgccAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGACACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGTTCATCAACATCATCCAGCCCTTGCCGa  >  2:244461/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GGTGGTACTTGAGGATCGCCGACAGCACCGAGGGCTTTTCGCCCAGGTCCACGGCCTTGGCGGTCAGCAGGCGGGCACTGTCCATCAGCCAGGCGTTGCCGCCAATGCGTGCCAGCGATTCCTGGATGCCTTCGAAGGCCGCCAATGGCACGTTGAACTGCTCGCGGATATTGGCGTACTGGCCGGTCACCAGGCTGGTGTACTTGGCGGCACCGGTGCCGACCGCCGGCAGTGAGATGGAGCGGCCCACCGACAGGCAGTTCATCAACATCATCCAGCCCTTGCCGA  >  NC_002947/2135422‑2135709

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