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

CAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093622‑1093751
                                                                |                                                                 
cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                              <  1:426269/70‑1 (MQ=255)
cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                              <  1:2861712/70‑1 (MQ=255)
cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                              <  1:2842373/70‑1 (MQ=255)
cAAAATGACATTGCAAACTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                             >  1:1862497/1‑71 (MQ=255)
   aaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                          >  1:4102121/1‑71 (MQ=255)
     tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                        <  1:804181/71‑1 (MQ=255)
     tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                        <  1:493496/71‑1 (MQ=255)
     tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                        <  1:4726273/71‑1 (MQ=255)
       aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                       <  1:2659738/70‑1 (MQ=255)
            gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                  <  1:3492026/70‑1 (MQ=255)
            gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                                 >  1:2722763/1‑71 (MQ=255)
                 cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                >  1:1635364/1‑67 (MQ=255)
                  tattatTTGTTTCGTCGGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:4561399/1‑71 (MQ=255)
                  tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                              >  1:3291607/1‑52 (MQ=255)
                  tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                              >  1:1578734/1‑52 (MQ=255)
                  tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                  >  1:4347819/1‑64 (MQ=255)
                  tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:1021410/1‑71 (MQ=255)
                    ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             <  1:4787632/67‑1 (MQ=255)
                     tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    <  1:3834906/59‑1 (MQ=255)
                      attTGTTTCGTCATATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:2205959/1‑67 (MQ=255)
                      attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:2180672/1‑67 (MQ=255)
                      attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:2459773/1‑67 (MQ=255)
                      attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                        >  1:4487415/1‑70 (MQ=255)
                      attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                        >  1:1963928/1‑70 (MQ=255)
                       tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                                 >  1:2216309/1‑60 (MQ=255)
                       tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGAGAAAAcaca                                       >  1:4280409/1‑70 (MQ=255)
                        ttGTTTCGTCAGATGAGTTTGTGCATTTTAATATTTCTGGCATTGGTGAATTAAGTTa                                                  <  1:10259/58‑1 (MQ=255)
                           tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGATAACACAAAAc                                   <  1:413881/70‑1 (MQ=255)
                              cGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgc                                <  1:4711314/70‑1 (MQ=255)
                                   gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                          >  1:320558/1‑71 (MQ=255)
                                     tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                        <  1:793046/71‑1 (MQ=255)
                                      gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                       >  1:1533869/1‑39 (MQ=255)
                                      gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                  >  1:2550287/1‑44 (MQ=255)
                                      gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                  >  1:1342174/1‑44 (MQ=255)
                                         tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                    <  1:4167963/71‑1 (MQ=255)
                                           tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                   >  1:4597525/1‑70 (MQ=255)
                                            gtgCATATTAATATTTCTGGCATTGGTTAATAAAGTTAGTGaa                                             <  1:2143232/43‑1 (MQ=38)
                                                atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa             <  1:3176725/71‑1 (MQ=255)
                                                atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa             <  1:1112625/71‑1 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaaa           <  1:4109743/70‑1 (MQ=255)
                                                   ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt          >  1:2898783/1‑71 (MQ=255)
                                                    tAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:82356/1‑37 (MQ=255)
                                                    tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt          >  1:4370961/1‑70 (MQ=255)
                                                    tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTACAAAAAAAt          >  1:2885057/1‑70 (MQ=255)
                                                     aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa            <  1:3118770/67‑1 (MQ=255)
                                                      atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAc       >  1:271533/1‑71 (MQ=255)
                                                       tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:4090243/39‑1 (MQ=255)
                                                           tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGTAATTTATTTTCCAAAAAAATGACTTtgt  >  1:4279690/1‑71 (MQ=255)
                                                                |                                                                 
CAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093622‑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: