Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,866,753 A→C 13.0% intergenic (+20/+37) PP_5140 → / ← thyA MerR family transcriptional regulator/thymidylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,866,7530AC13.0% 118.8 / 13.5 54intergenic (+20/+37)PP_5140/thyAMerR family transcriptional regulator/thymidylate synthase
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base C (4/3);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

GTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCA  >  NC_002947/5866605‑5866892
                                                                                                                                                    |                                                                                                                                           
ggggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGAGCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTGGCGGGTGaa                                                                                                                                            >  2:1015538/3‑150 (MQ=255)
           gcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAagag                                                                                                                                 >  1:591669/1‑150 (MQ=255)
             ggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGt                                                                                                                               <  2:230621/150‑1 (MQ=255)
             ggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGt                                                                                                                               <  2:85041/150‑1 (MQ=255)
                  tGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGtt                                                                                                                          >  1:1089617/1‑150 (MQ=255)
                   gCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTc                                                                                                                         <  1:741976/150‑1 (MQ=255)
                     aaCGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTa                                                                                                                       >  2:979819/1‑150 (MQ=255)
                       cgcgCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAAc                                                                                                                     <  2:646645/150‑1 (MQ=255)
                             ggggAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga                                                                                                               <  1:1078250/150‑1 (MQ=255)
                             ggggAGGAAAGGCGGGATCTGGCAGAAACTGGAGACAAACGGGGCGGTATAGGTACCGGATATCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga                                                                                                               <  2:107739/150‑1 (MQ=255)
                                gAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGAt                                                                                                            >  1:1061616/1‑150 (MQ=255)
                                 aGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATag                                                                                                                <  1:1173034/145‑1 (MQ=255)
                                 aGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATag                                                                                                                >  2:1173034/1‑145 (MQ=255)
                                        gcgcgcTCTTGCAGCAACTGGAGGCAAACGTGGCGGTATCGGTACCGGATACCGGACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCGTCGCGGGTGACCCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGa                                                                                                    <  1:891446/150‑1 (MQ=255)
                                               cTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAt                                                                                             <  1:448430/150‑1 (MQ=255)
                                                     gcaACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGgcgc                                                                                       <  2:1061616/150‑1 (MQ=255)
                                                      caACTGGAGACAAACGGGGCGGGATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCt                                                                                      <  2:488422/150‑1 (MQ=255)
                                                                      gggCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTAAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGACGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGt                                                                      <  1:428820/150‑1 (MQ=255)
                                                                           gTATCCGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCa                                                                 <  1:1022944/150‑1 (MQ=255)
                                                                            tATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCAt                                                                >  2:1179270/1‑150 (MQ=255)
                                                                                ggTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa                                                                                                                                    <  1:694499/78‑1 (MQ=255)
                                                                                ggTAACGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa                                                                                                                                    >  2:694499/1‑78 (MQ=39)
                                                                                 gTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCt                                                           <  2:1089617/150‑1 (MQ=255)
                                                                                            cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcca                                                       <  2:195585/143‑1 (MQ=255)
                                                                                            cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTcca                                                       >  1:195585/1‑143 (MQ=255)
                                                                                              aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgc                                                                                                                                       <  2:758816/61‑1 (MQ=39)
                                                                                              aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgc                                                                                                                                       <  2:758793/61‑1 (MQ=39)
                                                                                              aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgc                                                                                                                                       >  1:758816/1‑61 (MQ=255)
                                                                                              aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCcgc                                                                                                                                       >  1:758793/1‑61 (MQ=255)
                                                                                              aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga                                                                                                                                  <  1:796363/66‑1 (MQ=255)
                                                                                              aaCATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga                                                                                                                                  >  2:796363/1‑66 (MQ=255)
                                                                                                     ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagcca                                                   <  1:38532/138‑1 (MQ=255)
                                                                                                     ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagcca                                                   >  2:38532/1‑138 (MQ=255)
                                                                                                     ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCgg                                       >  1:171555/1‑150 (MQ=255)
                                                                                                        tgtgGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTc                                    >  2:329631/1‑150 (MQ=255)
                                                                                                           ggggCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAt                                 <  1:329631/150‑1 (MQ=255)
                                                                                                                   gTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACTCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCGTCGGCGCTGTCATCGGCGCATGGTGCTCATCGCCTTCCAGCCAGAAATCACTCGGCTCAAGCTTCTCCa                         <  1:1015538/150‑1 (MQ=255)
                                                                                                                    tCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAg                        >  1:253139/1‑150 (MQ=255)
                                                                                                                     cACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGc                       <  1:979819/150‑1 (MQ=255)
                                                                                                                      aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAGCCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAtcttc                             <  1:67111/143‑1 (MQ=255)
                                                                                                                      aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAtcttc                             >  2:67111/1‑143 (MQ=255)
                                                                                                                      aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGcc                      >  2:1271967/1‑150 (MQ=255)
                                                                                                                      aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTAGAGcc                      >  1:983337/1‑150 (MQ=255)
                                                                                                                            aTTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCg                >  2:814057/1‑150 (MQ=255)
                                                                                                                             ttGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCgg               >  1:180333/1‑150 (MQ=255)
                                                                                                                                    aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc        >  2:208873/1‑150 (MQ=255)
                                                                                                                                    aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc        >  2:208786/1‑150 (MQ=255)
                                                                                                                                      aaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACtt      <  2:253139/150‑1 (MQ=255)
                                                                                                                                       ctcTTCGCGGGTTCACCCGCGAAGCGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTg     >  1:232300/3‑150 (MQ=255)
                                                                                                                                        gCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc    >  1:528879/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg                                        <  2:1028941/112‑1 (MQ=255)
                                                                                                                                          ttCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCg                                        >  1:1028941/1‑112 (MQ=255)
                                                                                                                                          ttCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCACTCTTCTCCAGCCACTCCGGCTCATACTTGCCa  >  2:25593/1‑150 (MQ=255)
                                                                                                                                                 gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt                                                                                                        <  1:256408/41‑1 (MQ=255)
                                                                                                                                                 gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt                                                                                                        >  2:256408/1‑41 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
GTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCATCGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCA  >  NC_002947/5866605‑5866892

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