Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,097,638 C→G 10.0% G28R (GGG→CGG)  PP_3605 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,097,6380CG10.0% 115.0 / 4.9 40G28R (GGG→CGG) PP_3605LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (18/18);  new base G (2/2);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

TCGAGAACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCT  >  NC_002947/4097493‑4097784
                                                                                                                                                 |                                                                                                                                                  
tCGAGAACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGGACTCGGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCcggcg                                                                                                                                                >  2:1051300/1‑150 (MQ=255)
   agaACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGgca                                                                                                                                             <  1:898914/150‑1 (MQ=255)
    gaACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCGGGCGgcaa                                                                                                                                            <  2:1040815/150‑2 (MQ=255)
              tACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGcc                                                                                                                                  >  2:177411/1‑150 (MQ=255)
                          gCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACaa                                                                                                                      <  1:1052430/150‑1 (MQ=255)
                            ccTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAgg                                                                                                                    >  1:403297/1‑150 (MQ=255)
                                    gTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAgg                                                                                                            >  2:1082924/1‑150 (MQ=255)
                                      aaaGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGt                                                                                                          <  2:454217/150‑1 (MQ=255)
                                       aaGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACTAGCCCTGCTCATACAAGGCCAGCAGGGTg                                                                                                         <  2:89561/150‑1 (MQ=255)
                                       aaGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTg                                                                                                         <  2:333061/150‑1 (MQ=255)
                                                          tGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCAAAAGCCTCGTAGGGTGGCGCGCCATGTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGcgac                                                                                      >  1:136128/1‑150 (MQ=255)
                                                              cTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACt                                                                                  <  1:696438/150‑1 (MQ=255)
                                                              cTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACt                                                                                  >  2:696438/1‑150 (MQ=255)
                                                               tCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTg                                                                                 <  1:729342/150‑1 (MQ=255)
                                                                  accaGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCGGGCGgcaactg                                                                                                                                         <  1:718654/91‑5 (MQ=255)
                                                                  accaGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCGGGCGgcaactg                                                                                                                                         >  2:718654/1‑87 (MQ=255)
                                                                  accaGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCg                                                                              >  1:251533/1‑150 (MQ=255)
                                                                   ccaGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCCAAAATTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGa                                                                             <  1:1081875/150‑1 (MQ=255)
                                                                     aGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAg                                                                           <  1:177411/150‑1 (MQ=255)
                                                                     aGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAg                                                                           <  2:309998/150‑1 (MQ=255)
                                                                      ggCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAgg                                                                          >  1:146282/1‑150 (MQ=255)
                                                                       gCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGt                                                                         >  1:819299/1‑150 (MQ=255)
                                                                           aTGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATcc                                                                     <  2:403297/150‑1 (MQ=255)
                                                                               ccAGCCGCGTCTCTAAAGCACGCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCTAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATgg                                                                 >  2:572240/1‑150 (MQ=255)
                                                                                    cgcgTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAgg                                                            <  2:146282/150‑1 (MQ=255)
                                                                                     gcgTCTCTTAAGCACTCAGGCGTTTCGACAGGACCGTAGGGTGGNACAANAGTTGCAGCCCGGCGGCAGCAAGCAAGCCCTGCTCATACAAGGCAAGAAGGTTGGCAAGTTCATCAGCGCGACGACTGTCGCAGGTATCCGTgggcggga                                                           >  1:248516/1‑144 (MQ=255)
                                                                                     gcgTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATAACTGGCAAggc                                                           <  2:103500/150‑1 (MQ=255)
                                                                                               aGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAACAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAgg                                                 >  2:112448/1‑150 (MQ=255)
                                                                                                    tCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGt                                            >  2:815331/1‑150 (MQ=255)
                                                                                                          gTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATc                                      >  1:948780/1‑150 (MQ=255)
                                                                                                                       cGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTCACATCAGCAAGCGGAtga                         <  2:15880/150‑3 (MQ=255)
                                                                                                                        gTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGAGCGGATCAg                        >  1:597001/1‑150 (MQ=255)
                                                                                                                         tAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGTCGGATCAgg                       <  1:1082924/150‑1 (MQ=255)
                                                                                                                            ggTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCa                     <  1:881991/149‑1 (MQ=255)
                                                                                                                            ggTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCa                     >  2:881991/1‑149 (MQ=255)
                                                                                                                             gTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGt                   <  2:473995/150‑1 (MQ=255)
                                                                                                                                    tCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCAcc            >  2:601113/1‑150 (MQ=255)
                                                                                                                                    gtttgctGCCGCCGGGCGGCAACAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCAcc            >  1:670613/7‑150 (MQ=255)
                                                                                                                                        gTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAgcg        >  2:1051532/1‑150 (MQ=255)
                                                                                                                                              gccCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCt  <  2:434312/150‑1 (MQ=255)
                                                                                                                                              gccCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCt  >  2:708068/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCGAGAACCTTGGCTACCAGCCGCTCGCCCTCATCGGTAAAGCGCAGCTGCCGCGTACTGCGCTCCACCAGGCGAATGCCCAGCCGCGTCTCTAAAGCACTCAGGCGTTTCGACAGCACCGTAGGGTGGCGCTCCAGTTGCCGCCCGGCGGCAGCAAACGAGCCCTGCTCATACAAGGCCAGCAGGGTGGCAAGTTCATCAGCGCGACGACTGTCGAAGGTATCCATGGCAAGGCTCGGCTCAGGGCAGTGAAATCAGGGGGCGGATCAGGCAGTCAGCACCAGCGCGCCCT  >  NC_002947/4097493‑4097784

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