Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 701,210 A→G 18.2% A269A (GCT→GCC nagC ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb701,2100AG18.2% 91.3 / 13.6 44A269A (GCT→GCCnagCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (16/20);  new base G (5/3);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGC  >  W3110S.gb/701142‑701280
                                                                    |                                                                      
aGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGc                                                                      <  1:3071297/71‑1 (MQ=255)
   ggCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGa                                                                    >  1:137996/1‑70 (MQ=255)
   ggCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGa                                                                    >  1:206907/1‑70 (MQ=255)
    gCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATa                                                                  <  1:3286573/71‑1 (MQ=255)
    gCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATa                                                                  <  1:54008/71‑1 (MQ=255)
                  tAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCATTGGCGGCGATa                                                                  >  1:2364428/1‑57 (MQ=255)
                  tAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCATTGGCGGCGATAg                                                                 >  1:3327231/1‑58 (MQ=255)
                  tAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCATTGGCGGCGATAg                                                                 >  1:3322712/1‑58 (MQ=255)
                  tAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCATTGGCGGCGATAg                                                                 >  1:2574518/1‑58 (MQ=255)
                  tAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCATTGGCGGCGATAg                                                                 >  1:2448050/1‑58 (MQ=255)
                      ccTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCaaa                                                 <  1:967092/70‑1 (MQ=255)
                       cTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGt                                                                <  1:1310696/54‑1 (MQ=255)
                       cTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCaa                                                  <  1:1748948/68‑1 (MQ=255)
                        tGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAg                                                      <  1:547753/63‑1 (MQ=255)
                        tGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCaa                                                  <  1:3087723/67‑1 (MQ=255)
                          cTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGAt                                                                   <  1:2349899/48‑1 (MQ=255)
                              aaCAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGc                                         <  1:427992/70‑1 (MQ=255)
                              aaCAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGc                                         <  1:3304949/70‑1 (MQ=255)
                              aaCAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGc                                         <  1:1016942/70‑1 (MQ=255)
                              aaCAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGc                                         <  1:304476/70‑1 (MQ=255)
                              aaCAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGCTGCCGc                                         <  1:492044/70‑1 (MQ=255)
                                 aGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAg                                                >  1:2099501/1‑60 (MQ=255)
                                 aGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTg                                     >  1:3189489/1‑71 (MQ=255)
                                    ttCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTcc                                                            >  1:30193/1‑45 (MQ=255)
                                    ttCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGAAAGTTTCCAGGCAGCCAAAGTTGCCGCAGTgg                                    >  1:3675562/1‑69 (MQ=255)
                                     tCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAg                                 >  1:3798019/1‑71 (MQ=255)
                                          acCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCa                                                           >  1:2043156/1‑40 (MQ=255)
                                          acCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCa                                                           >  1:709455/1‑40 (MQ=255)
                                             cGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCa                                                       >  1:1438726/1‑41 (MQ=255)
                                              gttgttCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCAcc                         >  1:2739123/1‑70 (MQ=255)
                                                tgttCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCa                       <  1:3720863/70‑1 (MQ=255)
                                                  ttCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGc                     >  1:540461/1‑70 (MQ=255)
                                                   tCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCgg                   >  1:1668865/1‑71 (MQ=255)
                                                            gCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCgg                   <  1:3810597/62‑1 (MQ=255)
                                                             cattGGCGGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGa         <  1:287546/69‑1 (MQ=255)
                                                             cattGGCGGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGa         <  1:2596305/69‑1 (MQ=255)
                                                             cattGGCGGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGa         <  1:2334773/69‑1 (MQ=255)
                                                              gTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGaa        <  1:2350999/71‑1 (MQ=255)
                                                              gTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGaa        <  1:1301376/71‑1 (MQ=255)
                                                               ttGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGt                  <  1:626275/60‑1 (MQ=255)
                                                                 ggCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAatat     >  1:137110/1‑71 (MQ=255)
                                                                  gCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAatat     <  1:1394981/70‑1 (MQ=255)
                                                                  gCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATg    >  1:321807/1‑71 (MQ=255)
                                                                  gCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATg    >  1:1136035/1‑71 (MQ=255)
                                                                    aGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGc  >  1:154780/1‑71 (MQ=255)
                                                                    |                                                                      
AGCGGCACGCGGCTCTGGTAGCCCTGCTTTAACAGATTCAACACCCGTTGTTCAATGGCAGCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGC  >  W3110S.gb/701142‑701280

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: