Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,100,828 A→G 7.1% P234P (CCA→CCG PP_3608 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,100,8280AG 7.1% 161.7 / 2.1 56P234P (CCA→CCGPP_3608LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/2);  total (28/28)
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.97e-01

GTCCTGCCTGCCAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACTATCAGCAAG  >  NC_002947/4100682‑4100973
                                                                                                                                                  |                                                                                                                                                 
gTCCTGCCTGCCAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCcaca                                                                                                                                                 >  8:263793/1‑149 (MQ=255)
   ctgcctgcCAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTg                                                                                                                                              >  3:370663/1‑149 (MQ=255)
        tgccAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTg                                                                                                                                         >  8:155403/1‑149 (MQ=255)
              tGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAg                                                                                                                                   <  3:88026/149‑1 (MQ=255)
                  gATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCGCACTCCAGTGTGGCaaggcc                                                                                                                               <  4:137645/149‑6 (MQ=255)
                     atGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGt                                                                                                                            >  7:4446/1‑149 (MQ=255)
                       ggCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGc                                                                                                                          >  3:290954/1‑149 (MQ=255)
                                  aCTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCAcc                                                                                                               >  4:5391/1‑149 (MQ=255)
                                      gctggcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCGCACTCCAGTGTGGCaa                                                                                                                                   >  7:186515/1‑124 (MQ=255)
                                      gctggcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCGCACTCCAGTGTGGCaa                                                                                                                                   <  8:186515/125‑2 (MQ=255)
                                      gctggcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGAcc                                                                                                           <  2:296482/149‑1 (MQ=255)
                                          gcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGa                                                                                                        <  4:96538/148‑1 (MQ=255)
                                          gcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGa                                                                                                        <  4:290954/148‑1 (MQ=255)
                                          gcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGAt                                                                                                       >  8:59093/1‑149 (MQ=255)
                                            aCGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATcc                                                                                                     >  3:163688/1‑149 (MQ=255)
                                               cccACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGa                                                                                                  <  7:263793/149‑1 (MQ=255)
                                                ccACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGa                                                                                                  >  2:9392/1‑148 (MQ=255)
                                                    aTGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTg                                                                                             <  6:127835/149‑1 (MQ=255)
                                                      gTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGc                                                                                           >  3:341933/1‑149 (MQ=255)
                                                       tATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGc                                                                                           <  7:342873/148‑1 (MQ=255)
                                                         ttGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCa                                                                                         <  4:163688/148‑1 (MQ=255)
                                                            aTGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCg                                                                                                                    <  3:137552/118‑1 (MQ=255)
                                                            aTGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCg                                                                                                                    >  4:137552/1‑118 (MQ=255)
                                                            aTGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCcg                                                                                     >  2:54309/1‑149 (MQ=255)
                                                             tGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCg                                                                                                                    <  4:140502/117‑1 (MQ=255)
                                                             tGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCg                                                                                                                    >  3:140502/1‑117 (MQ=255)
                                                                 cATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGc                                                                                                                      >  6:199713/1‑111 (MQ=255)
                                                                   tGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCtgtg                                                                              <  5:39386/149‑1 (MQ=255)
                                                                     cgcCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTggg                                                                            <  4:341933/149‑1 (MQ=255)
                                                                           gATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGtt                                                                      <  4:370663/149‑1 (MQ=255)
                                                                             tGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTgc                                                                    <  1:361087/149‑1 (MQ=255)
                                                                                tttCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTgcgc                                                                  >  3:22703/1‑148 (MQ=255)
                                                                                 ttCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGcc                                                                 >  6:331301/1‑148 (MQ=255)
                                                                                    gATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAaca                                                             >  3:239171/1‑149 (MQ=255)
                                                                                      tGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAg                                                                                               >  4:201617/1‑113 (MQ=255)
                                                                                      ggAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAAcaca                                                           <  7:338157/148‑1 (MQ=255)
                                                                                         aaTCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACaa                                                          >  2:231748/1‑147 (MQ=255)
                                                                                         aaTCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACaa                                                          <  1:231748/147‑1 (MQ=255)
                                                                                         aaTCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACATGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGa                                                        >  3:218885/1‑149 (MQ=255)
                                                                                                 gTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTgcgc                                                                  <  5:250827/131‑1 (MQ=255)
                                                                                                 gTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTgcgc                                                                  >  6:250827/1‑131 (MQ=255)
                                                                                                      cTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGtt                                           <  1:9392/149‑1 (MQ=255)
                                                                                                         gccgccCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTcac                                        <  7:155403/149‑1 (MQ=255)
                                                                                                         gccgccCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTca                                         <  1:341700/148‑1 (MQ=255)
                                                                                                            gccCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGtt                                           >  4:226867/1‑143 (MQ=255)
                                                                                                            gccCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGtt                                           <  3:226867/143‑1 (MQ=255)
                                                                                                              ccGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACaa                                                          >  6:32541/1‑126 (MQ=255)
                                                                                                              ccGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACaa                                                          <  5:32541/126‑1 (MQ=255)
                                                                                                                    cTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTgg                             <  5:331301/149‑1 (MQ=255)
                                                                                                                           cgcCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATccc                      <  3:5391/149‑1 (MQ=255)
                                                                                                                            gcCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATcccc                     <  1:54309/149‑1 (MQ=255)
                                                                                                                                 aTTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGtt                >  5:50366/1‑149 (MQ=255)
                                                                                                                                  ttGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTc               >  2:117973/1‑149 (MQ=255)
                                                                                                                                      ccTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGAc            >  6:383133/1‑148 (MQ=255)
                                                                                                                                      ccTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACt           <  4:239171/149‑1 (MQ=255)
                                                                                                                                            tctGCCGCACTCCAGTGTGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCACTTCACACCGCCACTGGCGATCCCCTCGATCGACTATCAg      >  8:296248/3‑148 (MQ=255)
                                                                                                                                               gCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACTATCAGCAAg  >  4:342764/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTCCTGCCTGCCAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACTATCAGCAAG  >  NC_002947/4100682‑4100973

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: