Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,532,412 T→C 38.6% M236T (ATG→ACG)  cpxA → sensory histidine kinase in two‑component regulatory system with CpxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,532,4120TC38.6% 28.3 / 48.9 44M236T (ATG→ACG) cpxAsensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base T (17/10);  new base C (7/10);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CGGGGCCACAGGAATTCCTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGCAACT  >  W3110S.gb/3532343‑3532481
                                                                     |                                                                     
cGGGGCCACAGGAATTCCTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCatgatg                                                                      <  1:5000403/71‑1 (MQ=255)
       acaGGAATTCCTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACctct                                                                >  1:3887654/1‑70 (MQ=255)
         aGGAATTCCTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACctct                                                                >  1:175452/1‑68 (MQ=255)
          ggAATTCCTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTcag                                                             >  1:1919880/1‑70 (MQ=255)
              ttCCTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCAtgatga                                                                     <  1:797441/58‑1 (MQ=255)
                ccTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTc                                                      >  1:829221/1‑71 (MQ=255)
                 cTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCt                                                     >  1:233529/1‑71 (MQ=255)
                    gCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCt                                                  >  1:4084551/1‑71 (MQ=255)
                    gCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCt                                                  >  1:2593329/1‑71 (MQ=255)
                       gCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTc                                               <  1:1525810/71‑1 (MQ=255)
                        cAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCt                                              <  1:4563514/71‑1 (MQ=255)
                        cAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCt                                              <  1:3054382/71‑1 (MQ=255)
                           gTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGat                                           <  1:4000070/71‑1 (MQ=255)
                            tGCCAGTTTTAACCAGATGGTAACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGata                                          <  1:2068572/71‑1 (MQ=255)
                             gCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTcag                                                             >  1:2890948/1‑51 (MQ=255)
                              ccAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGat                                           <  1:1084530/68‑1 (MQ=255)
                                aGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGatat                                         >  1:4748861/1‑68 (MQ=255)
                                 gTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGATAtctct                                     <  1:2236716/71‑1 (MQ=255)
                                 gTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGATAtctct                                     <  1:2484688/71‑1 (MQ=255)
                                  ttttAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTca                                                              <  1:1356762/45‑1 (MQ=255)
                                  ttttAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTca                                                              <  1:4912334/45‑1 (MQ=255)
                                  ttttAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTca                                                              <  1:2077234/45‑1 (MQ=255)
                                   tttAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCa                                   >  1:647695/1‑71 (MQ=255)
                                   tttAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCa                                   >  1:3553624/1‑71 (MQ=255)
                                    ttAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGa                                            >  1:1739026/1‑61 (MQ=255)
                                    ttAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCAc                                  >  1:245513/1‑71 (MQ=255)
                                    ttAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCAc                                  >  1:2877921/1‑71 (MQ=255)
                                    ttAACCAGATGGTCACCGCGCTGGAGCGCATGACGACCTCTAAGCAGCGTCTGCTTTCTGATATCTCTCAc                                  >  1:2616037/1‑71 (MQ=255)
                                     tAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCAc                                  >  1:3024800/1‑70 (MQ=255)
                                      aaCCAGATGGTCACCGCGCTGGAGCGCATGATGAc                                                                    >  1:3244905/1‑35 (MQ=255)
                                                tCACCGCGCTGTAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCAcc                      >  1:4208003/1‑71 (MQ=255)
                                                tCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCAcc                      >  1:4473343/1‑71 (MQ=255)
                                                tCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCAcc                      >  1:4126506/1‑71 (MQ=255)
                                                 cACCGCGCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCAccc                     >  1:1605237/1‑71 (MQ=255)
                                                  aCCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCAccc                     <  1:2272869/70‑1 (MQ=255)
                                                    cgcgCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGc                  >  1:4858327/1‑71 (MQ=255)
                                                    cgcgCTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGc                  >  1:860215/1‑71 (MQ=255)
                                                     gcgcTGGAGCGCATGACGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCt                 <  1:1897007/71‑1 (MQ=255)
                                                     gcgcTCGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCt                 <  1:1778872/71‑1 (MQ=255)
                                                           gAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGGTGCGCACCCCGCTGAcgc            <  1:2708634/70‑1 (MQ=255)
                                                           gAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGAcgcg           <  1:1149263/71‑1 (MQ=255)
                                                                cATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCa                                   >  1:3137048/1‑42 (MQ=255)
                                                                cATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGc      <  1:3949084/71‑1 (MQ=255)
                                                                cATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGc      <  1:3795452/71‑1 (MQ=255)
                                                                cATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGc      <  1:1706442/71‑1 (MQ=255)
                                                                cATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGc      <  1:1214548/71‑1 (MQ=255)
                                                                  tgacgaCCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACg                                 >  1:543085/1‑42 (MQ=255)
                                                                  tgacgaCCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGAc              <  1:3359150/61‑1 (MQ=255)
                                                                  tgacgaCCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGCaa    >  1:120981/1‑71 (MQ=255)
                                                                   gatgaCCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGCAAc   >  1:2047280/1‑71 (MQ=255)
                                                                    atgaCCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGa               >  1:4090670/1‑58 (MQ=255)
                                                                    atgaCCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGa               >  1:2683200/1‑58 (MQ=255)
                                                                    atgaCCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGCAACt  >  1:1463370/1‑71 (MQ=255)
                                                                     |                                                                     
CGGGGCCACAGGAATTCCTTGCCGCAGGTGCCAGTTTTAACCAGATGGTCACCGCGCTGGAGCGCATGATGACCTCTCAGCAGCGTCTGCTTTCTGATATCTCTCACGAGCTGCGCACCCCGCTGACGCGTCTGCAACT  >  W3110S.gb/3532343‑3532481

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: