Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,863,781 A→C 10.8% M902L (ATG→CTG) ‡ 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,7810AC10.8% 102.8 / 5.9 37M902L (ATG→CTG) ‡PP_3413two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (14/19);  new base C (2/2);  total (16/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 = 9.43e-01

CCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAA  >  NC_002947/3863640‑3863929
                                                                                                                                             |                                                                                                                                                    
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)
                                                                                                                                             aTGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCGGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCGTGGCGTGACCGCCACTGCCCCGGCCGAGGAGCGGACTCGCGGACTGGCCAGCGGCATGGACGAATGGCTGTTctc  >  4:243252/1‑147 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCGGCCAACCTGCTTTTGCTGGAACGGCAACTGCAAACCCTGGGGCACCATGTGACACTGGCGGAGAATGGCGAGATCGCCCTTGCACGCTGGCAGGAAGCGCGCTTCGACCTGGTGATCACCGATTGCAGCATGCCGGTCATGGACGGGCATGAACTCACCCGGCGCATTCGCAGCCTGGAGGGCGAGCGCGGCCTGCCAGCAGGCCGTATTCTTGGCGTGACCGCCAATGCCCAGGCCGAGGAGCGTACTCGCTGCCTGGCCAGCGGCATGGACGAATGCCTGTTCAA  >  NC_002947/3863640‑3863929

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