Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 680,138 G→A 9.7% C342Y (TGC→TAC)  PP_0582 → thiolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947680,1380GA 9.7% 111.0 / 2.4 41C342Y (TGC→TAC) PP_0582thiolase family protein
Reads supporting (aligned to +/- strand):  ref base G (20/17);  new base A (2/2);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

ACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGAT  >  NC_002947/679990‑680285
                                                                                                                                                    |                                                                                                                                                   
agcGCTGCGGGGGAGTTTGTCTCGGGGCGGGATGGGTTGTTGATGGCGGCGGTGTATGCCGTTCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCtgtg                                                                                                                                                     <  8:258887/147‑1 (MQ=255)
   gCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCAc                                                                                                                                                  <  6:172928/149‑1 (MQ=255)
                ttGTCACGGGGCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGTACGCTCAAGGCCTg                                                                                                                                      <  2:251132/148‑1 (MQ=255)
                 tGTCACGGGGCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTgg                                                                                                                                     >  4:207477/1‑148 (MQ=255)
                    cACGGGGCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTggg                                                                                                                                    >  5:130884/1‑146 (MQ=255)
                      cGGGGCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGAt                                                                                                                               <  3:207477/149‑1 (MQ=255)
                       ggggCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAg                                                                                                                                           >  3:242328/1‑136 (MQ=255)
                       ggggCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAg                                                                                                                                           <  4:242328/136‑1 (MQ=255)
                         ggCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGcc                                                                                                                            >  4:54509/1‑149 (MQ=255)
                           cgcgAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGa                                                                                                                          <  6:64836/149‑1 (MQ=255)
                            gcgAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGACACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGc                                                                                                                                                    >  2:28973/1‑122 (MQ=255)
                            gcgAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGACACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGc                                                                                                                                                    <  1:28973/122‑1 (MQ=255)
                                          aTGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACg                                                                                                           <  6:237814/149‑1 (MQ=255)
                                                 cGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTACACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGc                                                                                                    >  4:243639/1‑149 (MQ=255)
                                                        tgccgtgccgCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAg                                                                                              <  6:130884/148‑1 (MQ=255)
                                                              gccgCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGcgccgc                                                                                       >  8:195866/1‑149 (MQ=255)
                                                                       gCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGgcgc                                                                                <  3:54509/147‑1 (MQ=255)
                                                                           gCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGa                                                                          >  2:161121/1‑149 (MQ=255)
                                                                               gCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAg                                                                                                                                           >  6:193825/1‑80 (MQ=255)
                                                                               gCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAg                                                                                                                                           <  5:193825/80‑1 (MQ=255)
                                                                               gCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTACACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGc                                                                                                    >  8:303059/1‑119 (MQ=255)
                                                                               gCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTACACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGc                                                                                                    <  7:303059/119‑1 (MQ=255)
                                                                                 aaTGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCgcag                                                                    >  3:259754/1‑149 (MQ=255)
                                                                                     gCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCaaa                                                                <  7:195866/149‑1 (MQ=255)
                                                                                           cGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTACACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAAt                                                          <  3:243639/149‑1 (MQ=255)
                                                                                            gCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATg                                                         <  1:161121/149‑1 (MQ=255)
                                                                                             cTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGt                                                        <  4:259754/149‑1 (MQ=255)
                                                                                                      tttGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAg                                               >  1:7502/1‑149 (MQ=255)
                                                                                                            tactacGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATTCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGAt                                                                             >  5:315635/1‑113 (MQ=255)
                                                                                                            tactacGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGAt                                                                             <  6:315635/113‑1 (MQ=255)
                                                                                                            tactacGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCt                                         >  6:22733/1‑149 (MQ=255)
                                                                                                             actacGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCaa                                                           >  5:296795/1‑130 (MQ=255)
                                                                                                             actacGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCaa                                                           <  6:296795/130‑1 (MQ=255)
                                                                                                                        cACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCa                             >  2:190365/1‑149 (MQ=255)
                                                                                                                         aCGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCAc                            >  7:134600/1‑149 (MQ=255)
                                                                                                                              gCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGtt                       <  1:190365/149‑1 (MQ=255)
                                                                                                                                cTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTg                     <  4:331958/149‑1 (MQ=255)
                                                                                                                                cTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCACTGGCGGCCGGGCACCCGTTTg                     >  1:40390/1‑149 (MQ=255)
                                                                                                                                      ggCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGAGCGACCGCCGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCaa               >  4:324710/1‑149 (MQ=255)
                                                                                                                                       gCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAAc              >  6:197948/1‑149 (MQ=255)
                                                                                                                                         ccAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCg            >  4:179602/1‑149 (MQ=255)
                                                                                                                                             gTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTgg        >  6:365880/1‑149 (MQ=255)
                                                                                                                                                   tgCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGat  <  5:22733/149‑1 (MQ=255)
                                                                                                                                                   tgCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGa   <  2:7502/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCGCTGCGGTGGATTTTGTCACGGGGCGCGAAGGGTTGTTGATGGCGCCGGTGTATGCCGTGCCGCGCTTGCTGGCACGCAATGGCCTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGAT  >  NC_002947/679990‑680285

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: ‑