Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 651,886 G→C 10.9% G143G (GGG→GGC citA → sensory histidine kinase in two‑component regulatory system with citB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb651,8860GC10.9% 141.2 / 7.0 55G143G (GGG→GGCcitAsensory histidine kinase in two‑component regulatory system with citB
Reads supporting (aligned to +/- strand):  ref base G (21/28);  new base C (2/4);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

CAGTTCACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGC  >  W3110S.gb/651818‑651955
                                                                    |                                                                     
cAGTTCACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGa                                                                      <  1:532320/70‑1 (MQ=255)
cAGTTCACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGAt                                                                     >  1:59396/1‑71 (MQ=255)
     cACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCg                                                                 <  1:3936498/70‑1 (MQ=255)
       ccAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATg                                                              >  1:4327820/1‑71 (MQ=255)
         aaaCAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATgcg                                                            >  1:2594893/1‑71 (MQ=255)
          aaCAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATgcgc                                                           >  1:4143098/1‑71 (MQ=255)
                     cTGGAGAAAGGGGAGAGCTACTTCATTACCTGGAAAGGGTCAATGGGGATGGc                                                                  <  1:298342/53‑1 (MQ=255)
                      tGGAGAATGGGGAGAGCTTCTTCATTACCGGGAAAGGGTCAATGGGGATGGc                                                                  <  1:3869464/52‑1 (MQ=39)
                       ggAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCaa                                                      <  1:4278288/63‑1 (MQ=255)
                       ggAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCaa                                                      <  1:1158716/63‑1 (MQ=255)
                        gagaAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCaa                                              >  1:2677812/1‑70 (MQ=255)
                        gagaAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCa                                               <  1:2613115/69‑1 (MQ=255)
                        gagaAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAAt                                             >  1:1392204/1‑71 (MQ=255)
                            aaGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCttt                                         <  1:1929513/71‑1 (MQ=255)
                                 gagaGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                    >  1:1359440/1‑71 (MQ=255)
                                 gagaGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                    <  1:1646040/71‑1 (MQ=255)
                                 gagaGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                    >  1:3591208/1‑71 (MQ=255)
                                  agagCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGAc                                   <  1:1865159/71‑1 (MQ=255)
                                  agagCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGAc                                   <  1:3328053/71‑1 (MQ=255)
                                      cTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGa                                                                >  1:2893388/1‑38 (MQ=255)
                                       tACTTCATTACCGGGAAACGGTCAATGGGGATGGc                                                                  <  1:268302/35‑1 (MQ=255)
                                        aCTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtg                                        >  1:2765414/1‑60 (MQ=255)
                                         cTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCttt                                         <  1:1817628/58‑1 (MQ=255)
                                         cTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATgg                              <  1:1153923/69‑1 (MQ=255)
                                         cTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGaa                            <  1:1530816/71‑1 (MQ=255)
                                         cTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGaa                            <  1:1464835/71‑1 (MQ=255)
                                          ttCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCaa                                                      <  1:2279313/44‑1 (MQ=255)
                                            cATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgat                                      >  1:2582243/1‑58 (MQ=255)
                                            cATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAg                         >  1:2909059/1‑71 (MQ=255)
                                             aTTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCtt                                          >  1:2502901/1‑53 (MQ=255)
                                               tACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCa                      >  1:894166/1‑71 (MQ=255)
                                                aCCGGGAAAGTGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCAt                     >  1:440737/1‑71 (MQ=255)
                                                aCCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtga                                       >  1:2353999/1‑53 (MQ=255)
                                                  cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                       <  1:613953/51‑1 (MQ=255)
                                                  cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                       <  1:1401705/51‑1 (MQ=255)
                                                  cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                       <  1:2429994/51‑1 (MQ=255)
                                                  cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                       <  1:474005/51‑1 (MQ=255)
                                                   gggAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCgg                  >  1:574689/1‑71 (MQ=255)
                                                      aaaGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCgt               >  1:3586125/1‑71 (MQ=255)
                                                      aaaGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCgt               <  1:607385/71‑1 (MQ=255)
                                                       aaGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                    <  1:4073854/49‑1 (MQ=255)
                                                       aaGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCgtg              <  1:3087684/71‑1 (MQ=255)
                                                            tCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTc         <  1:514749/71‑1 (MQ=255)
                                                            tCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTc         <  1:1530638/71‑1 (MQ=255)
                                                             cAATGGGGATGGCGATGCGCGCCAAAACGCCAATc                                            <  1:1203702/35‑1 (MQ=255)
                                                             cAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGt                        >  1:813670/1‑55 (MQ=255)
                                                              aaTTGGCATGGCGATGCGTGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGa       >  1:4360123/1‑71 (MQ=255)
                                                              aaTTGGCATGGCGATGCGTGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGa       >  1:2965587/1‑71 (MQ=255)
                                                              aaTGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGaa                            <  1:2307800/50‑1 (MQ=255)
                                                              aaTGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTc         <  1:4142023/69‑1 (MQ=255)
                                                              aaTGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCg        <  1:3850860/70‑1 (MQ=255)
                                                              aaTGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGa       >  1:206015/1‑71 (MQ=255)
                                                               aTGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGaaa                           <  1:1116836/50‑1 (MQ=255)
                                                                tGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGa       <  1:276059/69‑1 (MQ=255)
                                                                tGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGAtt     >  1:989430/1‑71 (MQ=255)
                                                                 ggggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATc                    >  1:4258512/1‑55 (MQ=255)
                                                                  gtgATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCgtg              >  1:3784698/3‑60 (MQ=255)
                                                                  gggATGGCGATGCGCTCCAAAACGCCAATCTTTGATGACGATg                               >  1:1468003/1‑43 (MQ=255)
                                                                  gggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATg                               >  1:2746529/1‑43 (MQ=255)
                                                                  gggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATg                               >  1:518541/1‑43 (MQ=255)
                                                                  gggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGaaaa                          >  1:1916827/1‑48 (MQ=255)
                                                                  gggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCAt                     >  1:4183670/1‑53 (MQ=255)
                                                                  gggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTgg   >  1:3047458/1‑71 (MQ=255)
                                                                   ggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGc  <  1:1575724/71‑1 (MQ=255)
                                                                   ggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGc  <  1:3953917/71‑1 (MQ=255)
                                                                    |                                                                     
CAGTTCACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGC  >  W3110S.gb/651818‑651955

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