Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,093,686 T→C 100% V130A (GTA→GCA)  ycdT → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,093,6860TC100.0% 218.9 / NA 63V130A (GTA→GCA) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/32);  total (31/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTT  >  W3110S.gb/1093620‑1093749
                                                                  |                                                               
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:1950684/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                             <  1:2510470/71‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                             <  1:741072/71‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                             <  1:1726880/71‑1 (MQ=255)
  cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:4601169/1‑71 (MQ=255)
    aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                          >  1:3753135/1‑70 (MQ=255)
      aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGAGaa                                                        >  1:4584049/1‑70 (MQ=255)
       tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                      <  1:5139346/71‑1 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                      >  1:796920/1‑69 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                      >  1:2833247/1‑69 (MQ=255)
           aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:530876/71‑1 (MQ=255)
           aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:1309304/71‑1 (MQ=255)
             tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  >  1:628782/1‑69 (MQ=255)
              gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                             >  1:4357522/1‑57 (MQ=255)
              gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGa                                                         >  1:3862327/1‑61 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                    >  1:4473630/1‑60 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                    >  1:763020/1‑60 (MQ=255)
                     attattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                 >  1:3463127/1‑62 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAAt                                                       <  1:2323284/55‑1 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAAt                                                       <  1:2397619/55‑1 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                          <  1:1536856/68‑1 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                          <  1:2977973/68‑1 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:3862751/67‑1 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:2447046/67‑1 (MQ=255)
                        attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                    >  1:870793/1‑56 (MQ=255)
                        attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                      >  1:1768109/1‑70 (MQ=255)
                        attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                      >  1:1053628/1‑70 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:2990312/71‑1 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:655261/71‑1 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:2200525/71‑1 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:624532/71‑1 (MQ=255)
                          ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:4084969/70‑1 (MQ=255)
                           tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAaca                                       >  1:4103745/1‑66 (MQ=255)
                           tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaa                                   <  1:1251829/70‑1 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  <  1:4076243/70‑1 (MQ=255)
                               tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgc                              >  1:1826287/1‑71 (MQ=255)
                                 gTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgcg                             >  1:2792636/1‑70 (MQ=255)
                                 gTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgcg                             >  1:530356/1‑70 (MQ=255)
                                  tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCa                           <  1:2472883/71‑1 (MQ=255)
                                    aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCaa                          <  1:4899362/70‑1 (MQ=255)
                                    aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                         >  1:374780/1‑71 (MQ=255)
                                     gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                        <  1:1486743/71‑1 (MQ=255)
                                     gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                        <  1:4778497/71‑1 (MQ=255)
                                       tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      <  1:2636231/71‑1 (MQ=255)
                                        gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  >  1:2637746/1‑58 (MQ=255)
                                         aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCaa                          <  1:47800/65‑1 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                  <  1:4245052/71‑1 (MQ=255)
                                                 cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaa            <  1:4243368/71‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                 >  1:1244782/1‑64 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                 >  1:157155/1‑64 (MQ=255)
                                                    atTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      <  1:2315755/58‑1 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      >  1:4015788/1‑57 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:4111202/1‑71 (MQ=255)
                                                      tAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:2044667/1‑37 (MQ=255)
                                                      tAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:3312355/1‑37 (MQ=255)
                                                          atTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  <  1:4104535/40‑1 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  >  1:1478967/1‑39 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  >  1:1867934/1‑39 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  >  1:317141/1‑39 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGa      <  1:2194788/67‑1 (MQ=255)
                                                            ttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt  <  1:787920/70‑1 (MQ=255)
                                                              cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTcc                >  1:1573632/1‑54 (MQ=255)
                                                              cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTcc                >  1:4627482/1‑54 (MQ=255)
                                                                  |                                                               
TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTT  >  W3110S.gb/1093620‑1093749

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: