Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,638,429 A→C 8.1% F118C (TTT→TGT)  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,4290AC 8.1% 154.8 / 3.7 57F118C (TTT→TGT) PP_1437heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (21/31);  new base C (2/3);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATGGG  >  NC_002947/1638282‑1638576
                                                                                                                                                   |                                                                                                                                                   
aaTGTCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAg                                                                                                                                                   <  1:252670/150‑1 (MQ=255)
    tCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTTTCACGCCTGGCGTGACAgcggc                                                                                                                                               <  2:141351/150‑5 (MQ=255)
    tCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCa                                                                                                                                               >  1:481425/1‑150 (MQ=255)
        gTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGgc                                                                                                                                            >  2:535977/1‑149 (MQ=255)
        gTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGgc                                                                                                                                            <  1:535977/149‑1 (MQ=255)
             ggTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCAt                                                                                                                                      >  2:1144506/1‑150 (MQ=255)
              gTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATc                                                                                                                                     >  2:491299/1‑150 (MQ=255)
                      gcagcaGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCACCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCa                                                                                                                             <  1:93526/150‑1 (MQ=255)
                         gcagGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGtt                                                                                                                          >  2:875475/1‑150 (MQ=255)
                                            gcgGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAg                                                                                                       <  2:490537/150‑1 (MQ=255)
                                                ttGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAgggcc                                                                                                   <  2:1269385/150‑1 (MQ=255)
                                                       aCATCAGCACCTTGAGAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGa                                                                                            <  1:17442/150‑1 (MQ=255)
                                                        catcatCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAg                                                                                           <  2:133377/150‑1 (MQ=255)
                                                         atcatcACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCa                                                                                                                             >  2:684854/1‑115 (MQ=255)
                                                         atcatcACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCa                                                                                                                             <  1:684854/115‑1 (MQ=255)
                                                         atcatcACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAgg                                                                                          <  1:466489/150‑1 (MQ=255)
                                                          tcatcaCCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGc                                                                                         >  1:803191/1‑150 (MQ=255)
                                                             tcaCCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGa                                                                                      >  1:1302130/1‑150 (MQ=255)
                                                                     aCAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTgg                                                                              <  1:364075/150‑1 (MQ=255)
                                                                         gggTGTTGTCGCGCAGGCGGTTCAACCGGGTAACGGTGAGGTTCTGGTTACTGCCGCTGTTACGGCAGGCGTGAAAGGCGAGGAGCGTATTGGTGTGCGGTTCACCGATGCGCGATTCCAGGGCGGGGGCGAGGCGGGGCAGGTGGGGGt                                                                          >  2:135451/1‑150 (MQ=255)
                                                                             gttgtCGCGCAGGCGCATCAACCGTGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTggcg                                                                      <  1:457621/150‑1 (MQ=255)
                                                                             gttgtCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTggcg                                                                      <  1:808361/150‑1 (MQ=255)
                                                                             attgtCGCGCAGACGCATCAACCGCGTAACGGTGAGGATCTCGTTACCGCCGCTGTCACGCCAGGCGCGAAAGGCCAGGCGCGCATCGGTGCGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTggcg                                                                      <  1:135451/149‑1 (MQ=255)
                                                                                 tCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCg                                                                  >  2:836861/1‑150 (MQ=255)
                                                                                      cAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgag                                                             <  1:817794/150‑1 (MQ=255)
                                                                                       aGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                            <  1:146917/150‑1 (MQ=255)
                                                                                         gcgcATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGa                                                                                      >  1:1032403/1‑122 (MQ=255)
                                                                                         gcgcATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTCTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGa                                                                                      <  2:1032403/122‑1 (MQ=255)
                                                                                            cATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGAc                                                                                                                   >  1:84672/1‑90 (MQ=255)
                                                                                            cATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGAc                                                                                                                   <  2:84672/90‑1 (MQ=255)
                                                                                             aTCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGc                                                      >  1:91251/1‑150 (MQ=255)
                                                                                              tCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGca                                                     >  1:1373562/1‑150 (MQ=255)
                                                                                                     ggTAACGGGGAGGATCTGGTTACCGCCGCTGTCACGCCAGCCGTGAAAGGCCAGGCGCGCATCGGTGTGCCGATCACCGACCCGCTATTCCAGGGCCGGTCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAg                                              <  1:454981/150‑1 (MQ=255)
                                                                                                          cGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCGGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTa                                         <  2:481425/150‑1 (MQ=255)
                                                                                                                 gATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                            >  2:743904/1‑124 (MQ=255)
                                                                                                                 gATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                            <  1:743904/124‑1 (MQ=255)
                                                                                                                 gATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTc                                  <  2:803191/150‑1 (MQ=255)
                                                                                                                   tCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCAt                                <  1:836861/150‑1 (MQ=255)
                                                                                                                    cTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATGGGTGTGCAGTTCACCGACCAGCGATTCCAGGGCCGGGCCGAGGATGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCAAGAACAATGAGACTTAGGTTGTCATg                               >  2:205357/1‑150 (MQ=255)
                                                                                                                       gTTACCGACGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGccc                            >  1:609854/1‑150 (MQ=255)
                                                                                                                             gccgcTGTCACGCCAGGCGTGAAAGGCCAGACGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGc                                                          <  2:852200/114‑1 (MQ=255)
                                                                                                                             gccgcTGTCACGCCAGGCGTGAAAGGCCAGACGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGc                                                          >  1:852200/1‑114 (MQ=255)
                                                                                                                               cgcTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGCGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCa                    >  2:359602/1‑150 (MQ=255)
                                                                                                                               cgcTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCa                    <  1:491299/150‑1 (MQ=255)
                                                                                                                                gcTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAg                   <  1:875475/150‑1 (MQ=255)
                                                                                                                                 ctttCACGCCTGGCGTGACAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGa                                                                                            >  2:93283/4‑76 (MQ=39)
                                                                                                                                 ctttCACGCCTGGCGTGACAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGa                                                                                            <  1:93283/73‑1 (MQ=38)
                                                                                                                                   gTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTc                <  1:1014100/150‑1 (MQ=255)
                                                                                                                                      aCGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTc                                  >  2:653046/1‑129 (MQ=255)
                                                                                                                                      aCGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTc                                  <  1:653046/129‑1 (MQ=255)
                                                                                                                                         ccAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCcagcag          >  2:1274958/1‑150 (MQ=255)
                                                                                                                                          cttgCGTGACAGGCCAGGCGCGCGTCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                            <  1:797660/96‑1 (MQ=255)
                                                                                                                                          ctggCGTGACAGGCCAGGCGCGCGTCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCgaga                                                            >  2:797660/3‑99 (MQ=255)
                                                                                                                                          cAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGt         <  2:1376198/150‑1 (MQ=255)
                                                                                                                                                tGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATgg   <  1:205357/150‑1 (MQ=255)
                                                                                                                                                tGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATgg   <  2:1281872/150‑1 (MQ=255)
                                                                                                                                                 gAAAGGCCAGGCGCGCATCGGTGTGCAGTTCCACGACCGGCGTTTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATggg  <  2:1302130/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AATGTCGAGTGCAGGTAGGCTTGCAGCAGGGCGTTGTCGTCGGCGCGGTTGAGCGACATCATCACCTTGACAGGGGTGTTGTCGCGCAGGCGCATCAACCGGGTAACGGTGAGGATCTGGTTACCGCCGCTGTCACGCCAGGCGTGAAAGGCCAGGCGCGCATCGGTGTGCAGTTCACCGACCCGCGATTCCAGGGCCGGGCCGAGGCTGAGCAGGTGGGGGTGGCGGCCGTCGAGAGCCAGCACACTGAGGCTTAGGTTGTCATGCCCCATCACCAGGTCCAGCAGTGGATGGG  >  NC_002947/1638282‑1638576

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