Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 441,232 A→C 38.6% Q37P (CAG→CCG)  bioF → 8‑amino‑7‑oxononanoate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947441,2320AC38.6% 16.2 / 3.8 23Q37P (CAG→CCG) bioF8‑amino‑7‑oxononanoate synthase
Reads supporting (aligned to +/- strand):  ref base A (6/8);  new base C (7/2);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCC  >  NC_002947/441106‑441360
                                                                                                                              |                                                                                                                                
cgccgcGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGcg                                                                                                                         >  2:385896/1‑135 (MQ=255)
cgccgcGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCCCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCt                                                                                                                         >  2:424714/1‑136 (MQ=255)
 gccgcGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTg                                                                                                                        <  2:282279/136‑1 (MQ=255)
     cgACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCCGCGGGTGCTgggcg                                                                                                                    >  2:149779/1‑133 (MQ=255)
     cgACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCCCCCGTTGCGGGccg                                                                                                                    >  2:301413/1‑135 (MQ=255)
     cgACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGccg                                                                                                                    >  2:8384/1‑135 (MQ=255)
        cAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCCGCGGGTGCGGGCCGTCt                                                                                                                 >  2:27686/1‑136 (MQ=255)
               aCATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCaa                                                                                                          <  1:45907/136‑1 (MQ=255)
               aCATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCaa                                                                                                          <  1:27686/136‑1 (MQ=255)
                  tGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGAGGGCCAGCGGTTGCTGGCCGTCTGGAGCaaagt                                                                                                       >  2:321307/1‑133 (MQ=255)
                            ccTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGc                                                                                             >  2:327195/1‑136 (MQ=255)
                                 cAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGc                                                                                        <  2:517421/136‑1 (MQ=255)
                                       gcCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCCGCGGGTTCTGGGCGTCTGCCGCAATGACTACCAGGGGCTGGGCCaccc                                                                                  >  2:256258/1‑135 (MQ=255)
                                                          cgcAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCgccgcc                                                               <  2:411103/136‑1 (MQ=255)
                                                                  tGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCCGCCGGTGCTGGGCTTATGCCGCAATGACGAACAGGGGCTGGGCCACCACCCCGAGGTGTTAGCCCCCAGGCCGGc                                                       >  1:505430/1‑136 (MQ=255)
                                                                  tGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCCCCCGTTGCGGGGCGTGAGGGGCCATGACGACCAGGGCCAGGGCCACCCCCCCCCGGTGATCGCCCCCAGGCCGGc                                                       >  2:490609/1‑136 (MQ=255)
                                                                      tCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGgc                                                   <  1:455335/136‑1 (MQ=255)
                                                                      tCGGCAGCGGCCATTGCTGGAAAGCCCCCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGgc                                                   <  1:424714/136‑1 (MQ=255)
                                                                       cGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCCGCGGTTTCTGGGCGTGTGGAGCAATGACTACCAGGGCCTGGCCAACCCCCCCCAGGTGATCGCCGCCAGGCAGGCCCgcc                                                  >  2:311011/1‑135 (MQ=255)
                                                                                     gCTGGACACCCCGCACGGGCCGGAAGGGGGGGGCGCCGGCCCGCGGTTGCTGGCCTTCTGCACCAATGCCTCCCTGGGCCTGGCCAACCCCCCCGAGGTGATCCCCGCCGGGCAGGCCGGCGCCGCGCGCGggggg                                    <  2:396721/136‑2 (MQ=255)
                                                                                                      ggCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCCGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCGGGCAGGCCGGCGCCGCGCGGTGGGGGGTGGGGGGGGGGGcggg                   >  1:133833/1‑133 (MQ=255)
                                                                                                             aGTGGGGGGCGACGGCCCGCGGTGGGTGGCCTTCTGCACCAATGCCTCCCCGGGCCCGGCCAACCACCCCGAGGGTCTCGCCGCCTGGCCCGCCGGCGCCGCGCGCGGGGGGGGGGGGGGGGGGGCCTCGCACCTg            <  2:524640/136‑1 (MQ=255)
                                                                                                                       caCGGCCAGCGGTTGCGGGCCTTCTGCAGCAATGACTACCTGGGCCGGGCCAACCACCCCGAGGTGATCGCCGCCTGGCCGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGcc  <  2:143071/135‑1 (MQ=255)
                                                                                                                              |                                                                                                                                
CGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCC  >  NC_002947/441106‑441360

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