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

TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAC  >  W3110S.gb/1093620‑1093746
                                                                  |                                                            
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                           <  1:1578766/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                          <  1:1361119/71‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                          <  1:806869/71‑1 (MQ=255)
      aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAAt                                                    <  1:938289/71‑1 (MQ=255)
           aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                               <  1:575126/71‑1 (MQ=255)
            ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                        <  1:1333517/61‑1 (MQ=255)
              gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGa                                                      >  1:569312/1‑61 (MQ=255)
              gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                             <  1:3167563/70‑1 (MQ=255)
              gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                             <  1:3251940/70‑1 (MQ=255)
              gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                            >  1:2892464/1‑71 (MQ=255)
                 aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTg                                          >  1:853326/1‑70 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                      >  1:1973851/1‑71 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                      >  1:216407/1‑71 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                      >  1:2053100/1‑71 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                      >  1:879043/1‑71 (MQ=255)
                        attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                   >  1:3287780/1‑70 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                >  1:2916645/1‑56 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                               >  1:169868/1‑70 (MQ=255)
                             tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                        <  1:3157782/60‑1 (MQ=255)
                                    aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                      >  1:1156982/1‑71 (MQ=255)
                                     gATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                     <  1:111902/71‑1 (MQ=255)
                                       tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                               >  1:3050821/1‑59 (MQ=255)
                                          gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                               <  1:2166112/40‑1 (MQ=255)
                                          gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                             >  1:3231772/1‑42 (MQ=255)
                                                  atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                     <  1:2167487/58‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                        <  1:181708/38‑1 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaaa      >  1:1040974/1‑70 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt     <  1:782665/71‑1 (MQ=255)
                                                     ttAATATTTCTGGCAGTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt     >  1:2230380/1‑71 (MQ=255)
                                                      tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt     >  1:3193720/1‑70 (MQ=255)
                                                        atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGa   <  1:3385055/70‑1 (MQ=255)
                                                        atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAc  <  1:1456244/71‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                   <  1:1838406/37‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                  <  1:30074/38‑1 (MQ=255)
                                                             tCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                               <  1:949222/37‑1 (MQ=255)
                                                                  |                                                            
TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAC  >  W3110S.gb/1093620‑1093746

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: