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% 125.7 / NA 37V130A (GTA→GCA) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/15);  total (22/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTG  >  W3110S.gb/1093622‑1093750
                                                                |                                                                
cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                              >  1:3623939/1‑69 (MQ=255)
cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                              <  1:1819847/69‑1 (MQ=255)
            gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                           >  1:1561662/1‑60 (MQ=255)
            gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                  >  1:198933/1‑69 (MQ=255)
                 cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                   >  1:2115903/1‑63 (MQ=255)
                 cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTg                                              >  1:2008984/1‑68 (MQ=255)
                  tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            >  1:506471/1‑69 (MQ=255)
                      attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAac                                         >  1:3599689/1‑68 (MQ=255)
                      attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAac                                         >  1:1899606/1‑68 (MQ=255)
                         tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:2041386/69‑1 (MQ=255)
                         tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:1735333/69‑1 (MQ=255)
                                      gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                      >  1:285418/1‑39 (MQ=255)
                                      gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                      >  1:2396962/1‑39 (MQ=255)
                                      gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                      >  1:564595/1‑39 (MQ=255)
                                      gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattt                        >  1:392303/1‑69 (MQ=255)
                                              gCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                  <  1:1286146/35‑1 (MQ=255)
                                                  atTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa            <  1:1479014/69‑1 (MQ=255)
                                                  atTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa            <  1:1029143/69‑1 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          <  1:1322993/38‑1 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          <  1:1310312/38‑1 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     >  1:570551/1‑43 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     >  1:3557202/1‑43 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     >  1:2995078/1‑43 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     >  1:2128800/1‑43 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     >  1:2558648/1‑43 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa           >  1:1675379/1‑69 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa           >  1:491253/1‑69 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa           >  1:1310417/1‑69 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa           >  1:1302294/1‑69 (MQ=255)
                                                        atTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:3283674/38‑1 (MQ=255)
                                                        atTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:1938347/38‑1 (MQ=255)
                                                        atTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:1412935/38‑1 (MQ=255)
                                                         tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGAGCAATTTATTTTCCAAAAAAATGACt     >  1:540891/1‑69 (MQ=255)
                                                           tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt   <  1:2766753/69‑1 (MQ=255)
                                                           tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt    <  1:347889/68‑1 (MQ=255)
                                                            cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt   <  1:896943/68‑1 (MQ=255)
                                                             tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg  <  1:2527025/68‑1 (MQ=255)
                                                                |                                                                
CAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTG  >  W3110S.gb/1093622‑1093750

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: