Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,357,816 C→T 10.3% E180K (GAA→AAA)  PP_2072 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,357,8160CT10.3% 100.2 / 3.7 39E180K (GAA→AAA) PP_2072AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (18/17);  new base T (2/2);  total (20/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 = 6.87e-01

GTCGAGGCGGCGCAGGGTGACAAAGCGATAGGGGCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTCCTGTTCTTCCAATGGTTCCAGCGCGTGGTCCACT  >  NC_002947/2357682‑2357958
                                                                                                                                      |                                                                                                                                               
gTCGAGGCGGCGCAGGGTGACAAAGCGATAGGGGCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCg                                                                                                                                  <  2:1150468/150‑1 (MQ=255)
         gcgcAGGGTGACAAAGCGATAGGGGCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTc                                                                                                                         >  2:850513/1‑150 (MQ=255)
                   aCAAAGCGATAGGGGCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTgg                                                                                                               <  1:221572/150‑1 (MQ=255)
                    cAAATCGATAGGGGCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGt                                                                                                              <  1:642269/150‑1 (MQ=255)
                          gATAGGGGCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCa                                                                                                        <  1:850513/150‑1 (MQ=255)
                               gggCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGgcgc                                                                                                   >  1:427330/1‑150 (MQ=255)
                                ggCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCggg                                                                                                                                >  2:902528/1‑120 (MQ=255)
                                ggCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCggg                                                                                                                                <  1:902528/120‑1 (MQ=255)
                                   tGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGgc                                                                                                                               >  2:1175762/1‑118 (MQ=255)
                                   tGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGgc                                                                                                                               <  1:1175762/118‑1 (MQ=255)
                                       gCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGccc                                                                                           <  1:409926/150‑1 (MQ=255)
                                               gCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCt                                                                                   >  2:394941/1‑150 (MQ=255)
                                                  tGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTTCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTgc                                                                                >  1:301571/1‑150 (MQ=255)
                                                  tGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTTCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTgc                                                                                >  1:301885/1‑150 (MQ=255)
                                                  tGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTgc                                                                                >  2:821719/1‑150 (MQ=255)
                                                     gAAAGTCCCGCGACAGGCTCCAGCATTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTgccgc                                                                             >  2:417253/1‑150 (MQ=255)
                                                          tCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTTCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑g                                                                       <  2:301571/150‑1 (MQ=255)
                                                          tCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTTCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑g                                                                       <  2:301885/150‑1 (MQ=255)
                                                            ccGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGc                                                                     <  1:355262/150‑1 (MQ=255)
                                                            ccGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTACAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGc                                                                     >  1:310559/1‑150 (MQ=255)
                                                                     gCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGc                                                            <  1:1271691/150‑1 (MQ=255)
                                                                       tCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGa                                                          >  2:326278/1‑150 (MQ=255)
                                                                          aGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCg                                                       <  1:821719/150‑1 (MQ=255)
                                                                              ttCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAg                                                   <  1:326278/150‑1 (MQ=255)
                                                                                        cTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAgg                                         <  1:394941/150‑1 (MQ=255)
                                                                                        cTGGCCTGTTACAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAgg                                         <  2:310559/150‑1 (MQ=255)
                                                                                           gCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCtt                                                                                                 >  1:747638/1‑92 (MQ=255)
                                                                                           gCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCtt                                                                                                 <  2:747638/92‑1 (MQ=255)
                                                                                             cTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTc                                    <  1:417253/150‑1 (MQ=255)
                                                                                                ttCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTCCTg                                 >  1:110514/1‑150 (MQ=255)
                                                                                                   cAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTCCTGttc                              <  2:427330/150‑1 (MQ=255)
                                                                                                          tCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTCCTGTTCTTCCAAt                       >  2:527538/1‑150 (MQ=255)
                                                                                                             aGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTCCTGTTCTTCCAATGGt                    >  1:279997/1‑150 (MQ=255)
                                                                                                                 ttATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGAGCTTGCGCACGCGCGGCTTGCCGCCGACAAGCACGC‑ACACAAGCGCCAGGTCAGACGGGGGGTCCGGTTGTTGCAATGGTTCCa                >  2:402666/1‑150 (MQ=255)
                                                                                                                       ccACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGGCCTGTTCTTCCAATGGTTCCAGCGCGt          >  2:804231/1‑150 (MQ=255)
                                                                                                                          caTGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTCCTGTTCTTCCAATGGTTCCAGCGCGTGGt       >  2:27045/1‑150 (MQ=255)
                                                                                                                             gCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAg                                          <  2:413478/112‑1 (MQ=255)
                                                                                                                             gCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAg                                          >  1:413478/1‑112 (MQ=255)
                                                                                                                               aGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTCCTGTTCTTCCAATGGTTCCAGCGCGTGGTCCACt  >  1:753937/1‑150 (MQ=255)
                                                                                                                                      |                                                                                                                                               
GTCGAGGCGGCGCAGGGTGACAAAGCGATAGGGGCTGGTGCCGTACAGCGTGCGAAAGTCCCGCGACAGGCTCCAGCGTTCGCGGCCACTGGCCTGTTCCAGCATGTCCAGTGTTATGCCCACATGCAGGTGTTCCATGATGAAGGCTCGGGCGCGTTCGGCTGCCTGGTAATCCAGGCGCTTGCGCCCGCGCGGCTTGCCGCCGAC‑GGCCCGCAAGGCGATCGCCAGGTCATACAGGGCGTCCTGTTCTTCCAATGGTTCCAGCGCGTGGTCCACT  >  NC_002947/2357682‑2357958

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