Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 651,898 C→T 8.6% R147R (CGC→CGT 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,8980CT 8.6% 204.9 / 5.8 69R147R (CGC→CGTcitAsensory histidine kinase in two‑component regulatory system with citB
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base T (2/4);  total (29/40)
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

AACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGCTACCTGG  >  W3110S.gb/651828‑651962
                                                                      |                                                                
aaCAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATgcgc                                                                  >  1:4143098/1‑71 (MQ=255)
             ggAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCaa                                                             <  1:1158716/63‑1 (MQ=255)
             ggAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCaa                                                             <  1:4278288/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:3591208/1‑71 (MQ=255)
                       gagaGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                           >  1:1359440/1‑71 (MQ=255)
                       gagaGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                           <  1:1646040/71‑1 (MQ=255)
                        agagCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGAc                                          <  1:3328053/71‑1 (MQ=255)
                        agagCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGAc                                          <  1:1865159/71‑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:474005/51‑1 (MQ=255)
                                        cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                              <  1:2429994/51‑1 (MQ=255)
                                        cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                              <  1:613953/51‑1 (MQ=255)
                                        cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                              <  1:1401705/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:518541/1‑43 (MQ=255)
                                                        gggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATg                                      >  1:2746529/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:3953917/71‑1 (MQ=255)
                                                         ggATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGc         <  1:1575724/71‑1 (MQ=255)
                                                           aTGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGa                                        <  1:3769980/38‑1 (MQ=255)
                                                           aTGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGa                                        <  1:2103711/38‑1 (MQ=255)
                                                           aTGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGa                                        <  1:3416957/38‑1 (MQ=255)
                                                           aTGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGa                                        <  1:1343716/38‑1 (MQ=255)
                                                            tGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGCTAc      <  1:2580262/71‑1 (MQ=255)
                                                            tGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGCTAc      <  1:4093609/71‑1 (MQ=255)
                                                            tGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGCTAc      <  1:1821899/71‑1 (MQ=255)
                                                             ggCGATGCGCGCCAAAACGCCAATCTTTGATGACTATGGAAAAGTCATCGGCGTGGTGTCGATTGGCt        >  1:163573/1‑68 (MQ=255)
                                                             ggCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGTAAAAGTCATCGGCGTGGTGTCGa              <  1:2665441/62‑1 (MQ=255)
                                                              gCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGCTACCt    <  1:1090680/71‑1 (MQ=255)
                                                               cGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGCTACCTg   >  1:3216652/1‑71 (MQ=255)
                                                                gATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGCTACCTgg  >  1:2566168/1‑71 (MQ=255)
                                                                 aTGCGCGCCAAAACGCCAATCTTTGATGACGATGGaaaa                                 <  1:3159771/39‑1 (MQ=255)
                                                                 aTGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGAtt            <  1:3136071/60‑1 (MQ=255)
                                                                      |                                                                
AACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATTGGCTACCTGG  >  W3110S.gb/651828‑651962

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: