Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,282,358 C→T 10.0% R291Q (CGG→CAG)  PP_3753 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,282,3580CT10.0% 108.0 / 2.7 40R291Q (CGG→CAG) PP_3753AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base T (2/2);  total (18/22)
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

ACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAAT  >  NC_002947/4282210‑4282504
                                                                                                                                                    |                                                                                                                                                  
aCAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGAGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGcc                                                                                                                                                    >  4:179498/1‑149 (MQ=255)
             aGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGAGCGGCTGCGCGAAGGcgcg                                                                                                                                       <  5:226926/149‑4 (MQ=255)
              ggCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATg                                                                                                                                      <  6:107616/149‑1 (MQ=255)
                            gcCCGGTAGCTCCCTGGGGTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGTCCGGCGCACtt                                                                                                                         >  5:278918/1‑148 (MQ=255)
                              ccGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCAt                                                                                                                      <  8:268303/149‑1 (MQ=255)
                              ccGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCAt                                                                                                                      >  5:209861/1‑149 (MQ=255)
                                         cTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATctgc                                                                                                           <  1:176474/149‑1 (MQ=255)
                                         cTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATctgc                                                                                                           <  3:303198/149‑1 (MQ=255)
                                            ggCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATctgctg                                                                                                         <  6:209861/148‑1 (MQ=255)
                                            agCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATctgct                                                                                                          <  3:179498/146‑1 (MQ=255)
                                                              tcttGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGc                                                                                               <  3:250635/140‑1 (MQ=255)
                                                              tcttGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGc                                                                                               >  4:250635/1‑140 (MQ=255)
                                                               cttGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGc                                                                                      <  8:267378/148‑1 (MQ=255)
                                                                 aataCGGCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTcc                                                                                   <  1:189178/146‑1 (MQ=255)
                                                                  gAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCt                                                                                     <  2:322348/146‑1 (MQ=255)
                                                                           gTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATc                                                                                                                     <  6:289788/105‑1 (MQ=255)
                                                                           gTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATc                                                                                                                     >  5:289788/1‑105 (MQ=255)
                                                                               aaGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGTCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGaaa                                                                     <  6:278918/149‑1 (MQ=255)
                                                                                    aCTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTAc                                                                >  2:126794/1‑149 (MQ=255)
                                                                                      tCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCg                                                                           <  6:91511/136‑1 (MQ=255)
                                                                                      tCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCg                                                                           >  5:91511/1‑136 (MQ=255)
                                                                                                   ccccGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCt                                                 >  6:181368/1‑149 (MQ=255)
                                                                                                         cTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTtcct                                           >  4:31812/1‑149 (MQ=255)
                                                                                                               ggCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTgcgc                                       >  7:292203/1‑147 (MQ=255)
                                                                                                               ggCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGcc                                      <  5:181368/148‑1 (MQ=255)
                                                                                                               ggCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGGAcc                                                               >  8:196358/1‑122 (MQ=255)
                                                                                                                    tttCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGGTCCCCTGTTGCTGCGCCAGTGCg                                >  4:280084/1‑149 (MQ=255)
                                                                                                                                tCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGc                    >  5:25415/1‑149 (MQ=255)
                                                                                                                                  aaCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCa                  <  8:292203/149‑1 (MQ=255)
                                                                                                                                   aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                       >  2:276257/1‑111 (MQ=255)
                                                                                                                                   aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                       >  2:276256/1‑111 (MQ=255)
                                                                                                                                   aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                       >  2:272178/1‑111 (MQ=255)
                                                                                                                                   aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                       <  1:276257/111‑1 (MQ=255)
                                                                                                                                   aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                       <  1:276256/111‑1 (MQ=255)
                                                                                                                                   aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                       <  1:272178/111‑1 (MQ=255)
                                                                                                                                       gcCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTc             <  1:126794/149‑1 (MQ=255)
                                                                                                                                        ccTTCGAGCGGCTGCGCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTc                                                                                          >  4:248463/1‑71 (MQ=38)
                                                                                                                                        ccTTCGAGCGGCTGCGCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTc                                                                                          <  3:248463/71‑1 (MQ=38)
                                                                                                                                         cTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCa           >  6:149134/1‑149 (MQ=255)
                                                                                                                                         cTTCGAGCGGCTGCGCGAAGGCCATGGCAAGGCGCACTTCATCCTTGATCTGATCGTACGAACGCACCTACAGATACAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTcaa           >  6:92255/1‑147 (MQ=255)
                                                                                                                                                  gCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAAt  <  6:25415/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAAT  >  NC_002947/4282210‑4282504

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