Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,351,047 (T)7→8 coding (616/1632 nt) dcuS ← sensory histidine kinase DcuS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,351,0401.T97.0% 114.3 / ‑3.3 33coding (623/1632 nt)dcuSsensory histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (20/12);  total (21/12)
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

GCCTTCTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAA  >  NC_000913/4350961‑4351115
                                                                                |                                                                            
gCCTTCTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGAt                                                                             >  2:826618/1‑81 (MQ=255)
   ttcttTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGAtttt                                                                         >  2:826629/1‑81 (MQ=255)
        tGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGAtttttt                                                                       <  2:826578/78‑1 (MQ=255)
         gATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGAttttttttc                                                                    <  2:826579/80‑3 (MQ=255)
         gATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGAtttttttt                                                                     >  1:826610/1‑78 (MQ=255)
          aTAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGAttttttttcag                                                                  <  1:826581/81‑5 (MQ=255)
           tAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGAttttttttcag                                                                  >  1:826627/1‑76 (MQ=255)
                 gCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCtt                                                             >  2:826620/1‑80 (MQ=255)
                  cAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTaa                                                           >  1:826617/1‑81 (MQ=255)
                  cAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTa                                                            <  1:826583/80‑1 (MQ=255)
                   aaCATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAAc                                                          >  2:826635/1‑81 (MQ=255)
                   aaCATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAAc                                                          <  1:826584/81‑1 (MQ=255)
                     cATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAAc                                                          >  2:826667/1‑79 (MQ=255)
                       tGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGa                                                      <  2:826584/81‑1 (MQ=255)
                                 ttGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGc                                            >  1:826628/1‑81 (MQ=255)
                                    cTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCa                                          <  2:826585/80‑1 (MQ=255)
                                      cAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAAt                                        >  2:826621/1‑80 (MQ=255)
                                       aaaCAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCa                                      <  2:826586/81‑1 (MQ=255)
                                         aCAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGt                                    <  2:826587/81‑1 (MQ=255)
                                           aGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTcc                                  <  1:826589/81‑1 (MQ=255)
                                             cGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGa                                >  2:826632/1‑81 (MQ=255)
                                                 gagaTTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAg                            >  1:826623/1‑81 (MQ=255)
                                                 gagaTTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAg                            >  2:826658/1‑81 (MQ=255)
                                                   gaTTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCa                          >  2:826648/1‑81 (MQ=255)
                                                     tttCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATg                        >  1:826641/1‑81 (MQ=255)
                                                       tCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGc                       <  2:826589/80‑1 (MQ=255)
                                                         gTAGGGTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCa                     >  1:826632/1‑80 (MQ=255)
                                                            gggTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAaat                  <  1:826591/80‑1 (MQ=255)
                                                             ggTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAaataa                <  1:826592/81‑1 (MQ=255)
                                                             ggTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAaataa                >  2:826634/1‑81 (MQ=255)
                                                             ggTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAaataa                >  2:826642/1‑81 (MQ=255)
                                                              gTTCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAaataat               <  1:826593/81‑1 (MQ=255)
                                                                tCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATAc             >  2:826646/1‑81 (MQ=255)
                                                                tCCAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATAc             >  1:826658/1‑81 (MQ=255)
                                                                  cAGGCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCg           <  2:826593/81‑1 (MQ=255)
                                                                    ggCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGAc         >  1:826661/1‑81 (MQ=255)
                                                                     gCCGAAAAGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGAc         >  1:826642/1‑80 (MQ=255)
                                                                     gCCGAAAAGGATTCTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGAc         >  1:826631/1‑80 (MQ=255)
                                                                         aaaaGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGa     >  1:826640/1‑80 (MQ=255)
                                                                           aaGGATTTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATaa  >  1:826638/1‑81 (MQ=255)
                                                                                |                                                                            
GCCTTCTTTGATAGACTGCAACATGGCCTGGCGTTGCTCAAACAGCGTGGAGATTTCGTAGGGTTCCAGGCCGAAAAGGATTTTTTTCAGTACCTTAACCAGAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAA  >  NC_000913/4350961‑4351115

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