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

TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093620‑1093751
                                                                  |                                                                 
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                               <  1:1624389/71‑1 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                       <  1:984040/70‑1 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                       <  1:1289854/70‑1 (MQ=255)
           aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    <  1:1577305/71‑1 (MQ=255)
                  tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                       <  1:888815/61‑1 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATAAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                            >  1:184945/1‑70 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                       >  1:51226/1‑70 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:316606/71‑1 (MQ=255)
                          ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:1343187/70‑1 (MQ=255)
                               tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAg                                 <  1:1136009/70‑1 (MQ=255)
                                   cAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgc                                >  1:299144/1‑67 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                    <  1:160009/71‑1 (MQ=255)
                                                 cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaa               >  1:1109957/1‑70 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:830917/47‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa            <  1:1709864/71‑1 (MQ=255)
                                                    atTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                       <  1:600109/59‑1 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt          <  1:5684/71‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:1744770/39‑1 (MQ=255)
                                                             tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt  >  1:52644/1‑71 (MQ=255)
                                                             tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg   <  1:1742617/70‑1 (MQ=255)
                                                             tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg   >  1:1164027/1‑70 (MQ=255)
                                                                  |                                                                 
TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093620‑1093751

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: