Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,638,421 A→T 11.1% W121R (TGG→AGG)  PP_1437 ← heavy metal sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,638,4210AT11.1% 125.5 / 8.2 45W121R (TGG→AGG) PP_1437heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (24/16);  new base T (3/2);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

ATGTCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCA  >  NC_002947/1638283‑1638567
                                                                                                                                          |                                                                                                                                                  
aTGTCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAgg                                                                                                                                         <  2:205638/150‑1 (MQ=255)
 tGTCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGc                                                                                                                                        <  2:622608/150‑1 (MQ=255)
                  cTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCCGTCACGCCAGGCGTGAGAGGCCAGGTGCGCATCGgtgt                                                                                                                       >  2:159729/1‑150 (MQ=255)
                   ttGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGtgtg                                                                                                                      >  2:773033/1‑150 (MQ=255)
                           gggCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCAc                                                                                                              >  2:729002/1‑150 (MQ=255)
                                        ggCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTc                                                                                                 <  1:154202/150‑1 (MQ=255)
                                                               ccTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCa                                                                          >  2:388529/1‑150 (MQ=255)
                                                                cTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAg                                                                         >  1:736196/1‑150 (MQ=255)
                                                                cTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAg                                                                         <  1:731812/150‑1 (MQ=255)
                                                                cTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAg                                                                         >  2:731812/1‑150 (MQ=255)
                                                                  tGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGt                                                                       <  2:714320/150‑1 (MQ=255)
                                                                      aGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCa                                                                                               >  2:53419/1‑122 (MQ=255)
                                                                      aGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTgggg                                                                   >  2:377716/1‑150 (MQ=255)
                                                                      aGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTCTGTATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCa                                                                                               <  1:53419/122‑1 (MQ=255)
                                                                         ggTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTg                                                                <  1:388529/150‑1 (MQ=255)
                                                                             tggtCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTggcgg                                                            <  2:684470/148‑1 (MQ=255)
                                                                                   cgcAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTc                                                       >  1:548868/1‑149 (MQ=255)
                                                                                   cgcAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTc                                                       <  2:548868/149‑1 (MQ=255)
                                                                                            aTCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTCAAAGGCCAGGCGCGCATGCGAGTTCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGc                                             >  1:829161/1‑150 (MQ=255)
                                                                                                aCCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAg                                                                                  <  1:332375/109‑1 (MQ=255)
                                                                                                aCCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAg                                                                                  >  2:332375/1‑109 (MQ=255)
                                                                                                aCCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                   <  1:628439/140‑1 (MQ=255)
                                                                                                aCCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                   >  2:628439/1‑140 (MQ=255)
                                                                                                 ccGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACt                                        >  1:102547/1‑150 (MQ=255)
                                                                                                          ggTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGt                                                                                                                                               >  1:387207/1‑38 (MQ=255)
                                                                                                          ggTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGt                                                                                                                                               <  2:387207/38‑1 (MQ=255)
                                                                                                          ggTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAg                               >  2:467112/1‑150 (MQ=255)
                                                                                                             gAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGtt                            >  2:822026/1‑150 (MQ=255)
                                                                                                                  tCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCga                                                     <  2:328704/120‑1 (MQ=255)
                                                                                                                  tCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCga                                                     >  1:328704/1‑120 (MQ=255)
                                                                                                                    tGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTc                                                       <  1:815613/116‑1 (MQ=255)
                                                                                                                    tGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTc                                                       >  2:815613/1‑116 (MQ=255)
                                                                                                                      gTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGccc                   <  1:385132/150‑1 (MQ=255)
                                                                                                                             ccgcTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATcccc            >  2:368910/1‑150 (MQ=255)
                                                                                                                                ctttCACGCCTGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAgg         >  1:681569/4‑150 (MQ=255)
                                                                                                                                 tttCACGCCTGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGccc                   <  2:142829/137‑1 (MQ=255)
                                                                                                                                 tttCACGCCTGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGccc                   >  1:142829/3‑139 (MQ=255)
                                                                                                                                 tttCACGCCTGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAgg         >  1:746091/3‑149 (MQ=255)
                                                                                                                                 tttCACGCCTGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAgg         <  2:746091/147‑1 (MQ=255)
                                                                                                                                   tCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTcc      >  1:577280/1‑150 (MQ=255)
                                                                                                                                     aCGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                   <  2:119105/103‑1 (MQ=255)
                                                                                                                                     aCGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                   >  1:119105/1‑103 (MQ=255)
                                                                                                                                     aCGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCcag    >  2:130548/1‑150 (MQ=255)
                                                                                                                                     aCGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCcag    >  1:382453/1‑150 (MQ=255)
                                                                                                                                       gCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCcagca  >  2:485382/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
ATGTCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCA  >  NC_002947/1638283‑1638567

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