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

GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT  >  W3110S.gb/1093619‑1093748
                                                                   |                                                              
gtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                              <  1:3012713/70‑1 (MQ=255)
gtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                              <  1:2619467/70‑1 (MQ=255)
              tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                 >  1:155586/1‑69 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                               <  1:1111352/70‑1 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                            >  1:3196177/1‑51 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                          <  1:3153410/67‑1 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                    <  1:1312542/70‑1 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                    >  1:2172680/1‑70 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                   <  1:132734/71‑1 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                   <  1:2109097/71‑1 (MQ=255)
                              tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                               <  1:1748552/71‑1 (MQ=255)
                                tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAaca                                      >  1:769258/1‑62 (MQ=255)
                                 cGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                               >  1:3241372/1‑68 (MQ=255)
                                  gTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgcgc                           <  1:126420/71‑1 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCaa                         <  1:1659599/70‑1 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                        >  1:379184/1‑71 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                        >  1:3827418/1‑71 (MQ=255)
                                      gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                       <  1:1510764/71‑1 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                              >  1:256392/1‑44 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                              >  1:142582/1‑44 (MQ=255)
                                              tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                >  1:161011/1‑70 (MQ=255)
                                               gtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTcc               >  1:3352443/1‑70 (MQ=255)
                                               gtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTcc               >  1:1332440/1‑70 (MQ=255)
                                                 gCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaa             >  1:2622021/1‑70 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                 >  1:3308253/1‑45 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt       <  1:2718174/71‑1 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt       >  1:3088829/1‑71 (MQ=255)
                                                       tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                   >  1:973648/1‑42 (MQ=255)
                                                       tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt       >  1:170096/1‑70 (MQ=255)
                                                       tAATATTACTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                 >  1:780695/1‑44 (MQ=38)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAGACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt  <  1:1769214/70‑1 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGAGCAATTTATTTTCCAAAAAAATGACtt  <  1:1446549/70‑1 (MQ=255)
                                                                   |                                                              
GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTT  >  W3110S.gb/1093619‑1093748

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: