Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,699,942 G→A 10.5% D210N (GAT→AAT)  PP_1494 → GGDEF domain‑containing response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,699,9420GA10.5% 102.7 / 2.7 38D210N (GAT→AAT) PP_1494GGDEF domain‑containing response regulator
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base A (2/2);  total (19/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 = 6.73e-01

GATTCCAACCTGATGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATGGCGGCGAAGAGTTTGCCCTGGTG  >  NC_002947/1699795‑1700085
                                                                                                                                                   |                                                                                                                                               
gATTCCAACCTGATGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTgg                                                                                                                                                 <  4:157505/148‑1 (MQ=255)
  ttCCAACCTGATGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGAt                                                                                                                                               >  8:7477/1‑148 (MQ=255)
      aaCCTGATGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACtt                                                                                                                                          >  2:7058/1‑149 (MQ=255)
      aaCCTGATGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACtt                                                                                                                                          >  8:2817/1‑149 (MQ=255)
         cTGATGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCaa                                                                                                                                       <  6:267740/149‑1 (MQ=255)
            aTGCTGCAGCGGCTGATGACCTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCaaa                                                                                                                                      <  7:7477/147‑1 (MQ=255)
            aTGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCaaa                                                                                                                                      <  7:2817/147‑1 (MQ=255)
            aTGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCaaa                                                                                                                                      <  4:31156/147‑1 (MQ=255)
                   aGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAAc                                                                                                                             >  4:281943/1‑149 (MQ=255)
                     cGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGa                                                                                                                           <  1:7058/149‑1 (MQ=255)
                       gCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACa                                                                                                                         <  5:188849/149‑1 (MQ=255)
                              aaCTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTg                                                                                                                   <  7:54137/148‑1 (MQ=255)
                               aCTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGc                                                                                                                 <  6:202333/149‑1 (MQ=255)
                                                  gCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAgg                                                                                              >  3:231655/1‑149 (MQ=255)
                                                   cTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGAATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAgga                                                                                             >  7:70826/1‑148 (MQ=255)
                                                   cTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGAATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAg                                                                                               >  7:70822/1‑147 (MQ=255)
                                                   cTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGAATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAg                                                                                               <  8:70826/147‑1 (MQ=255)
                                                   cTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGAATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAg                                                                                               <  8:70822/147‑1 (MQ=255)
                                                       ccAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTg                                                                                         >  4:123321/1‑149 (MQ=255)
                                                         aaCCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTgcg                                                                                       >  3:4477/1‑149 (MQ=255)
                                                                  cACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGc                                                                              <  3:281943/149‑1 (MQ=255)
                                                                   aCTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGcc                                                                             >  3:186644/1‑149 (MQ=255)
                                                                                 ctggagctggagTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACa                                                                                                                               >  8:275008/1‑85 (MQ=255)
                                                                                 ctggagctggagTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACa                                                                                                                               <  7:275008/85‑1 (MQ=255)
                                                                                      gctggagTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGt                                                          >  6:247893/1‑149 (MQ=255)
                                                                                          gagTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCt                                                       <  3:160318/148‑1 (MQ=255)
                                                                                                       ccATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCcgaccgac                                         <  4:231655/149‑1 (MQ=255)
                                                                                                         aTGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCt                                       <  4:4477/149‑1 (MQ=255)
                                                                                                         aTGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCt                                       >  3:29188/1‑149 (MQ=255)
                                                                                                                aaCAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGgcg                                >  7:161319/1‑149 (MQ=255)
                                                                                                                         agTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATggc                       <  4:186644/149‑1 (MQ=255)
                                                                                                                          gTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATggcg                      >  7:186750/1‑149 (MQ=255)
                                                                                                                              tCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATGgcggcg                   <  4:29188/148‑1 (MQ=255)
                                                                                                                                gTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATGgcggcg                   >  6:170248/1‑146 (MQ=255)
                                                                                                                                 ttGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGcc                                                                                     >  3:87230/1‑79 (MQ=255)
                                                                                                                                 ttGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGcc                                                                                     <  4:87230/79‑1 (MQ=255)
                                                                                                                                 ttGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATGGCGGCGaa                 >  8:193962/1‑147 (MQ=255)
                                                                                                                                    cTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATGGCGGCGAAGAGt             <  5:247893/148‑1 (MQ=255)
                                                                                                                                              aCGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATGGCGGCGAAGAGTTTGCCCTGGtg  >  2:95740/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GATTCCAACCTGATGCTGCAGCGGCTGATGAACTCCGATGGCCTGACCGGGCTGTCCAACCGTCGCCACTTTGACGAGTACCTGGAGCTGGAGTGGCGCCGGGCCATGCGTGAACAGCAGCAGTTGTCGTTGCTGATGATCGACGTGGATTACTTCAAAGTCTACAACGACAGCTTTGGCCATCTGGCCGGTGACGAGGCCTTGCGCCAGGTGGCCGAAGCCATCCGTAGCTCGTGCTCGCGGCCGACCGACCTGCCGGCGCGCTATGGCGGCGAAGAGTTTGCCCTGGTG  >  NC_002947/1699795‑1700085

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: