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

GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093619‑1093751
                                                                   |                                                                 
gtCCAAAATGTCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtc                                                                >  1:2109061/1‑70 (MQ=255)
gtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                                 <  1:1828774/70‑1 (MQ=255)
gtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                                 <  1:1555668/70‑1 (MQ=255)
 tCCAAAATGACTTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:1082694/70‑1 (MQ=255)
 tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:1157548/70‑1 (MQ=255)
 tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:214970/70‑1 (MQ=255)
 tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:4436644/70‑1 (MQ=255)
 tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                                 <  1:4492229/69‑1 (MQ=255)
 tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                               <  1:1512400/71‑1 (MQ=255)
  ccAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:3955572/69‑1 (MQ=255)
   cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                             >  1:186221/1‑71 (MQ=255)
    aaaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                             >  1:3456025/1‑70 (MQ=255)
     aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                              >  1:1384176/1‑68 (MQ=255)
     aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                            >  1:2720040/1‑70 (MQ=255)
      aaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                          >  1:4572790/1‑71 (MQ=255)
      aaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                          >  1:785032/1‑71 (MQ=255)
       aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                          >  1:490160/1‑70 (MQ=255)
         gACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                <  1:1062820/62‑1 (MQ=255)
            aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    <  1:4681076/71‑1 (MQ=255)
            aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    >  1:432342/1‑71 (MQ=255)
             ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    <  1:1957738/70‑1 (MQ=255)
             ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    <  1:2026885/70‑1 (MQ=255)
             ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    >  1:1630700/1‑70 (MQ=255)
              tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                      >  1:1194083/1‑67 (MQ=255)
              tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    >  1:756616/1‑69 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    >  1:4587461/1‑68 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    >  1:2133838/1‑68 (MQ=255)
                cAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                <  1:582496/71‑1 (MQ=255)
                 aaTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                               >  1:670671/1‑55 (MQ=255)
                      attattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:2617084/1‑70 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                <  1:2780294/64‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAaca                                         >  1:3131113/1‑71 (MQ=255)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                        >  1:220500/1‑70 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGTATAAAg                                                     >  1:882299/1‑56 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                     >  1:1801323/1‑56 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:2392552/71‑1 (MQ=255)
                           ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                              >  1:2172004/1‑46 (MQ=255)
                             gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:4096942/70‑1 (MQ=255)
                             gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:2195980/70‑1 (MQ=255)
                             gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:2080430/70‑1 (MQ=255)
                              tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                   <  1:2831903/70‑1 (MQ=255)
                              tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                   <  1:2964009/70‑1 (MQ=255)
                              tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                   <  1:1454516/70‑1 (MQ=255)
                               ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                             >  1:2916711/1‑43 (MQ=255)
                                tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgc                                >  1:1312516/1‑71 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:2246783/60‑1 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                  >  1:3985133/1‑46 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                  >  1:1105789/1‑46 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    >  1:1964838/1‑60 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                          >  1:1140558/1‑70 (MQ=255)
                                         gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                  >  1:3972044/1‑60 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                      <  1:1244442/71‑1 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                      <  1:2790265/71‑1 (MQ=255)
                                           gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                           >  1:1873338/1‑49 (MQ=255)
                                                tgCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                    <  1:601899/35‑1 (MQ=255)
                                                 gCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaa                >  1:2563945/1‑70 (MQ=255)
                                                   atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa             <  1:2744255/71‑1 (MQ=255)
                                                    tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:3628950/47‑1 (MQ=255)
                                                    tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCa                 <  1:1928759/66‑1 (MQ=255)
                                                     atTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaaa           >  1:4124083/1‑71 (MQ=255)
                                                     atTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaa              <  1:3186572/68‑1 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                             <  1:264796/36‑1 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    >  1:1614028/1‑45 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                        <  1:2815068/57‑1 (MQ=255)
                                                       tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt          >  1:919611/1‑70 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                  >  1:2553873/1‑45 (MQ=255)
                                                         atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAc       >  1:4450221/1‑71 (MQ=255)
                                                         atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGAc       >  1:4017673/1‑71 (MQ=255)
                                                          tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                    <  1:3913086/41‑1 (MQ=255)
                                                          tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                      <  1:3088223/39‑1 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                    >  1:551725/1‑55 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt    >  1:23704/1‑71 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt    >  1:1966297/1‑71 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt     >  1:3977268/1‑70 (MQ=255)
                                                              tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt  >  1:3081631/1‑71 (MQ=255)
                                                              tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt  >  1:69049/1‑71 (MQ=255)
                                                               cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                       <  1:309898/49‑1 (MQ=255)
                                                                   |                                                                 
GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGT  >  W3110S.gb/1093619‑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: