Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,094,325 A→G 9.4% N343S (AAT→AGT)  ycdT → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,094,3250AG 9.4% 135.0 / 5.4 53N343S (AAT→AGT) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (19/29);  new base G (3/2);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTT  >  W3110S.gb/1094255‑1094394
                                                                      |                                                                     
gACCACTTCAAATAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGaa                                                                       >  1:1020103/1‑71 (MQ=255)
   cACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAatat                                                                    <  1:685150/71‑1 (MQ=255)
       tCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCAtt                                                                >  1:968629/1‑71 (MQ=255)
       tCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCAtt                                                                >  1:1801468/1‑71 (MQ=255)
        cAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCa                                                                  >  1:2223569/1‑68 (MQ=255)
        cAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCa                                                                  >  1:2776653/1‑68 (MQ=255)
        cAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCa                                                                  >  1:308498/1‑68 (MQ=255)
        cAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCa                                                                  >  1:3666927/1‑68 (MQ=255)
        cAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTg                                                               >  1:3793978/1‑71 (MQ=255)
            aaaGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGa                                                            >  1:2035508/1‑70 (MQ=255)
            aaaGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGa                                                            >  1:3092937/1‑70 (MQ=255)
               gTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAg                                                        <  1:3302734/71‑1 (MQ=255)
                  aaCGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCa                                                      <  1:1705091/70‑1 (MQ=255)
                  aaCGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCa                                                      <  1:1906417/70‑1 (MQ=255)
                  aaCGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTAAATATCATTGGGAAAAGCa                                                      <  1:510124/70‑1 (MQ=255)
                    cGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCa                                                                  <  1:765936/56‑1 (MQ=255)
                    cGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCacac                                                   <  1:1651540/71‑1 (MQ=255)
                    cGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATAc                                                   <  1:1094717/71‑1 (MQ=255)
                        acCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAGTATCATTGGGAAAAGCATACg                                                  >  1:1042860/1‑68 (MQ=255)
                        acCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAGTATCATTGGGAAAAGCATACg                                                  >  1:1622244/1‑68 (MQ=255)
                         ccTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCa                                              <  1:3544606/71‑1 (MQ=255)
                             ggggCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTTAATATCATTGGGaaa                                                          >  1:3166051/1‑55 (MQ=255)
                              gggCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGa                                         <  1:1923872/71‑1 (MQ=255)
                                  aTCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAGTATCAt                                                                 >  1:475360/1‑43 (MQ=255)
                                   tCCGTTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCttt                                    <  1:1952505/71‑1 (MQ=255)
                                   tCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGAcc                                               <  1:1511434/60‑1 (MQ=255)
                                   tCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGAcc                                               <  1:935354/60‑1 (MQ=255)
                                   tCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCttt                                    <  1:626388/71‑1 (MQ=255)
                                    ccGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCttt                                    >  1:3642843/1‑70 (MQ=255)
                                     cGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCa                                              <  1:1708465/59‑1 (MQ=255)
                                     cGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCttt                                    >  1:1116878/1‑69 (MQ=255)
                                       gTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCttt                                    <  1:1175742/67‑1 (MQ=255)
                                       gTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAgc                                <  1:438456/71‑1 (MQ=255)
                                        ttGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGc                                                       <  1:1740455/47‑1 (MQ=255)
                                        ttGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGc                                                       <  1:1710309/47‑1 (MQ=255)
                                        ttGGCGATCAGGTGATAAAAACAGTGGGGAATATCATTGGGAAAAGc                                                       <  1:425853/47‑1 (MQ=255)
                                         tGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAgcgc                              >  1:2847596/1‑71 (MQ=255)
                                           gCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAgcgcg                             <  1:724448/70‑1 (MQ=255)
                                           gCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAgcgcg                             <  1:535544/70‑1 (MQ=255)
                                              aTCAGGTGATAAAAACAGTGGTGAGTATCATTGGGAAAAGCATACGGCCAGATGACCTTTTTGCGCGCGTc                         <  1:3521520/71‑1 (MQ=255)
                                              aTCAGGTGATAAAAACAGTGGTGAGTATCATTGGGAAAAGCATACGGCCAGATGACCTTTTTGCGCGCGTc                         <  1:2483407/71‑1 (MQ=255)
                                                 aGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTcggc                      >  1:494450/1‑71 (MQ=255)
                                                      aTAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTc                         >  1:1872817/1‑63 (MQ=255)
                                                         aaaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGa                 >  1:2314648/1‑68 (MQ=255)
                                                         aaaaCAGTGGTGAATATCATTGGGAAAAGCATACGACAAGATGACCTTTTAGCGCGCGTCGGCGGCGaa                >  1:65765/1‑69 (MQ=255)
                                                             cAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGttt          <  1:939175/71‑1 (MQ=255)
                                                                 tgtgAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGtt           <  1:283045/65‑1 (MQ=255)
                                                                  gtgAATATCATTGGGAAAAGCATACGACCAGATGACCttt                                    <  1:2075991/40‑1 (MQ=255)
                                                                  gtgAATATCATTGGGAAAAGCATACGACCAGATGACCttt                                    <  1:3026908/40‑1 (MQ=255)
                                                                  gtgAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAg                                 >  1:2422627/1‑43 (MQ=255)
                                                                   tgAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTc    <  1:657144/71‑1 (MQ=255)
                                                                   tgAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTc    <  1:1251442/71‑1 (MQ=255)
                                                                   tgAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTc    <  1:2670284/71‑1 (MQ=255)
                                                                   tgAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTc    <  1:2215432/71‑1 (MQ=255)
                                                                    gAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCt   >  1:3615764/1‑71 (MQ=255)
                                                                     aaTATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTcggc                      <  1:2278759/51‑1 (MQ=255)
                                                                     aaTATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGtt           <  1:1985585/62‑1 (MQ=255)
                                                                      atatCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCtt  <  1:366483/70‑1 (MQ=255)
                                                                      |                                                                     
GACCACTTCAAAAAAGTCAACGACACCTGGGGGCATCCGGTTGGCGATCAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTT  >  W3110S.gb/1094255‑1094394

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: