Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 441,601 A→T 12.2% D161V (GAC→GTC)  bioF → 8‑amino‑7‑oxononanoate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947441,6010AT12.2% 102.7 / 6.2 41D161V (GAC→GTC) bioF8‑amino‑7‑oxononanoate synthase
Reads supporting (aligned to +/- strand):  ref base A (20/16);  new base T (3/2);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

CATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGCCTGGGCACCCTGGGT  >  NC_002947/441458‑441749
                                                                                                                                                 |                                                                                                                                                    
cATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAgg                                                                                                                                                 <  3:109270/149‑1 (MQ=255)
    aCCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGt                                                                                                                                             <  2:80780/149‑1 (MQ=255)
                 ggCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCtg                                                                                                                                <  1:17374/149‑1 (MQ=255)
                      ggCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTggtggt                                                                                                                           >  6:18658/1‑149 (MQ=255)
                                     cAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCCATGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCgtgt                                                                                                             <  1:62957/149‑1 (MQ=255)
                                     cAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCgtgt                                                                                                             <  3:253697/148‑1 (MQ=255)
                                                    cATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGTGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGa                                                                                             >  4:235224/1‑149 (MQ=255)
                                                         cTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGAc                                                                                            <  6:210934/145‑1 (MQ=255)
                                                         cTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGAc                                                                                            >  5:210934/1‑145 (MQ=255)
                                                         cTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTgg                                                                                        >  5:210956/1‑149 (MQ=255)
                                                          tCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCgg                                                                                                                                          <  4:154590/98‑1 (MQ=255)
                                                          tCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCgg                                                                                                                                          >  3:154590/1‑98 (MQ=255)
                                                            gTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCTACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGcc                                                                                      <  7:140148/148‑1 (MQ=255)
                                                              tGCTGGATGGCGGGTTGCTCAGAGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCAAGCCT‑GGCCAGCCGCC‑TTGACAAGGCGGTCGGCAATACACTGCTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGAc                                                                                   >  5:215821/1‑149 (MQ=255)
                                                               gCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGAc                                                                                   <  1:261074/148‑1 (MQ=255)
                                                                  ggATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGc                                                                               >  2:111523/1‑149 (MQ=255)
                                                                        cGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGgcgc                                                                         <  5:18658/149‑1 (MQ=255)
                                                                           gtcgCTCAGTGGCGCCCGCTTCTACCGCTACCTGCACAACCACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGCACGGTGACCTGGCCGACCTGCCGGCGCTgg                                                                      <  6:253230/146‑1 (MQ=255)
                                                                                  aGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCTGGGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATg                                                                >  1:71991/1‑149 (MQ=255)
                                                                                      gtgCCCGGGTCAACCGCTACCTGCACTACGAGCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACTGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGc                                                            <  6:215821/148‑1 (MQ=255)
                                                                                      gtgCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGcc                                                           >  2:30038/1‑149 (MQ=255)
                                                                                         cccGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGt                                                        >  8:156180/1‑149 (MQ=255)
                                                                                          ccGCTTCAACTGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTg                                                       >  2:307858/1‑149 (MQ=255)
                                                                                             cTTCAACCGCTACCTGCACACCGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGccc                                                    >  5:54713/1‑149 (MQ=255)
                                                                                                  aCCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGt                                               >  8:14859/1‑149 (MQ=255)
                                                                                                          cTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCt                                       >  6:254368/1‑149 (MQ=255)
                                                                                                             cacaACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGa                                     <  6:210956/148‑1 (MQ=255)
                                                                                                                 ggcggctggCCAGGCT‑GGCCGGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTc                                >  2:259265/8‑149 (MQ=255)
                                                                                                                               cT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATgg                  >  3:208603/1‑149 (MQ=255)
                                                                                                                               at‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATgg                  <  1:258928/148‑1 (MQ=255)
                                                                                                                                t‑GGCCAGCAGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATTGACGGTGACCTGGCCGACATGCCGGCGCTGGCCAATGTCGCCCGACCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATgg                  >  2:265083/1‑148 (MQ=255)
                                                                                                                                 ttgccgaccgcc‑TTGTCCAGGCGGCTGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGc                 >  6:57015/8‑149 (MQ=255)
                                                                                                                                    ccgaccgcc‑TTGTCCAGGCGGCTGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGcc                                            >  2:144948/5‑119 (MQ=255)
                                                                                                                                    ccgaccgcc‑TTGTCCAGGCGGCTGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGcc                                            <  1:144948/115‑1 (MQ=255)
                                                                                                                                    ccgaccgcc‑TTGTCCAGGCGGCTGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGAt                                    <  7:83426/123‑1 (MQ=255)
                                                                                                                                    ccgaccgcc‑TTGTCCAGGCGGCTGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGAt                                    >  8:83426/5‑127 (MQ=255)
                                                                                                                                     cAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGCCTgg             <  7:156180/149‑1 (MQ=255)
                                                                                                                                     cAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGCCTgg             >  8:286158/1‑149 (MQ=255)
                                                                                                                                          gcc‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGCCTGGGCAcc        >  6:184363/1‑149 (MQ=255)
                                                                                                                                           cc‑TGGACAAGGTGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGCCTGGGCAccc       <  3:235224/149‑1 (MQ=255)
                                                                                                                                           cc‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGCCTGGGCAccc       >  2:51607/1‑149 (MQ=255)
                                                                                                                                                 aCAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGCCTGGGCACCCTGGGt  <  5:254368/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCT‑GGCCAGCCGCC‑TGGACAAGGCGGTCGGCAATACCCTGGTGGTGACCGATGGCGTGTTCAGCATGGACGGTGACCTGGCCGACCTGCCGGCGCTGGCCAATGTCGCCCGTGCCCGCGGTGCCTGGCTGATGGTCGACGACGCCCATGGCCTGGGCACCCTGGGT  >  NC_002947/441458‑441749

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