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% 147.7 / NA 44V130A (GTA→GCA) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/27);  total (17/27)
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:12006/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:2781289/70‑1 (MQ=255)
  cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                             <  1:731945/71‑1 (MQ=255)
  cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                             >  1:3515161/1‑71 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                       <  1:2410148/70‑1 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                       <  1:2822395/70‑1 (MQ=255)
           aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    <  1:146403/71‑1 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:1294427/1‑71 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:2173700/1‑71 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                              <  1:2655411/68‑1 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             <  1:3313341/67‑1 (MQ=255)
                       tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             <  1:968822/66‑1 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                                 >  1:1834042/1‑60 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:304021/1‑66 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:500909/1‑66 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:1603572/1‑66 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:2588158/1‑66 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:1197048/71‑1 (MQ=255)
                            gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                             >  1:2820213/1‑45 (MQ=255)
                                      aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                            <  1:3494604/36‑1 (MQ=255)
                                      aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:3296386/58‑1 (MQ=255)
                                       tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                        <  1:1747749/71‑1 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                     <  1:562271/38‑1 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                        <  1:1606243/67‑1 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                    <  1:2197068/71‑1 (MQ=255)
                                             tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                   >  1:287859/1‑70 (MQ=255)
                                               tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                      <  1:2822461/65‑1 (MQ=255)
                                               tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                      <  1:196499/65‑1 (MQ=255)
                                               tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTcc                  >  1:2581026/1‑69 (MQ=255)
                                               tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCa                 <  1:762993/70‑1 (MQ=255)
                                               tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCa                 <  1:758783/70‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           <  1:2480189/40‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             <  1:1312942/38‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa            <  1:1385032/71‑1 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaa                                     >  1:1980000/1‑44 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa             >  1:3363527/1‑68 (MQ=255)
                                                       aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    >  1:1711654/1‑43 (MQ=255)
                                                        atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGa        <  1:2149172/70‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                   <  1:854673/42‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAc       <  1:3600683/70‑1 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                       >  1:3095558/1‑36 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt    <  1:973765/71‑1 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGAGAAAACACaaaa                                    >  1:3045324/1‑39 (MQ=38)
                                                             tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt  >  1:787030/1‑71 (MQ=255)
                                                                  |                                                                 
TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093620‑1093751

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