Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,655,568 C→T 12.9% G58S (GGT→AGT)  fdx ← [2Fe‑2S] ferredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,655,5680CT12.9% 170.2 / 14.7 62G58S (GGT→AGT) fdx[2Fe‑2S] ferredoxin
Reads supporting (aligned to +/- strand):  ref base C (33/21);  new base T (3/5);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

AGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGA  >  W3110S.gb/2655504‑2655633
                                                                |                                                                 
aGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacg                                                             <  1:699014/71‑1 (MQ=255)
aGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacg                                                             <  1:190360/71‑1 (MQ=255)
    cccAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacgaac                                                          >  1:528606/1‑70 (MQ=255)
    cccAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacgaac                                                          >  1:650107/1‑70 (MQ=255)
     ccAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacgaacga                                                        <  1:3698752/71‑1 (MQ=255)
        ggCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATg                                                      >  1:2236295/1‑70 (MQ=255)
        ggCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCATCACGAACGATg                                                      >  1:1703125/1‑70 (MQ=255)
             tGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACTTTCACGAACGATGCAGt                                                  <  1:258515/69‑1 (MQ=255)
             tGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACTTTCACGAACGATGCAGt                                                  <  1:768248/69‑1 (MQ=255)
             tGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACTTTCACGAACGATGCAGt                                                  <  1:1352100/69‑1 (MQ=255)
               tCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGc                                               <  1:1784140/70‑1 (MQ=255)
                      tGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacgaacg                                                         <  1:2660210/53‑1 (MQ=255)
                      tGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacgaacg                                                         <  1:2546190/53‑1 (MQ=255)
                      tGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATg                                                      <  1:889494/56‑1 (MQ=255)
                      tGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATg                                                      <  1:429781/56‑1 (MQ=255)
                      tGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAggt                                           <  1:2612273/67‑1 (MQ=255)
                        tcgtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAggtgg                                         >  1:3134382/1‑67 (MQ=255)
                        tcgtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAggtgg                                         >  1:2307952/1‑67 (MQ=255)
                        tcgtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGgtggtg                                       >  1:2026937/1‑69 (MQ=255)
                        tcgtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGc                                      >  1:3923236/1‑70 (MQ=255)
                        tcgtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCa                                     <  1:300637/71‑1 (MQ=255)
                        tcgtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCa                                     <  1:1416477/71‑1 (MQ=255)
                        tcgtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCa                                     <  1:1947872/71‑1 (MQ=255)
                        tcgtcTACCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAggtgg                                         >  1:214671/1‑67 (MQ=255)
                          gtcTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacgaac                                                          <  1:2807685/48‑1 (MQ=255)
                             ttCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACTTTCacgaacga                                                        >  1:2192265/1‑47 (MQ=255)
                             ttCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACTTTCacgaacga                                                        >  1:2925935/1‑47 (MQ=255)
                             ttCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACTTTCacgaacga                                                        >  1:1361293/1‑47 (MQ=255)
                             ttCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACTTTCACGAACGATGCAGTGGCa                                              <  1:933078/57‑1 (MQ=255)
                             ttCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACTTTCACGAACGATGCAGTGGCa                                              <  1:2585239/57‑1 (MQ=255)
                               ccTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCa                                                               >  1:3531580/1‑38 (MQ=255)
                               ccTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCa                                                               >  1:3583203/1‑38 (MQ=255)
                               ccTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCa                                                               >  1:618770/1‑38 (MQ=255)
                               ccTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCa                                                               >  1:2382305/1‑38 (MQ=255)
                               ccTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGca                                 <  1:701306/68‑1 (MQ=255)
                               ccTGCTCTGAGCTTTACGGCAGTGAGTCAAAACCTTCa                                                               >  1:3177334/1‑38 (MQ=255)
                                 tGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGGACGATGCAGTGGCAGGTGGTGCAAGCACAGGa                            >  1:2201937/1‑71 (MQ=255)
                                 tGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGa                            >  1:1029180/1‑71 (MQ=255)
                                  gatctGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGAt                                                       >  1:1059308/3‑43 (MQ=255)
                                  gCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGAt                                                       >  1:2768302/1‑43 (MQ=255)
                                      tGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAATCACAGGAttt                         <  1:564962/69‑1 (MQ=255)
                                      tGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAttttt                       >  1:885245/1‑71 (MQ=255)
                                      tGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAttttt                       >  1:2518969/1‑71 (MQ=255)
                                        aGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGt                                                  >  1:2027664/1‑42 (MQ=255)
                                                      gAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCaa        >  1:2279869/1‑70 (MQ=255)
                                                       aGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTcaca                   >  1:808806/1‑58 (MQ=255)
                                                       aGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCa         <  1:3487647/68‑1 (MQ=255)
                                                         tCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGc                <  1:435159/59‑1 (MQ=255)
                                                         tCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGc                <  1:26241/59‑1 (MQ=255)
                                                         tCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGc                <  1:2804244/59‑1 (MQ=255)
                                                         tCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAAtctc    >  1:213017/1‑71 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAt                           >  1:482944/1‑46 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAt                           >  1:386980/1‑46 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAt                           >  1:3518578/1‑46 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAt                           >  1:1506209/1‑46 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCgt              >  1:2426381/1‑59 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCgt              >  1:949206/1‑59 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCgt              >  1:991543/1‑59 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCg               >  1:2008457/1‑58 (MQ=255)
                                                           aaaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGa  >  1:3819673/1‑71 (MQ=255)
                                                            aaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGcac                                <  1:3963476/40‑1 (MQ=255)
                                                            aaaCCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGcac                                <  1:1299552/40‑1 (MQ=255)
                                                                cTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAttttt                       >  1:2179205/1‑45 (MQ=255)
                                                                |                                                                 
AGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGTGTTCAATCTCGA  >  W3110S.gb/2655504‑2655633

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: