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

CCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTG  >  W3110S.gb/1093621‑1093750
                                                                 |                                                                
ccAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTg                                                              >  1:1243847/1‑70 (MQ=255)
   aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGa                                                          <  1:260456/71‑1 (MQ=255)
     aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                         >  1:908183/1‑70 (MQ=255)
     aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                         >  1:192139/1‑70 (MQ=255)
           ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                  <  1:289283/71‑1 (MQ=255)
           ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                  <  1:3709686/71‑1 (MQ=255)
                aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTATTATTTCTGGCATTGGTGAATAAAGTTAGTGa                                             <  1:4378343/71‑1 (MQ=255)
                aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                             <  1:2772180/71‑1 (MQ=255)
                aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                             <  1:162722/71‑1 (MQ=255)
                aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                             <  1:1358054/71‑1 (MQ=255)
                   tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          >  1:3145062/1‑71 (MQ=255)
                   tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                           >  1:5014498/1‑70 (MQ=255)
                     ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:285133/67‑1 (MQ=255)
                     ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:3686390/67‑1 (MQ=255)
                     ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:2788822/67‑1 (MQ=255)
                     ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:2462517/67‑1 (MQ=255)
                     ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:2480796/67‑1 (MQ=255)
                     ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCCGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:5127494/67‑1 (MQ=255)
                       attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                       >  1:4641271/1‑70 (MQ=255)
                        tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTTGTGAATAAAGTTAGTGAAAACACaa                                     <  1:2696927/71‑1 (MQ=255)
                        tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:1284947/71‑1 (MQ=255)
                        tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:807149/71‑1 (MQ=255)
                         ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:1332002/70‑1 (MQ=255)
                           gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   <  1:1172957/70‑1 (MQ=255)
                            tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   <  1:1241244/69‑1 (MQ=255)
                                 tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:221281/62‑1 (MQ=255)
                                 tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                     <  1:2467625/62‑1 (MQ=255)
                                   aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   >  1:2837511/1‑62 (MQ=255)
                                   aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   >  1:994363/1‑62 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATaaa                                                     >  1:1299986/1‑35 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:2467101/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:1019049/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:574983/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:5025492/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:1397905/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:3746258/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:2425490/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:1988340/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                    >  1:3045224/1‑36 (MQ=255)
                                            tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                  >  1:1728669/1‑70 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            >  1:2635292/1‑40 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            >  1:5205156/1‑40 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            >  1:799045/1‑40 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            >  1:4832733/1‑40 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            >  1:4102881/1‑40 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            >  1:360758/1‑40 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            >  1:1173831/1‑40 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaa             <  1:3683057/71‑1 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaa             <  1:3448316/71‑1 (MQ=255)
                                                cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaa             <  1:4726385/71‑1 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:3983421/36‑1 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                            <  1:3115942/36‑1 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                   >  1:3371203/1‑45 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaaa          <  1:3272897/70‑1 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt         >  1:1404196/1‑71 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt         >  1:5329926/1‑71 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt         >  1:1911536/1‑71 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt         >  1:4511538/1‑71 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt         >  1:528044/1‑71 (MQ=255)
                                                    ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt         >  1:1297751/1‑71 (MQ=255)
                                                      aaTATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          >  1:5477546/1‑36 (MQ=255)
                                                      aaTATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          >  1:3281577/1‑36 (MQ=255)
                                                      aaTATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          >  1:4129913/1‑36 (MQ=255)
                                                      aaTATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                          >  1:3639860/1‑36 (MQ=255)
                                                      aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATg        >  1:4424576/1‑70 (MQ=255)
                                                          tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt   >  1:2286916/1‑71 (MQ=255)
                                                              tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg  <  1:1281111/68‑1 (MQ=255)
                                                              tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg  <  1:1265460/68‑1 (MQ=255)
                                                              tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg  <  1:1827514/68‑1 (MQ=255)
                                                              tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtg  <  1:3143721/68‑1 (MQ=255)
                                                                 |                                                                
CCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTG  >  W3110S.gb/1093621‑1093750

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: