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% 209.2 / NA 61V130A (GTA→GCA) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (42/19);  total (42/19)
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
                                                                   |                                                               
gtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                               <  1:2067486/70‑1 (MQ=255)
 tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:288240/70‑1 (MQ=255)
   cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:30090/1‑71 (MQ=255)
      aaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                        <  1:2417840/71‑1 (MQ=255)
             ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  >  1:1990085/1‑70 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:744155/1‑59 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:2959400/1‑59 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                <  1:770178/70‑1 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                <  1:1674554/70‑1 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                               >  1:346798/1‑71 (MQ=255)
               gCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                               >  1:1342872/1‑71 (MQ=255)
                cAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                <  1:2780339/69‑1 (MQ=255)
                   tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:2870289/64‑1 (MQ=255)
                    cTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           >  1:236973/1‑70 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGTCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:3999552/1‑71 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                >  1:3826561/1‑64 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                >  1:3901166/1‑64 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                >  1:2501319/1‑64 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAgagatag                                         >  1:3796498/1‑68 (MQ=255)
                     tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:2331258/1‑71 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                    >  1:28491/1‑58 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                    >  1:1116423/1‑58 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                    >  1:2681414/1‑58 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                    >  1:2562876/1‑58 (MQ=255)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                   >  1:1669496/1‑57 (MQ=255)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:1497947/1‑67 (MQ=255)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                      >  1:2287730/1‑70 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:3859195/71‑1 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:3448358/71‑1 (MQ=255)
                          tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:1989572/71‑1 (MQ=255)
                                tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                               >  1:2616975/1‑54 (MQ=255)
                                tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAg                                               >  1:1970386/1‑54 (MQ=255)
                                tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgc                              >  1:1654028/1‑71 (MQ=255)
                                  gTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgcgc                            <  1:2597595/71‑1 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         <  1:3127030/53‑1 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         <  1:2908234/53‑1 (MQ=255)
                                           gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:1850976/47‑1 (MQ=255)
                                           gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:1262260/47‑1 (MQ=255)
                                           gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:3758395/47‑1 (MQ=255)
                                           gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:1471254/47‑1 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                 >  1:3752892/1‑46 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                 >  1:2982879/1‑46 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                 >  1:443721/1‑46 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                 >  1:663634/1‑46 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAg                               >  1:1098069/1‑48 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAg                               >  1:1626701/1‑48 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAg                               >  1:4050074/1‑48 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:2969586/1‑71 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:2893746/1‑71 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:3985157/1‑71 (MQ=255)
                                                       tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                  >  1:2590691/1‑44 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                 >  1:3524406/1‑44 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                 >  1:1349350/1‑44 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATg       <  1:3465557/70‑1 (MQ=255)
                                                         atatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATg       <  1:3163412/69‑1 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAg                               >  1:2488081/1‑42 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt  >  1:2585929/1‑71 (MQ=255)
                                                              tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                     >  1:3434583/1‑50 (MQ=255)
                                                              tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                     >  1:2710508/1‑50 (MQ=255)
                                                              tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                     >  1:270927/1‑50 (MQ=255)
                                                              tCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                     >  1:1171371/1‑50 (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: