Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 651,883 G→T 10.8% M142I (ATG→ATT 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,8830GT10.8% 166.2 / 4.4 55M142I (ATG→ATTcitAsensory histidine kinase in two‑component regulatory system with citB
Reads supporting (aligned to +/- strand):  ref base G (21/28);  new base T (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 = 1.00e+00

TGCAGTTCACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATT  >  W3110S.gb/651816‑651952
                                                                   |                                                                     
tGCAGTTCACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATgg                                                                      <  1:1697490/69‑1 (MQ=255)
tGCAGTTCACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATgg                                                                      <  1:1560769/69‑1 (MQ=255)
  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)
                               aGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATg                                                                       >  1:3468288/1‑37 (MQ=255)
                                   gagaGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                 >  1:3591208/1‑71 (MQ=255)
                                   gagaGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                 <  1:1646040/71‑1 (MQ=255)
                                   gagaGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTtgatga                                 >  1:1359440/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:1464835/71‑1 (MQ=255)
                                           cTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGaa                         <  1:1530816/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:2429994/51‑1 (MQ=255)
                                                    cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                    <  1:613953/51‑1 (MQ=255)
                                                    cGGGAAAGGGTCAATTGGCATGGCGATGCGTGCCAAAACGCCAATCTTtga                                    <  1:474005/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)
                                                                   |                                                                     
TGCAGTTCACCAAACAGGGCGCGCTGGAGAAAGGGGAGAGCTACTTCATTACCGGGAAAGGGTCAATGGGGATGGCGATGCGCGCCAAAACGCCAATCTTTGATGACGATGGAAAAGTCATCGGCGTGGTGTCGATT  >  W3110S.gb/651816‑651952

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