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

TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTT  >  W3110S.gb/1093620‑1093749
                                                                  |                                                               
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:4310527/70‑1 (MQ=255)
tCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:2173409/70‑1 (MQ=255)
  cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           <  1:1677539/71‑1 (MQ=255)
  cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:23560/1‑71 (MQ=255)
  cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:4736751/1‑71 (MQ=255)
  cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                           >  1:223751/1‑71 (MQ=255)
     aaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                        <  1:5115249/71‑1 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                     <  1:2357202/70‑1 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                     <  1:666272/70‑1 (MQ=255)
         aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                     <  1:4468597/70‑1 (MQ=255)
            ttGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                 <  1:2957368/71‑1 (MQ=255)
                 aTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                            <  1:3124696/71‑1 (MQ=255)
                  tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:4970478/64‑1 (MQ=255)
                  tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                            >  1:888435/1‑70 (MQ=255)
                    tattatTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTa                                                >  1:2131433/1‑64 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:4665612/1‑69 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:3435988/1‑69 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:4763992/67‑1 (MQ=255)
                      ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAaca                                       >  1:1879974/1‑71 (MQ=255)
                       tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:2780070/47‑1 (MQ=255)
                       tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:27226/47‑1 (MQ=255)
                       tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:4515117/47‑1 (MQ=255)
                       tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                              <  1:4933340/47‑1 (MQ=255)
                        attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                      >  1:3614264/1‑70 (MQ=255)
                        attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                      >  1:114476/1‑70 (MQ=255)
                         tttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                      <  1:3280166/69‑1 (MQ=255)
                          ttGTTTCGTCAGATGAGTTTGTTCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         <  1:4726673/65‑1 (MQ=255)
                          ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:4999601/56‑1 (MQ=255)
                          ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:141903/56‑1 (MQ=255)
                          ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:2637484/56‑1 (MQ=255)
                          ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:2698425/56‑1 (MQ=255)
                           tGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:653624/1‑64 (MQ=255)
                             tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:980237/67‑1 (MQ=255)
                               tCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgc                              >  1:2332014/1‑71 (MQ=255)
                                  tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:4155444/55‑1 (MQ=255)
                                  tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:2307838/55‑1 (MQ=255)
                                   cAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCaa                          >  1:431377/1‑71 (MQ=255)
                                    aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                         >  1:4676274/1‑71 (MQ=255)
                                    aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                         >  1:4623049/1‑71 (MQ=255)
                                    aGATGAGTTTGTGCATATTAATATTTCTGGCATTGATGAATAAAGTTAGTGAAAACACAAAACAGCGCAat                         >  1:4670459/1‑71 (MQ=255)
                                       tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                >  1:1587500/1‑61 (MQ=255)
                                       tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                      >  1:5072334/1‑71 (MQ=255)
                                        gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                     >  1:3122616/1‑71 (MQ=255)
                                        gAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                     >  1:4489986/1‑71 (MQ=255)
                                         aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                         >  1:3462954/1‑50 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:1836601/39‑1 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                  <  1:2379572/39‑1 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaa                                    <  1:3945029/53‑1 (MQ=255)
                                           tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                  <  1:4585157/71‑1 (MQ=255)
                                             tgtgCATATTAATATTTCTGGCATTGGTGAATAAAg                                                   >  1:3011669/1‑36 (MQ=255)
                                             tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                 >  1:3378443/1‑70 (MQ=255)
                                                 cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                    <  1:5108977/63‑1 (MQ=255)
                                                 cATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaa             >  1:2381816/1‑70 (MQ=255)
                                                  atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa           <  1:633890/71‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:2540911/38‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaa                                   <  1:3741561/46‑1 (MQ=255)
                                                   tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa          <  1:2556718/71‑1 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                           <  1:2115100/36‑1 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:1000870/1‑71 (MQ=255)
                                                     ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:1761069/1‑71 (MQ=255)
                                                      tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                 <  1:241528/61‑1 (MQ=255)
                                                      tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                 <  1:2543642/61‑1 (MQ=255)
                                                      tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:1891651/1‑70 (MQ=255)
                                                      tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt        >  1:1885662/1‑70 (MQ=255)
                                                       aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                <  1:2480497/45‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                <  1:2432175/43‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACt    <  1:4513221/71‑1 (MQ=255)
                                                         tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACt    <  1:676343/71‑1 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAATACACaa                                    <  1:662630/37‑1 (MQ=38)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt  <  1:1907306/71‑1 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt  <  1:2382351/71‑1 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt  <  1:2713272/71‑1 (MQ=255)
                                                           tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt   <  1:2884294/70‑1 (MQ=255)
                                                              cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa          <  1:1210940/60‑1 (MQ=255)
                                                                  |                                                               
TCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTT  >  W3110S.gb/1093620‑1093749

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: