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% 135.8 / 4.9 50R291Q (CGG→CAG) PP_3753AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base T (3/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

ACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAATG  >  NC_002947/4282210‑4282505
                                                                                                                                                     |                                                                                                                                                   
aCAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGcc                                                                                                                                                     >  2:14982/1‑149 (MQ=255)
  aGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGc                                                                                                                                                   >  7:268709/1‑149 (MQ=255)
   ggTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCt                                                                                                                                                  >  7:99722/1‑149 (MQ=255)
       gAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCg                                                                                                                                                >  1:151382/1‑147 (MQ=255)
       gAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGaa                                                                                                                                              <  1:148473/149‑1 (MQ=255)
       gAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGAGCGGCTGcgcg                                                                                                                                                >  7:217083/1‑147 (MQ=255)
          cGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGAGCGGCTGCGCGAAGGc                                                                                                                                           <  8:217083/149‑1 (MQ=255)
            cAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGAGCGGCTGCGCGAAGGcgc                                                                                                                                         <  7:134392/149‑3 (MQ=255)
             aGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCAt                                                                                                                                        <  2:151382/149‑1 (MQ=255)
                 cGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGcc                                                                                                                                    >  2:2484/1‑149 (MQ=255)
                  gACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCAt                                                                                                                                        <  8:68185/144‑1 (MQ=255)
                  gACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCAt                                                                                                                                        >  7:68185/1‑144 (MQ=255)
                                ggTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATcc                                                                                                                     <  8:268709/149‑1 (MQ=255)
                                ggTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCGCGAAGGCCATGGCCCGGCGCACTTCATcc                                                                                                                     >  7:369124/1‑149 (MQ=255)
                                 gTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCt                                                                                                                    <  2:79578/149‑1 (MQ=255)
                                        ccTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATTCTCAAAGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGAGCAATTCATCCTTCATct                                                                                                              >  4:260547/1‑148 (MQ=255)
                                             gCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGt                                                                                                        <  8:369124/149‑1 (MQ=255)
                                                   cccTGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCTGTTCCTTAAACCCCGCCTGCTCGGCATTTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACg                                                                                                     <  2:341571/146‑1 (MQ=255)
                                                     cGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGc                                                                                                <  6:161991/149‑1 (MQ=255)
                                                         ccACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATc                                                                                                                      <  2:16165/123‑1 (MQ=255)
                                                         ccACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATc                                                                                                                      >  1:16165/1‑123 (MQ=255)
                                                         ccACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGcc                                                                                               >  6:259772/1‑146 (MQ=255)
                                                         ccACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGcc                                                                                               <  5:259772/146‑1 (MQ=255)
                                                          cACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTc                                                                                           <  8:99722/149‑1 (MQ=255)
                                                           aCTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTcc                                                                                          <  1:14982/149‑1 (MQ=255)
                                                                       cccGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGc                                                                              >  7:24689/1‑149 (MQ=255)
                                                                                aGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAg                                                                     >  5:84073/1‑149 (MQ=255)
                                                                                aGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAg                                                                     <  4:214145/149‑1 (MQ=255)
                                                                                     cTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACt                                                                >  2:71957/1‑149 (MQ=255)
                                                                                          ttCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCAc                                                           <  8:24689/149‑1 (MQ=255)
                                                                                             cTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGAt                                                                                                                >  1:44984/1‑93 (MQ=255)
                                                                                             cTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGAt                                                                                                                <  2:44984/93‑1 (MQ=255)
                                                                                              tGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCg                                                                         <  3:372417/131‑1 (MQ=255)
                                                                                              tGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCg                                                                         >  4:372417/1‑131 (MQ=255)
                                                                                              tGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCAc                                                           <  5:137144/145‑1 (MQ=255)
                                                                                              tGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCAc                                                           >  6:137144/1‑145 (MQ=255)
                                                                                              tGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGAt                                                       >  7:6080/1‑149 (MQ=255)
                                                                                               gAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGAt                                                       >  1:96076/1‑149 (MQ=255)
                                                                                                 aaCCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATccc                                                    <  3:327149/149‑1 (MQ=255)
                                                                                                  aCCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTg                                                         <  6:126135/143‑1 (MQ=255)
                                                                                                  aCCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTg                                                         >  5:126135/1‑143 (MQ=255)
                                                                                                        ccTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTtgc                                             <  1:2484/149‑1 (MQ=255)
                                                                                                            cTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTtgctgc                                          <  8:34079/148‑1 (MQ=255)
                                                                                                               ggCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTgcgc                                        >  3:168693/1‑147 (MQ=255)
                                                                                                                        ggCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCa                             <  1:71957/149‑1 (MQ=255)
                                                                                                                                        ccTTCGAGCGGCTGCGCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTc              >  7:238878/1‑148 (MQ=255)
                                                                                                                                         cTTCGAGCGGCTGCGCGAAGGCCATGGCCCGGCGCACTTCATGCTTGATCTGCTGGTACCAACGCCCCTCCCGCTCCAACTTCCGGCGAAAGGTACAGGCACTGATCCCCTGTTGCTGAGCCAGTGCCGCCAATGTCCGCCACTGTcc             >  4:373247/1‑148 (MQ=255)
                                                                                                                                            cgcgCAGCCGCTCGCAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGTTACTGGCACTGATCCGCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATa          <  3:260547/148‑1 (MQ=255)
                                                                                                                                                gcAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCGAACTGGCGGCGAAAGGTACTGGCACTGAt                                                       <  5:339706/100‑1 (MQ=255)
                                                                                                                                                gcAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAgcg       >  4:211008/1‑148 (MQ=255)
                                                                                                                                                gcAGCCGCTCGAAGGCCATCGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTCCCGGCGAAAGGTACTCGCACTGAt                                                       >  6:339706/1‑100 (MQ=255)
                                                                                                                                                    ccGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAATg  >  6:263713/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACAGGTAGAACGCAGGCCGACCAGCCGCGCCCGGTAGCTCCCTGGGCTCTGCCCGGTCCACTTCTTGAACGCCCGGTGGAAGGCACTCGGTTCCTGAAACCCCGCCTGCTCGGCAATTTCGGCAATGCTCAACGCGCCTTCGCGCAGCCGCTCGAAGGCCATGGCCCGGCGCACTTCATCCTTGATCTGCTGGTACGAACGCCCCTCCCGCTCCAACTGCCGGCGAAAGGTACTGGCACTGATCCCCTGTTGCTGCGCCAGTGCGGCCAATGTCGGCCACTGTCCATAGCGCAATG  >  NC_002947/4282210‑4282505

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: