Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,232,470 A→G 16.0% L111L (TTG→CTG)  ygjM ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,232,4700AG16.0% 118.7 / 16.3 50L111L (TTG→CTG) ygjMpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/13);  new base G (4/4);  total (33/17)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATA  >  W3110S.gb/3232400‑3232540
                                                                      |                                                                      
gCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:3628299/71‑1 (MQ=255)
gCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:200728/71‑1 (MQ=255)
gCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:2074074/71‑1 (MQ=255)
   ggAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:3163305/68‑1 (MQ=255)
   ggAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:289903/68‑1 (MQ=255)
   ggAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:270749/68‑1 (MQ=255)
   ggAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:2215239/68‑1 (MQ=255)
     agaAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgc                                                                   >  1:2621398/1‑71 (MQ=255)
       aaaTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAGAACGCGt                                                                 <  1:2797410/71‑1 (MQ=255)
       aaaTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGt                                                                 <  1:2379356/71‑1 (MQ=255)
         aTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAGAACGCGTGa                                                               >  1:3413388/1‑71 (MQ=255)
         aTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAGAACGCGTGa                                                               >  1:2832249/1‑71 (MQ=255)
         aTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAGAACGCGTGa                                                               >  1:2729260/1‑71 (MQ=255)
         aTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGa                                                               >  1:988304/1‑71 (MQ=255)
         aTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGa                                                               >  1:721868/1‑71 (MQ=255)
          tGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGac                                                              <  1:2369153/71‑1 (MQ=255)
            ccGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGa                                                               >  1:845499/1‑68 (MQ=255)
            ccGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGa                                                               >  1:1355883/1‑68 (MQ=255)
              gAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCa                                                          <  1:259109/71‑1 (MQ=255)
                   gcgTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGa                                                      <  1:2585470/70‑1 (MQ=255)
                        gCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACt                                                >  1:1264670/1‑71 (MQ=255)
                             ttttttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:2449962/42‑1 (MQ=255)
                             ttttttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                                                        <  1:1266858/42‑1 (MQ=255)
                               ttttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAGAACGCGTGACACCATAGATTTACTGCCAAtt                                         >  1:2348728/1‑71 (MQ=255)
                                       tgtTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCgg                                    >  1:573687/1‑68 (MQ=255)
                                         ttCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCgg                                    >  1:2385823/1‑66 (MQ=255)
                                         ttCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCa                                  <  1:3361629/68‑1 (MQ=255)
                                           ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGa                                                      >  1:2965301/1‑46 (MQ=255)
                                           ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGa                                                      >  1:1039700/1‑46 (MQ=255)
                                           ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGa                                                      >  1:1218078/1‑46 (MQ=255)
                                           ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGa                                                      >  1:555117/1‑46 (MQ=255)
                                           ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCgg                                    >  1:2269366/1‑64 (MQ=255)
                                           ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCgg                                    >  1:1398808/1‑64 (MQ=255)
                                           ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCgg                                    >  1:1578252/1‑64 (MQ=255)
                                           ccAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCg                             >  1:2388467/1‑71 (MQ=255)
                                             aGCGTTAATTTCCTCTTCCCGCTCAGAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGa                           <  1:2287024/71‑1 (MQ=255)
                                              gCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCa                                            <  1:3768293/53‑1 (MQ=255)
                                              gCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGaa                          <  1:297347/71‑1 (MQ=255)
                                                 ttAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc                                   >  1:183346/1‑59 (MQ=255)
                                                 ttAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc                                   >  1:809088/1‑59 (MQ=255)
                                                 ttAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc                                   >  1:3339800/1‑59 (MQ=255)
                                                  tAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAggg                      >  1:778534/1‑71 (MQ=255)
                                                   aaTTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGt                     >  1:3960535/1‑71 (MQ=255)
                                                    aTTTCCTCTTCCCGCTCAGAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGt                     <  1:3613891/70‑1 (MQ=255)
                                                     tttCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTa                   <  1:764164/71‑1 (MQ=255)
                                                       tCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTaaa                 >  1:2595217/1‑71 (MQ=255)
                                                            ttCCCGCTCAGAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCa              <  1:2557027/69‑1 (MQ=255)
                                                            ttCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCa              <  1:3756712/69‑1 (MQ=255)
                                                              cccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGAt                               >  1:2927751/1‑50 (MQ=255)
                                                              cccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCat             <  1:372140/68‑1 (MQ=255)
                                                              cccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCat             <  1:2223104/68‑1 (MQ=255)
                                                              cccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATAtt          >  1:649512/1‑71 (MQ=255)
                                                              cccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATAtt          >  1:2110289/1‑71 (MQ=255)
                                                              cccGCTCAAAACGCGTGACACCATAGATTTACTGCCAACTTCCGGCAGATCGGAAAGGGTTAAACCat             <  1:1964875/68‑1 (MQ=255)
                                                                  cTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATc      <  1:2753943/71‑1 (MQ=255)
                                                                    cAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTg         >  1:1219002/1‑66 (MQ=255)
                                                                    cAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGa        >  1:3496540/1‑67 (MQ=255)
                                                                    cAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGa        >  1:3491800/1‑67 (MQ=255)
                                                                    cAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGa        >  1:1236356/1‑67 (MQ=255)
                                                                      aaaCGCGTGACACCATAGATTTACTTCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATa  <  1:2649259/71‑1 (MQ=255)
                                                                      |                                                                      
GCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATA  >  W3110S.gb/3232400‑3232540

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: