Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,863,787 G→C 11.2% G904R (GGG→CGG)  PP_3413 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,863,7870GC11.2% 94.9 / 5.1 36G904R (GGG→CGG) PP_3413two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (15/16);  new base C (2/2);  total (17/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 = 7.67e-01

CCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGAT  >  NC_002947/3863640‑3863935
                                                                                                                                                   |                                                                                                                                                    
ccGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACgg                                                                                                                                                     <  8:89166/149‑1 (MQ=255)
    ccAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTTGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCAt                                                                                                                                                 <  3:206329/149‑1 (MQ=255)
     cAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCAt                                                                                                                                                 >  2:134929/1‑148 (MQ=255)
       aCCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgct                                                                                                                                              <  3:226450/149‑3 (MQ=255)
               ttGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCAcc                                                                                                                                        <  3:240361/147‑1 (MQ=255)
                  cTGGAACGGCAACTGCAAACCCTGGTGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGgc                                                                                                                                    <  2:264265/148‑1 (MQ=255)
                     gAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCAt                                                                                                                                >  4:60076/1‑149 (MQ=255)
                     gAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgctg                                                                                                                                             <  6:324754/136‑4 (MQ=255)
                     gAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCCGTCCATGACCGGCATgctg                                                                                                                                             >  5:324754/1‑133 (MQ=255)
                      aaCGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCAtt                                                                                                                               <  5:64933/149‑1 (MQ=255)
                       aCGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTc                                                                                                                              <  7:24330/149‑1 (MQ=255)
                                            gCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGGGAGCACCGATTGCAGCATGCCGGCCATGGACTGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGcgcg                                                                                                         <  1:7795/149‑1 (MQ=255)
                                                        acTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGCCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCcagcag                                                                                             >  2:245903/1‑149 (MQ=255)
                                                                       cGAGATCGCCCTTGCACTCTGGCAGGAAGCGCGCTTCGACCTTGTGATCAACGATTGCAGCTTGGCGGTCATGGAGGGGCATGAACTCACCCGGCGCATTCGCAGCCTGTAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGc                                                                               <  5:36001/148‑1 (MQ=255)
                                                                       cGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTAGACCTGGTGATCACCGATTGCAGCATGCCGTTCATGGACGGGCATGAACTCACCAGGCGCATTCGCAGACTGGAGGGAGAGCGCGGCCTGCCAGCAGGCCCTATTATTGGc                                                                               >  6:36001/1‑148 (MQ=255)
                                                                          gaTCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGa                                                                           <  1:134929/149‑1 (MQ=255)
                                                                              gCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGc                                                                       <  3:301636/149‑1 (MQ=255)
                                                                                       cGCTGGCAGGAAGCGCGCTTCGACCCGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGcc                                                                >  5:62948/1‑147 (MQ=255)
                                                                                            gCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGcc                                                          >  5:140738/1‑148 (MQ=255)
                                                                                                 aaGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCAAGAACTCACCCGGCGCATTCACAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCgaggag                                                    >  1:21331/1‑149 (MQ=255)
                                                                                                  aGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGTCCGTATTCTTGGCGTGACCGCCAATTCCCAGGCCgaggag                                                    <  4:152472/148‑1 (MQ=255)
                                                                                                  aGCGCGCTGCGACCTGGTGATCACCGATTGCAGCAGGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGGGACCGCCAATGCCCAGGCCGAGGAGc                                                   <  8:273030/149‑1 (MQ=255)
                                                                                                   gcgcgcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCt                                                                                                     >  1:240178/1‑98 (MQ=255)
                                                                                                   gcgcgcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCt                                                                                                     <  2:240178/98‑1 (MQ=255)
                                                                                                            gACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCca                                                                                                 >  8:238327/1‑93 (MQ=255)
                                                                                                            gACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCca                                                                                                 <  7:238327/93‑1 (MQ=255)
                                                                                                             aCCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGc                                        <  3:60076/149‑1 (MQ=255)
                                                                                                             aCCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGc                                        <  3:143635/149‑1 (MQ=255)
                                                                                                                      tCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGc                               >  8:33635/1‑149 (MQ=255)
                                                                                                                       cACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCg                              >  1:76126/1‑149 (MQ=255)
                                                                                                                         ccGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCgg                             >  3:183872/1‑148 (MQ=255)
                                                                                                                           gATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCAt                          <  7:33635/149‑1 (MQ=255)
                                                                                                                            aTTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATg                         <  6:140738/149‑1 (MQ=255)
                                                                                                                                cagcaTGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGAc                      >  4:251047/1‑148 (MQ=255)
                                                                                                                                   caTGCCCGTCCATGACCGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAAt                  >  4:283768/1‑149 (MQ=255)
                                                                                                                                             aTGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCGGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCGTGGCGTGACCGCCACTGCCCCGGCCGAGGAGCGGACTCGCGGACTGGCCAGCGGCATGGACGAATGGCTGTTctc        >  4:243252/1‑147 (MQ=255)
                                                                                                                                              tGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAg       >  1:68938/1‑149 (MQ=255)
                                                                                                                                                   gggCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCAAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGAt  <  1:245903/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAAGCCGAT  >  NC_002947/3863640‑3863935

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