Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,221,043 T→G 12.5% D90A (GAC→GCC)  dctD‑II ← C4‑dicarboxylate transport transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,221,0430TG12.5% 117.7 / 9.2 48D90A (GAC→GCC) dctD‑IIC4‑dicarboxylate transport transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (19/23);  new base G (2/4);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GTTCGGCCAGCTGGCGGCGCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCGGCAAAATCGTCGCCGATGCGTTCCAGCG  >  NC_002947/1220894‑1221190
                                                                                                                                                     |                                                                                                                                                   
gTTCGGCCAGCTGGCGGCGCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGt                                                                                                                                                     <  2:491353/150‑1 (MQ=255)
   cGGCCAGCTGGCGGCGCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGa                                                                                                                                                  <  2:35381/150‑1 (MQ=255)
       cAGCTGGCGGCGCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATc                                                                                                                                              <  2:79083/150‑1 (MQ=255)
           tGGCGGCGCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCa                                                                                                                                          <  2:267099/150‑1 (MQ=255)
            ggcggcGCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCAt                                                                                                                                         <  1:492359/150‑1 (MQ=255)
            ggcggcGCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGGCGATATCGACAt                                                                                                                                         <  1:911150/150‑1 (MQ=255)
               ggcgCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGc                                                                                                                                      >  1:1000895/1‑150 (MQ=255)
               ggcgCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGc                                                                                                                                      <  2:429991/150‑1 (MQ=255)
                   cAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGt                                                                                                                                  <  1:685403/150‑1 (MQ=255)
                    aGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGGCGATATCGACATGGccgttg                                                                                                                                 <  1:771745/150‑5 (MQ=255)
                          accacTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGc                                                                                                                           >  1:538991/1‑150 (MQ=255)
                                    gTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCa                                                                                                                 >  1:107263/1‑150 (MQ=255)
                                          gCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGc                                                                                                           <  1:57729/150‑1 (MQ=255)
                                            ccGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCgg                                                                                                         >  1:95251/1‑150 (MQ=255)
                                                gcggTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCAc                                                                                                     >  1:1205728/4‑150 (MQ=255)
                                                  tGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATAGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACgg                                                                                                   <  1:580192/150‑1 (MQ=255)
                                                                gcgTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGAtt                                                                                     <  1:113128/150‑1 (MQ=255)
                                                                  gTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGa                                                                                   <  2:1000895/150‑1 (MQ=255)
                                                                   tACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGag                                                                                  <  2:95251/150‑1 (MQ=255)
                                                                    aCCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGagc                                                                                 >  2:1026450/1‑150 (MQ=255)
                                                                            gACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAg                                                                         <  1:1013554/150‑1 (MQ=255)
                                                                             aCCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAgg                                                                        >  2:1048778/1‑150 (MQ=255)
                                                                                 ggCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCAt                                                                    <  1:693173/150‑1 (MQ=255)
                                                                                 ggCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCAt                                                                    <  1:464765/150‑1 (MQ=255)
                                                                                  gcgcTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATc                                                                   >  1:1029145/1‑150 (MQ=255)
                                                                                      tCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCACTGCTGAGCTACAATCACTAGCGTCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATa                                                               >  1:1177868/1‑150 (MQ=255)
                                                                                          gggAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCg                                                           <  1:1179637/150‑1 (MQ=255)
                                                                                           ggAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATAccc                                                            <  1:234567/148‑1 (MQ=255)
                                                                                           ggAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATAccc                                                            >  2:234567/1‑148 (MQ=255)
                                                                                              gaACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGCGCCTTGAGGCGATTGAGCAGATCCCGTCGATCAATACCCGGCAg                                                       <  2:107263/150‑1 (MQ=255)
                                                                                               aaCGGTTTTTCTTTTTAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAgg                                                      >  2:1066298/1‑150 (MQ=255)
                                                                                                           aTGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTg                                          >  1:310357/1‑150 (MQ=255)
                                                                                                                       gcgcCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAgg                                                                        <  1:896537/108‑1 (MQ=255)
                                                                                                                       gcgcCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAgg                                                                        >  2:896537/1‑108 (MQ=255)
                                                                                                                         gcCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGc                                                                                                           <  1:35248/71‑1 (MQ=255)
                                                                                                                         gcCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGc                                                                                                           >  2:35248/1‑71 (MQ=255)
                                                                                                                         gcCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCGGc                            >  1:646178/1‑150 (MQ=255)
                                                                                                                                  cATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCGGCAAAAtcgtc                   <  2:1177868/150‑1 (MQ=255)
                                                                                                                                   aTGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCg                              <  2:1116773/138‑1 (MQ=255)
                                                                                                                                   aTGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCg                              >  1:1116773/1‑138 (MQ=255)
                                                                                                                                   aTGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCGGCAAAAtcgtcg                  >  2:256430/1‑150 (MQ=255)
                                                                                                                                    tGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCGGCAAAATCGTCGc                 <  1:315961/150‑1 (MQ=255)
                                                                                                                                        gcCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCGGCAAAATCGTCGCCGAt             >  2:1083158/1‑150 (MQ=255)
                                                                                                                                             tggCCATGGCGATATCGACATGGCCGGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATAccc                                                            <  2:1071212/97‑1 (MQ=39)
                                                                                                                                             tggCCATGGCGATATCGACATGGCCGGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATAccc                                                            >  1:1071212/2‑98 (MQ=39)
                                                                                                                                             tggCCATGGCGATATCGACATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTga                                         <  1:466168/116‑1 (MQ=255)
                                                                                                                                             tggCCATGGCGATATCGACATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTga                                         >  2:466168/2‑117 (MQ=255)
                                                                                                                                             cGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCGGCAAAATCGTCGCCGATGCGtt        <  2:310357/150‑1 (MQ=255)
                                                                                                                                                   tGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGAAGGCGGATATCGCTGACGACAATGCCGGCAAAATCGTCGCCGATGCGTTCCAgcg  >  1:986263/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTTCGGCCAGCTGGCGGCGCAGGGCCACCACTTCGCGTGACAGCCCGCGCTGTTCCAGCGCGCGGCGTACCACGTCGACCAGGCGCTCGGGGGAGAACGGTTTTTCCATGAAGTCGTAGGCGCCGTTACGCATGGCGCCAACGGCCATGTCGATATCGCCATGGCCAGTGATCAGCACCACCGGCAGGCTGCGGTCACGGGCCTTGAGGCGATTGAGCAGCTCCAGGCCATCAATACCCGGCAGGCGGATATCGCTGACGACAATGCCGGCAAAATCGTCGCCGATGCGTTCCAGCG  >  NC_002947/1220894‑1221190

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