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

GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTT  >  W3110S.gb/1093619‑1093749
                                                                   |                                                               
gtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              >  1:2483163/1‑71 (MQ=255)
   cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           <  1:3424205/71‑1 (MQ=255)
   cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:2722334/1‑71 (MQ=255)
    aaaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:3450593/1‑70 (MQ=255)
     aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                          >  1:2538520/1‑70 (MQ=255)
        tGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATGGGTGAATa                                                      <  1:1613567/71‑1 (MQ=255)
        caacacTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAAt                                                       <  1:2081620/68‑1 (MQ=255)
             ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTgg                                                            <  1:4360661/60‑1 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                               >  1:2096845/1‑71 (MQ=255)
                   tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                     <  1:4304114/61‑1 (MQ=255)
                   tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                     <  1:2917212/61‑1 (MQ=255)
                    cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           >  1:1052751/1‑70 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                >  1:837760/1‑64 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                >  1:2986598/1‑64 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:4028603/1‑71 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGAATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:2984452/1‑71 (MQ=255)
                      attattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:1974735/1‑70 (MQ=255)
                      attattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:2854497/1‑70 (MQ=255)
                      attattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:4120012/1‑70 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                   >  1:4102920/1‑59 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:3551814/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:530482/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:3757258/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:766461/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:2738574/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTATGGCAt                                                               <  1:2782030/47‑1 (MQ=39)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                      >  1:808128/1‑70 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                     >  1:2853538/1‑70 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:1178930/71‑1 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:3578864/71‑1 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:2381414/71‑1 (MQ=255)
                             gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  <  1:2345517/70‑1 (MQ=255)
                              tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:4275906/67‑1 (MQ=255)
                              tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                 <  1:3914869/70‑1 (MQ=255)
                                 cGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgcg                             >  1:344015/1‑71 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGTATAAAGTTAGTGAAAACACAAAACAGCGCAat                         <  1:738047/71‑1 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                    >  1:3584932/1‑44 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                         <  1:389805/71‑1 (MQ=255)
                                        tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      >  1:3350242/1‑71 (MQ=255)
                                        tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      >  1:2685074/1‑71 (MQ=255)
                                         gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      >  1:3353940/1‑70 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                    <  1:3188804/71‑1 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                    <  1:2964124/71‑1 (MQ=255)
                                           gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttattt                   <  1:4051151/71‑1 (MQ=255)
                                             ttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                  >  1:4286204/1‑70 (MQ=255)
                                              tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                >  1:2336768/1‑39 (MQ=255)
                                              tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                 >  1:3279693/1‑70 (MQ=255)
                                                    tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa          <  1:1535258/71‑1 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAACCAGCGCAATTTATTTTCCAAAAAAAt        >  1:1022998/1‑71 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:3785293/1‑71 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:4424667/1‑71 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  >  1:342814/1‑43 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt  >  1:4363302/1‑71 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt  <  1:2931917/71‑1 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt   >  1:201148/1‑70 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt   >  1:2628785/1‑70 (MQ=255)
                                                              tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                >  1:1998327/1‑39 (MQ=255)
                                                                   |                                                               
GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTT  >  W3110S.gb/1093619‑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: