Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,043,345 C→A 8.3% V133V (GTG→GTT)  PP_3564 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,043,3450CA 8.3% 151.2 / 4.3 48V133V (GTG→GTT) PP_3564AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base A (2/2);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

CAGTAGTGCAACGGCCAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACGTGGCGCCAGCGCCGCGGTGGCCATGGCACTGCGAC  >  NC_002947/4043204‑4043493
                                                                                                                                             |                                                                                                                                                    
cAGTAGTGCAACGGCCAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACtt                                                                                                                                               <  5:111400/149‑1 (MQ=255)
 agtagtGCAACGGCCAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACttt                                                                                                                                              >  2:52656/1‑149 (MQ=255)
        cAACGGCCAGATCTCTGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGcccc                                                                                                                                       <  6:120073/149‑1 (MQ=255)
         aaCGGCCAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCt                                                                                                                                      <  5:63129/149‑1 (MQ=255)
          aCGGCCAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTg                                                                                                                                     >  4:175263/1‑149 (MQ=255)
            ggCCAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACttt                                                                                                                                              >  5:25895/1‑138 (MQ=255)
            ggCCAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACttt                                                                                                                                              <  6:25895/138‑1 (MQ=255)
              ccAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACttt                                                                                                                                              <  6:283900/136‑1 (MQ=255)
              ccAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTgcag                                                                                                                                  >  4:52537/1‑148 (MQ=255)
              ccAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTgcag                                                                                                                                  >  6:317948/1‑148 (MQ=255)
                       tGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGcc                                                                                                                            <  3:137629/145‑1 (MQ=255)
                               gTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTgcagca                                                                                                                                <  4:113371/133‑1 (MQ=255)
                               gTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTgcagca                                                                                                                                >  3:113371/1‑133 (MQ=255)
                                                      gTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGAccc                                                                                         >  5:57889/1‑149 (MQ=255)
                                                            aTCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCccac                                                                                     >  4:280543/1‑147 (MQ=255)
                                                            aTCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCccac                                                                                     >  7:72949/1‑147 (MQ=255)
                                                             tCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAg                                                                                  >  7:152333/1‑149 (MQ=255)
                                                               gATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCt                                                                                <  3:228106/149‑1 (MQ=255)
                                                                   cAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAACCGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGac                                                                            <  8:125242/149‑1 (MQ=255)
                                                                         aaCTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCgca                                                                                                                         <  2:92347/98‑1 (MQ=255)
                                                                         aaCTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCgca                                                                                                                         >  1:92347/1‑98 (MQ=255)
                                                                         aaCTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGc                                                                      >  7:1854/1‑149 (MQ=255)
                                                                              cTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCGCACCAGCTTGACACCCGCCtt                                                                   <  6:57889/147‑1 (MQ=255)
                                                                              cTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCtt                                                                   <  3:52537/147‑1 (MQ=255)
                                                                                 gcgcAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTggcgg                                                              >  1:164494/1‑149 (MQ=255)
                                                                                    cAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGTCCCCCACCAGCTTGACACCCGCCTTGGCGgcct                                                           <  5:317948/149‑1 (MQ=255)
                                                                                    cAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACATCCGCCTTGGCGgcct                                                           <  8:72949/149‑1 (MQ=255)
                                                                                    cAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGgcct                                                           <  2:164494/149‑1 (MQ=255)
                                                                                    cAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGgcct                                                           <  1:52656/149‑1 (MQ=255)
                                                                                      aCAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGgcctgc                                                         <  6:89431/149‑1 (MQ=255)
                                                                                       cAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGgcc                                                            <  8:73136/145‑1 (MQ=255)
                                                                                       cAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACACGCCTTGGCGgcctag                                                         >  7:73136/1‑146 (MQ=255)
                                                                                           tGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAg                                                    >  5:101962/1‑149 (MQ=255)
                                                                                                 tGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCCGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGc                                              >  8:189534/1‑149 (MQ=255)
                                                                                                   cGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGc                                              <  3:175263/147‑1 (MQ=255)
                                                                                                   cGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCt                                            <  8:152333/149‑1 (MQ=255)
                                                                                                     gAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGa                                                     <  4:121176/138‑1 (MQ=255)
                                                                                                     gAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGa                                                     >  3:121176/1‑138 (MQ=255)
                                                                                                     gAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTgg                                          <  8:1854/149‑1 (MQ=255)
                                                                                                         cGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAAc                                      <  3:280543/149‑1 (MQ=255)
                                                                                                          gcagcaCGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACg                                     >  5:286317/1‑149 (MQ=255)
                                                                                                                cGAAACTTGGCTTCGTGACGCCAGTGCACAACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGcc                                             >  1:287327/1‑135 (MQ=255)
                                                                                                                cGAAACTTGGCTTCGTGACGCCAGTGCACAACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGcc                                             >  1:291971/1‑135 (MQ=255)
                                                                                                                cGAAACTTGGCTTCGTGACGCCAGTGCACAACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGcc                                             <  2:287327/135‑1 (MQ=255)
                                                                                                                cGAAACTTGGCTTCGTGACGCCAGTGCACAACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGcc                                             <  2:291971/135‑1 (MQ=255)
                                                                                                                        ggCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACGTGGCGCCAGcgccg                       >  2:17037/1‑149 (MQ=255)
                                                                                                                                   gCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCTCCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACGTGGCGCCAGCGCCGCGGTGGCCATg            >  2:226642/1‑149 (MQ=255)
                                                                                                                                        tGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACGTGGCGCCAGCGCCGCGGTGGCCATGGCACt       >  4:194949/1‑149 (MQ=255)
                                                                                                                                            ccaccacTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACGTGGCGCCAGCGCCGCGGTGGCCATgg           >  7:266332/1‑141 (MQ=255)
                                                                                                                                            ccaccacTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACGTGGCGCCAGCGCCGCGGTGGCCATgg           <  8:266332/141‑1 (MQ=255)
                                                                                                                                             caccacTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACGTGGCGCCAGCGCCGCGGTGGCCATGGCACTGCGAc  >  2:152989/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CAGTAGTGCAACGGCCAGATCTATGGCGGCGGTGCCACCGGCGCACGTGATGCGGTTGCCATCGATGCAGTACAACTGCTCGCGCAACAGCTGCAACTGCGGAAACGCAGCACGAAACTCGGCTTCGTGACGCCAGTGCACCACCACTTTGTGCCCCTGCAGCAGGCCGCAGGCAGCCAGCAGGAAGGCGCCATTGTCGACCCCCACCAGCTTGACACCCGCCTTGGCGGCCTGCCTGAGCAACGCCTGGTAACGTGGCGCCAGCGCCGCGGTGGCCATGGCACTGCGAC  >  NC_002947/4043204‑4043493

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