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% 126.3 / NA 39V130A (GTA→GCA) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/30);  total (9/30)
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
                                                                  |                                                                 
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:2583760/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:454969/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:5099/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:1428817/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:643328/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:1823434/70‑1 (MQ=255)
    aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCGTTGGTg                                                            <  1:4821/70‑1 (MQ=255)
     aaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                            <  1:580838/69‑1 (MQ=255)
              gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                  <  1:256485/70‑1 (MQ=255)
              gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                  <  1:1925934/70‑1 (MQ=255)
                 aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:2599417/53‑1 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                            <  1:1295078/68‑1 (MQ=255)
                        attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                        >  1:2201135/1‑70 (MQ=255)
                        attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                        >  1:255564/1‑70 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             <  1:1015170/64‑1 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:1204849/71‑1 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:421976/71‑1 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:2485999/71‑1 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:868634/71‑1 (MQ=255)
                          ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:434140/70‑1 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                                 >  1:127464/1‑57 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                                 >  1:1269453/1‑57 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:619345/70‑1 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:2243726/70‑1 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:2034723/70‑1 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:1923181/70‑1 (MQ=255)
                             tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                  <  1:2242092/71‑1 (MQ=255)
                             tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                  <  1:1905603/71‑1 (MQ=255)
                              ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                  <  1:25612/70‑1 (MQ=255)
                              ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                  <  1:1094276/70‑1 (MQ=255)
                                                 cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaa               >  1:1928501/1‑70 (MQ=255)
                                                 cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaa               >  1:1775110/1‑70 (MQ=255)
                                                  atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa             <  1:559205/71‑1 (MQ=255)
                                                  atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa             <  1:865211/71‑1 (MQ=255)
                                                      tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt          >  1:1754719/1‑70 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                          <  1:383656/51‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                          <  1:2142549/51‑1 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt    >  1:655220/1‑71 (MQ=255)
                                                             tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt  >  1:1729244/1‑71 (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: