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

AAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTT  >  W3110S.gb/1093624‑1093749
                                                              |                                                               
aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                          >  1:2035765/1‑70 (MQ=255)
aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                          >  1:1980992/1‑70 (MQ=255)
   tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                      <  1:1619491/71‑1 (MQ=255)
     aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                            >  1:2098006/1‑63 (MQ=255)
       aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  >  1:676915/1‑71 (MQ=255)
       aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:1712384/71‑1 (MQ=255)
          gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                               >  1:177203/1‑71 (MQ=255)
           cAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                      <  1:1303084/63‑1 (MQ=255)
               cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                 <  1:990823/64‑1 (MQ=255)
                tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:1609458/50‑1 (MQ=255)
                    attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                   >  1:1576249/1‑57 (MQ=255)
                     tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:874401/71‑1 (MQ=255)
                     tttGTTACGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:2049502/1‑66 (MQ=255)
                      ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:1876754/63‑1 (MQ=255)
                           tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAGACaaaa                                  >  1:871195/1‑67 (MQ=255)
                             gTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgcgc                            <  1:993497/71‑1 (MQ=255)
                               cAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:1743668/61‑1 (MQ=255)
                                aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                         >  1:466754/1‑71 (MQ=255)
                                   tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      <  1:366852/71‑1 (MQ=255)
                                    gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:296427/1‑51 (MQ=255)
                                    gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:1308679/1‑51 (MQ=255)
                                     aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                    <  1:1296873/71‑1 (MQ=255)
                                     aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                    <  1:1304524/71‑1 (MQ=255)
                                      gTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                      <  1:1975429/36‑1 (MQ=255)
                                       tttGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                     <  1:1809081/36‑1 (MQ=255)
                                          gtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAAAGCGCAATTTATTTTcc                >  1:374121/1‑70 (MQ=255)
                                               tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:565001/38‑1 (MQ=255)
                                                 ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                 >  1:1940006/1‑46 (MQ=255)
                                                 ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:654572/1‑71 (MQ=255)
                                                  tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAaca                                       >  1:779731/1‑39 (MQ=255)
                                                   aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  >  1:1043635/1‑43 (MQ=255)
                                                   aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGa      >  1:1768800/1‑71 (MQ=255)
                                                       tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      >  1:1628392/1‑51 (MQ=255)
                                                       tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      >  1:977280/1‑51 (MQ=255)
                                                       tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt  >  1:765906/1‑71 (MQ=255)
                                                              |                                                               
AAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTT  >  W3110S.gb/1093624‑1093749

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: