Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 701,092 T→C 23.2% S309G (AGT→GGT)  nagC ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb701,0920TC23.2% 106.7 / 31.1 56S309G (AGT→GGT) nagCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (18/25);  new base C (8/5);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGCTCTGG  >  W3110S.gb/701028‑701159
                                                                |                                                                   
aTAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGc                                                                <  1:781491/70‑1 (MQ=255)
     aaTGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTg                                                           <  1:1156591/70‑1 (MQ=255)
      aTGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGtt                                                         <  1:1704189/71‑1 (MQ=255)
        ggCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACCATCGCCTTTGTTcg                                                       >  1:3944361/1‑71 (MQ=255)
        ggCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACCATCGCCTTTGTTcg                                                       >  1:2514798/1‑71 (MQ=255)
         gCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTcgc                                                      >  1:2807877/1‑71 (MQ=255)
         gCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACCATCGCCTTTGTTcgc                                                      >  1:2423768/1‑71 (MQ=255)
               gTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCACGGCttt                                                >  1:4820426/1‑71 (MQ=255)
                ttttACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGc                                                                <  1:3624543/54‑1 (MQ=255)
                ttttACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGc                                                                <  1:4168942/54‑1 (MQ=255)
                 tttACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTg                                               >  1:820802/1‑70 (MQ=255)
                  ttACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTAt                                                                   >  1:4714532/1‑49 (MQ=255)
                   tACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACCATCGCCTTTGTTCGCGGCTTTg                                               >  1:950604/1‑68 (MQ=255)
                          aTGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtt                                      <  1:442230/70‑1 (MQ=255)
                          aTGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtt                                      <  1:1027698/70‑1 (MQ=255)
                           tGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTg                                               <  1:4236091/60‑1 (MQ=255)
                           tGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt                                    <  1:3790276/71‑1 (MQ=255)
                           tGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt                                    <  1:1984671/71‑1 (MQ=255)
                           tGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtt                                      <  1:3968552/69‑1 (MQ=255)
                            gacgacCGACATACTCAATTACTTCCGACGCCAGACTAt                                                                   >  1:4241902/1‑39 (MQ=255)
                             acgacCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGa                                  <  1:313560/71‑1 (MQ=255)
                                  cgacATACTCAATTACTTCCGACGCCAGACCATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTg                             <  1:930839/71‑1 (MQ=255)
                                  cgacATACTCAATTACTTCCGACCCCAGACCATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTg                             <  1:3344146/71‑1 (MQ=255)
                                   gacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGttt                                     >  1:3883404/1‑62 (MQ=255)
                                   gacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGatt                                     >  1:4298693/1‑62 (MQ=255)
                                      aTACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAgt                         >  1:1353268/1‑71 (MQ=255)
                                       tACTCAATTACTTCCGACGCCAGACTATCGCCTTTg                                                           <  1:3004720/36‑1 (MQ=255)
                                       tACTCAATTACTTCCGACGCCAGACTATCGCCTTTg                                                           <  1:3408494/36‑1 (MQ=255)
                                       tACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTg                                   <  1:939602/60‑1 (MQ=255)
                                       tACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTg                                   <  1:3030041/60‑1 (MQ=255)
                                       tACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTg                                   <  1:3035881/60‑1 (MQ=255)
                                        aCTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTg                                   <  1:554740/59‑1 (MQ=255)
                                            aaTTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAgtcgtc                     >  1:3926854/1‑69 (MQ=255)
                                             aTTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATa                                         >  1:1471242/1‑48 (MQ=255)
                                             aTTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATa                                         >  1:2186685/1‑48 (MQ=255)
                                              ttACTTCCGACGCCAGACCATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCa                           <  1:3626497/61‑1 (MQ=255)
                                                aCTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAg                                            >  1:3461657/1‑42 (MQ=255)
                                                 cTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCgg               <  1:5028822/70‑1 (MQ=255)
                                                 cTTCCGACGCCAGACCATCGCCTTTGTTCGCGGCtt                                                 >  1:322041/1‑36 (MQ=255)
                                                    ccGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATgg                               >  1:2425149/1‑51 (MQ=255)
                                                    ccGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATgg                               >  1:2033861/1‑51 (MQ=255)
                                                    ccGACGCCAGACCATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGAt                                 >  1:686151/1‑49 (MQ=255)
                                                    ccGACGCCAGACCATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCgcc              >  1:4173291/1‑66 (MQ=255)
                                                         gCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGct      >  1:2348517/1‑71 (MQ=255)
                                                          ccAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt                                    <  1:2094401/40‑1 (MQ=255)
                                                          ccAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt                                    <  1:1812914/40‑1 (MQ=255)
                                                          ccAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt                                    <  1:2801218/40‑1 (MQ=255)
                                                          ccAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt                                    <  1:2100298/40‑1 (MQ=255)
                                                          ccAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGctc     >  1:4736143/1‑71 (MQ=255)
                                                           cAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGctct    >  1:1426016/1‑71 (MQ=255)
                                                           cAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGctc     <  1:4352498/70‑1 (MQ=255)
                                                           cAGACCATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGctct    <  1:4683607/71‑1 (MQ=255)
                                                           cAGACCATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGctct    >  1:4961010/1‑71 (MQ=255)
                                                            aGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGCTCTg   <  1:179712/71‑1 (MQ=255)
                                                            aGACCATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGa                                  <  1:4254233/40‑1 (MQ=255)
                                                               cTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGCTCTgg  >  1:2448733/1‑69 (MQ=255)
                                                                tATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAg                  <  1:500326/52‑1 (MQ=255)
                                                                tATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAg                  <  1:814044/52‑1 (MQ=255)
                                                                tATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAg                  <  1:2060412/52‑1 (MQ=255)
                                                                |                                                                   
ATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTCGTCCAGCGGCACGCGGCTCTGG  >  W3110S.gb/701028‑701159

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: