Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,282,361 T→G 9.8% E290A (GAG→GCG)  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,3610TG 9.8% 116.1 / 5.5 41E290A (GAG→GCG) PP_3753AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (17/20);  new base G (2/2);  total (19/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.28e-01

AGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAATGCCC  >  NC_002947/4282223‑4282508
                                                                                                                                          |                                                                                                                                                   
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                                                                               >  5:91511/1‑136 (MQ=255)
                                                                         tCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCg                                                                               <  6:91511/136‑1 (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                                                           <  1:272178/111‑1 (MQ=255)
                                                                                                                      aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                           <  1:276256/111‑1 (MQ=255)
                                                                                                                      aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                           <  1:276257/111‑1 (MQ=255)
                                                                                                                      aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                           >  2:272178/1‑111 (MQ=255)
                                                                                                                      aCGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGa                                                           >  2:276256/1‑111 (MQ=255)
                                                                                                                          gcCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTc                 <  1:126794/149‑1 (MQ=255)
                                                                                                                           ccTTCGAGCGGCTGCGCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTc                                                                                              <  3:248463/71‑1 (MQ=38)
                                                                                                                           ccTTCGAGCGGCTGCGCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTc                                                                                              >  4:248463/1‑71 (MQ=38)
                                                                                                                            cTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCa               >  6:149134/1‑149 (MQ=255)
                                                                                                                            cTTCGAGCGGCTGCGCGAAGGCCATGGCAAGGCGCACTTCATCCTTGATCTGATCGTACGAACGCACCTACAGATACAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTcaa               >  6:92255/1‑147 (MQ=255)
                                                                                                                                     gCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAAt      <  6:25415/149‑1 (MQ=255)
                                                                                                                                         cTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAATGccc  >  3:152589/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
AGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAATGCCC  >  NC_002947/4282223‑4282508

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