Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,863,777 G→C 10.3% P900P (CCG→CCC)  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,7770GC10.3% 109.6 / 5.0 39P900P (CCG→CCC) PP_3413two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (16/19);  new base C (2/2);  total (18/21)
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.10e-01

TCGACGACTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAAT  >  NC_002947/3863629‑3863919
                                                                                                                                                    |                                                                                                                                              
tCGACGACTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCg                                                                                                                                                <  7:70186/149‑1 (MQ=255)
  gacgacTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGt                                                                                                                                              >  6:64933/1‑149 (MQ=255)
   acgacTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGt                                                                                                                                              >  4:143635/1‑148 (MQ=255)
           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)
                                     gCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCAt                                                                                                                                           >  6:290079/1‑117 (MQ=255)
                                     gCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCAt                                                                                                                                           <  5:290079/117‑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                                                                                     <  2:240178/98‑1 (MQ=255)
                                                                                                              gcgcgcTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCt                                                                                     >  1:240178/1‑98 (MQ=255)
                                                                                                                       gACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCca                                                                                 <  7:238327/93‑1 (MQ=255)
                                                                                                                       gACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCca                                                                                 >  8:238327/1‑93 (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)
                                                                                                                                                    |                                                                                                                                              
TCGACGACTACCCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAAT  >  NC_002947/3863629‑3863919

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