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

AAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093623‑1093751
                                                               |                                                                 
aaaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                             >  1:4215498/1‑70 (MQ=255)
    tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                        <  1:2180340/71‑1 (MQ=255)
      aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                       <  1:2371130/70‑1 (MQ=255)
        aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    <  1:1474528/71‑1 (MQ=255)
        aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    >  1:167944/1‑71 (MQ=255)
          tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                        <  1:4256093/65‑1 (MQ=255)
          tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    >  1:3223526/1‑69 (MQ=255)
           gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                                 >  1:3193455/1‑71 (MQ=255)
                cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             >  1:1451597/1‑70 (MQ=255)
                 tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:276856/1‑71 (MQ=255)
                 tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGAGAATAAAg                                                     >  1:1299477/1‑61 (MQ=255)
                   ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             <  1:835268/67‑1 (MQ=255)
                   ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             <  1:3537210/67‑1 (MQ=255)
                   ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAac                                          <  1:845942/70‑1 (MQ=255)
                   ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAac                                          <  1:3479304/70‑1 (MQ=255)
                    tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           <  1:2195498/68‑1 (MQ=255)
                    tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                              <  1:137383/65‑1 (MQ=255)
                     attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                        <  1:4215110/54‑1 (MQ=255)
                     attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                        >  1:1436525/1‑70 (MQ=255)
                     attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                        <  1:847062/70‑1 (MQ=255)
                      tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                              >  1:2034276/1‑47 (MQ=255)
                      tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                <  1:1401673/61‑1 (MQ=255)
                      tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTg                                               <  1:2102597/62‑1 (MQ=255)
                      tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                       >  1:349543/1‑70 (MQ=255)
                        tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                                 >  1:3751178/1‑42 (MQ=255)
                        tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                                 >  1:2324696/1‑42 (MQ=255)
                         gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:2500604/70‑1 (MQ=255)
                         gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:609502/70‑1 (MQ=255)
                         gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:3625403/68‑1 (MQ=255)
                          tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                               <  1:1873151/42‑1 (MQ=255)
                            tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                        <  1:1781415/47‑1 (MQ=255)
                            tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    >  1:3727696/1‑67 (MQ=255)
                            tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:4239473/65‑1 (MQ=255)
                                 aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                            >  1:2209284/1‑54 (MQ=255)
                                    tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                        <  1:3676065/71‑1 (MQ=255)
                                    tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                          >  1:2738731/1‑69 (MQ=255)
                                      aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                      <  1:2760487/71‑1 (MQ=255)
                                      aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                      <  1:2059005/71‑1 (MQ=255)
                                        tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                              >  1:2903014/1‑45 (MQ=255)
                                        tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                    <  1:3552540/71‑1 (MQ=255)
                                            tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:2479123/51‑1 (MQ=255)
                                            tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                      <  1:2728616/65‑1 (MQ=255)
                                             gCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                  <  1:25225/52‑1 (MQ=255)
                                               atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                        >  1:227824/1‑60 (MQ=255)
                                                  ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt          <  1:2621333/71‑1 (MQ=255)
                                                   tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    >  1:3426181/1‑44 (MQ=255)
                                                   tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaaa           >  1:895040/1‑69 (MQ=255)
                                                    aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGa        >  1:790992/1‑71 (MQ=255)
                                                        tttCTTGCATTGGTGAATAAAGTTAGTGAAAacac                                        <  1:1003401/35‑1 (MQ=38)
                                                        tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt    >  1:924845/1‑71 (MQ=255)
                                                          tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt  >  1:374100/1‑71 (MQ=255)
                                                           cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgcgc                              <  1:1382726/42‑1 (MQ=255)
                                                               |                                                                 
AAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093623‑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: