Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,816,637 T→A 12.9% I21N (ATT→AAT)  PP_2472 → MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,816,6370TA12.9% 103.9 / 6.9 40I21N (ATT→AAT) PP_2472MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (19/15);  new base A (2/3);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTGCG  >  NC_002947/2816489‑2816785
                                                                                                                                                    |                                                                                                                                                    
ccAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAtt                                                                                                                                                     <  1:186479/150‑1 (MQ=255)
 cAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAATg                                                                                                                                                    >  1:253910/1‑150 (MQ=255)
          ggAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAATGGTGAAGtag                                                                                                                                           <  2:740542/150‑3 (MQ=255)
           gAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAATGGTGAAGtagc                                                                                                                                          <  2:253910/150‑4 (MQ=255)
             agagATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCg                                                                                                                                        >  2:126031/1‑150 (MQ=255)
                   cccGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgt                                                                                                                                  <  2:230039/150‑1 (MQ=255)
                          aCCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGt                                                                                                                           >  1:281051/1‑150 (MQ=255)
                           ccGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTa                                                                                                                          >  2:711130/1‑150 (MQ=255)
                                   cgcgTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCcaca                                                                                                                  <  2:281051/150‑1 (MQ=255)
                                       tcgtcACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGtg                                                                                                              <  2:168273/150‑1 (MQ=255)
                                              cTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGt                                                                                                       >  1:734216/1‑150 (MQ=255)
                                               tttCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTGTGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTa                                                                                                      >  1:813014/1‑150 (MQ=255)
                                               tttCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTa                                                                                                      >  2:305430/1‑150 (MQ=255)
                                                   gTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAAGCGCAGCAAGCGTTGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                                  <  1:787265/150‑1 (MQ=255)
                                                   gTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                                  <  1:280083/150‑1 (MQ=255)
                                                      cAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAg                                                                                               <  2:813014/150‑1 (MQ=255)
                                                         ggCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAg                                                                                            <  2:766641/150‑1 (MQ=255)
                                                                 ttGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGa                                                                                    <  1:305430/150‑1 (MQ=255)
                                                                       gCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCt                                                                              >  1:83848/1‑150 (MQ=255)
                                                                         ccGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGGCGTAAAGCCACACGTGCTGGGGTATTGGGGGCAGGAGTTCGAGCAGCTGa                                                                            >  1:181078/1‑150 (MQ=255)
                                                                         ccGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGa                                                                            >  1:214562/1‑150 (MQ=255)
                                                                         ccGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGGTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCtca                                                                            >  2:751592/1‑148 (MQ=255)
                                                                             gTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAAccc                                                                        <  2:214562/150‑1 (MQ=255)
                                                                                          cTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGcggcgg                                                           >  2:142401/1‑150 (MQ=255)
                                                                                           tGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCggg                                                          >  2:258376/1‑150 (MQ=255)
                                                                                                      cATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACcggcgg                                               <  1:142401/150‑1 (MQ=255)
                                                                                                       aTAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGt                                              <  2:529574/150‑1 (MQ=255)
                                                                                                        tAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTa                                             >  1:25728/1‑150 (MQ=255)
                                                                                                            gacgaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTAc                                         <  2:734216/150‑1 (MQ=255)
                                                                                                                 aCTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGtt                                                                                       <  1:17557/99‑1 (MQ=255)
                                                                                                                 aCTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGtt                                                                                       >  2:17557/1‑99 (MQ=255)
                                                                                                                  ctccctcccATACCGGGCAAACGCTACTTCACCACCGGTGAAATGGGCTGGTTGTGTGCCGTAAAGCAACACCTAGGTTGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAgcgc                                   <  1:192193/144‑1 (MQ=255)
                                                                                                                       ccccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCggg                                                          >  2:368101/1‑122 (MQ=255)
                                                                                                                       ccccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCggg                                                          <  1:368101/122‑1 (MQ=255)
                                                                                                                         ccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGCTGTGTCCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGGATTACCAGCCCCAGGAtg                            >  2:569530/1‑150 (MQ=255)
                                                                                                                              ccGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGGTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtg                       >  2:845099/1‑150 (MQ=255)
                                                                                                                                     aaCGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATcc                >  1:280783/1‑150 (MQ=255)
                                                                                                                                       gctcaCTTCACCAATTGTGAAGTAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtga                      >  2:673088/5‑141 (MQ=255)
                                                                                                                                       gctcaCTTCACCAATGGTGAAGTAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtga                      <  1:673088/137‑1 (MQ=255)
                                                                                                                                        gCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGcc             >  1:801036/1‑150 (MQ=255)
                                                                                                                                                   aTTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTgcg  <  2:280783/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGTGCG  >  NC_002947/2816489‑2816785

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: