Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,532,597 G→C 9.5% A210G (GCC→GGC)  PP_3121 ← transcriptional regulator, LysR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,532,5970GC 9.5% 128.7 / 4.9 42A210G (GCC→GGC) PP_3121transcriptional regulator, LysR family
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base C (2/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATG  >  NC_002947/3532449‑3532743
                                                                                                                                                    |                                                                                                                                                  
gACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCAGGCGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGATTTCGCAAGCGGATGCTTGATTACCGgggcgggc                                                                                                                                                   >  1:866610/1‑150 (MQ=255)
        tgccGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGCCCGCcacgg                                                                                                                                           <  1:513768/147‑3 (MQ=255)
            gCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGgccgc                                                                                                                                                >  2:1071040/1‑141 (MQ=255)
            gCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGgccgc                                                                                                                                                <  1:1071040/141‑1 (MQ=255)
                      cccAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGTCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCTACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATcagc                                                                                                                             <  2:359316/150‑1 (MQ=255)
                              cTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACGATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCAc                                                                                                                     >  1:329897/1‑150 (MQ=255)
                              cTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCAc                                                                                                                     >  2:1018339/1‑150 (MQ=255)
                                     tGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGaa                                                                                                              >  1:398106/1‑150 (MQ=255)
                                      ggTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAg                                                                                                             >  1:150496/1‑150 (MQ=255)
                                       gTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAgg                                                                                                            >  1:277842/1‑150 (MQ=255)
                                       gTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAgg                                                                                                            >  1:595384/1‑150 (MQ=255)
                                          gATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCCACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGTCCCGCCGaacttcg                                                                                                         <  2:866610/150‑6 (MQ=255)
                                                caGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAg                                                                                                   <  2:750275/150‑1 (MQ=255)
                                                 aGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAgg                                                                                                  >  1:191454/1‑150 (MQ=255)
                                                 aGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAgg                                                                                                  <  2:191454/150‑1 (MQ=255)
                                                       ggcacggcacTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGgcgc                                                                                                        <  2:188501/138‑1 (MQ=255)
                                                       ggcacggcacTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGgcgc                                                                                                        >  1:188501/1‑138 (MQ=255)
                                                                 tGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCgg                                                                                  >  1:596838/1‑150 (MQ=255)
                                                                     ccccAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGgccgc                                                                                                                                                >  2:827369/1‑84 (MQ=255)
                                                                     ccccAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGgccgc                                                                                                                                                <  1:827369/84‑1 (MQ=255)
                                                                      cccAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATc                                                                             >  1:1197219/1‑150 (MQ=255)
                                                                            ccGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCATGGCCACGAAATCTTCGGCCATCAGCtgt                                                                       <  1:642001/150‑1 (MQ=255)
                                                                                tGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGgc                                                                                                                         <  2:647599/96‑1 (MQ=255)
                                                                                tGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGgc                                                                                                                         >  1:647599/1‑96 (MQ=255)
                                                                                  gCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCCACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCATGACCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAg                                                                 <  1:688472/150‑1 (MQ=255)
                                                                                    caccaTGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCTACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCt                                                               >  1:958543/1‑150 (MQ=255)
                                                                                    caccaTGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGGGCACAGCt                                                               >  2:882801/1‑150 (MQ=255)
                                                                                         tGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGGGCACAGCTGatgg                                                          <  1:882801/150‑4 (MQ=255)
                                                                                            ggCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTggg                                                       >  2:1173665/1‑150 (MQ=255)
                                                                                              ccGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCATGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGgc                                                     >  1:624665/1‑150 (MQ=255)
                                                                                                 aTGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGgc                                                                                                                         >  1:881836/1‑79 (MQ=255)
                                                                                                 aTGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGgc                                                                                                                         <  2:881836/79‑1 (MQ=255)
                                                                                                    gtcgtcAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGt                                               <  2:1197219/150‑1 (MQ=255)
                                                                                                    gtcgtcAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGt                                               <  2:29596/150‑1 (MQ=255)
                                                                                                     tcgtcAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTg                                              >  1:727639/1‑150 (MQ=255)
                                                                                                             tggtggCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGc                                      <  2:727639/150‑1 (MQ=255)
                                                                                                                 ggCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGt                                  >  1:903989/1‑150 (MQ=255)
                                                                                                                    tttCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTcgg                               >  2:760872/1‑150 (MQ=255)
                                                                                                                                  gCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTggg                                                       >  1:813924/1‑112 (MQ=255)
                                                                                                                                  gCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTggg                                                       <  2:813924/112‑1 (MQ=255)
                                                                                                                                   ccAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCaa                <  2:958543/150‑1 (MQ=255)
                                                                                                                                       cagaatgTGGCGGCCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGaa            >  1:512804/7‑150 (MQ=255)
                                                                                                                                        ttGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGt                                               <  2:966450/114‑1 (MQ=255)
                                                                                                                                        ttGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGt                                               >  1:966450/1‑114 (MQ=255)
                                                                                                                                         gaatgTGGCGGCCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGTAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGc          >  2:950787/5‑150 (MQ=255)
                                                                                                                                          aatgTGGCGGCCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCg         <  2:512804/147‑1 (MQ=255)
                                                                                                                                             gTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTc      <  2:903989/150‑1 (MQ=255)
                                                                                                                                                 ggggCCGCCACATTCTGCTCCATCGGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCAtt  <  2:1017993/149‑2 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGACCCCCAGGCCGCTGGCCACCATGCGGCCGATGGTCGTCAGCTGGTGGCTTTCGAACGCGACTGCCAGCTTGCCGTGGCGGGCCGCCACATTCTGCTCCATCAGCAGGCGCACCGCCGAAGGGCGCTGCAGGGCCACGAAATCTTCGGCCATCAGCTGTGCCCAGCTGACCTGGGGCAGCAGTGCCAAGGGCGAGTCGGCGGCCACTACGGCAACGAAGCGGTCCATG  >  NC_002947/3532449‑3532743

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