Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,007,501 A→C 35.0% N174K (AAT→AAG ycbX ← predicted 2Fe‑2S cluster‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,007,5010AC35.0% 11.0 / 48.2 40N174K (AAT→AAGycbXpredicted 2Fe‑2S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (14/12);  new base C (6/8);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGTTGTTGGAGATCACGTAACGAGGCT  >  W3110S.gb/1007432‑1007567
                                                                     |                                                                  
caccacATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGAt                                                                   <  1:1527603/71‑1 (MQ=255)
   cacaTCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGCTTgg                                                                <  1:525987/71‑1 (MQ=255)
       tCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGc                                                             <  1:4316079/70‑1 (MQ=255)
        cACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGCTTGGGGCg                                                            >  1:337911/1‑70 (MQ=255)
         aCCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGCTTGGGGCgg                                                           >  1:325283/1‑70 (MQ=255)
          ccAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCTCCGGAAACCACCAGCTTGGGGCGGaa                                                         <  1:4658556/71‑1 (MQ=255)
          ccAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGCTTGGGGCGGaa                                                         >  1:1356521/1‑71 (MQ=255)
          ccAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGCTTGGGGCGGaa                                                         <  1:4276563/71‑1 (MQ=255)
          ccAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGaa                                                         >  1:767905/1‑71 (MQ=255)
            aaTGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAAtt                                                       <  1:3102280/71‑1 (MQ=255)
             aTGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCTGAAACCACCAGATTGGGGCGGAATTg                                                      <  1:1452175/71‑1 (MQ=255)
             aTGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGAt                                                                   >  1:2198792/1‑58 (MQ=255)
             aTGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTg                                                      <  1:2765304/71‑1 (MQ=255)
             aTGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTg                                                      <  1:1154526/71‑1 (MQ=255)
                cGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGaa                                                         >  1:1311806/1‑65 (MQ=255)
                cGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGaa                                                         >  1:715272/1‑65 (MQ=255)
                cGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTc                                                   >  1:1494649/1‑71 (MQ=255)
                    tCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCAtt                                               <  1:3004337/71‑1 (MQ=255)
                    acACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCAtt                                               >  1:4917561/2‑71 (MQ=255)
                      aCTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAGACCACCAGCTTGGGGCGGAATTGCTCCAtttt                                             <  1:1811373/71‑1 (MQ=255)
                      aCTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCAtttt                                             <  1:1714794/71‑1 (MQ=255)
                          tCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTaca                                         >  1:2580060/1‑71 (MQ=255)
                            cAGCGATCTTCATCCCAGGCTGACGCGCCGGAAACCACCAGCTTGGGGCGGAATTGCTCCATTTTTaca                                         >  1:4674535/1‑69 (MQ=255)
                              gCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTACACTgg                                     >  1:4791218/1‑71 (MQ=255)
                               cGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTa                                           >  1:4842491/1‑64 (MQ=255)
                                   cttcttCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCAt                                                <  1:4001174/55‑1 (MQ=255)
                                    ttcttcCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTa                                           <  1:3357386/59‑1 (MQ=255)
                                    ttcttcCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTACACTg                                      <  1:2210236/64‑1 (MQ=255)
                                          ccAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACgt                          <  1:1069450/70‑1 (MQ=255)
                                           cAGGCTGACGCGCCGGAAACCACCAGCTTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACgt                          >  1:4425985/1‑69 (MQ=255)
                                              gCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTaca                                         >  1:1706811/1‑51 (MQ=255)
                                               cTGACGCGCCGGAAACCACCAGCTTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGTTGTTgg                    <  1:1512454/71‑1 (MQ=255)
                                               cTGACGCGCCGGAAACCACCAGCTTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGTTGTTgg                    <  1:1550265/71‑1 (MQ=255)
                                                tGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGttgttg                     >  1:440797/1‑69 (MQ=255)
                                                tGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGttgttg                     >  1:4173608/1‑69 (MQ=255)
                                                tGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTAAACTGGCCGGACAACGttgttg                     >  1:1809322/1‑69 (MQ=255)
                                                       ccGGAAACCACCAGCTTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACCACGTTGTTGGAGATCAc             <  1:5066995/70‑1 (MQ=255)
                                                            aaCCACCAGCTTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGTTGTTGGAGATCACGTAAc        >  1:2719194/1‑70 (MQ=255)
                                                               caccaGCTTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGTTGTTGgag                  <  1:4270289/57‑1 (MQ=255)
                                                                 ccaGATTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGTTGTTGGAGATCACGTAACGAGGCt  >  1:5010264/1‑71 (MQ=255)
                                                                     |                                                                  
CACCACATCACCAATGCGAATCACTTTCCAGCGATCTTCTTCCCAGGCTGACGCGCCGGAAACCACCAGATTGGGGCGGAATTGCTCCATTTTTACACTGGCCGGACAACGTTGTTGGAGATCACGTAACGAGGCT  >  W3110S.gb/1007432‑1007567

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: