Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,081,459 A→G 25.4% F36L (TTT→CTT)  yeeX ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,081,4590AG25.4% 115.3 / 43.8 67F36L (TTT→CTT) yeeXconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/30);  new base G (11/6);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 9.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACGGATCTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAATGTGAGAGTCGGA  >  W3110S.gb/2081391‑2081528
                                                                    |                                                                     
aCGGATCTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAAt                                                                      <  1:2565719/70‑1 (MQ=255)
   gATCTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCa                                                                  <  1:557383/71‑1 (MQ=255)
     tCTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTcc                                                                   <  1:324786/68‑1 (MQ=255)
      cTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTa                                                               >  1:1317883/1‑71 (MQ=255)
         tttCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAg                                                                 <  1:1067431/66‑1 (MQ=255)
         tttCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAg                                                                 <  1:1339754/66‑1 (MQ=255)
         tttCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAg                                                                 <  1:3116191/66‑1 (MQ=255)
         tttCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACg                                                             <  1:3425698/70‑1 (MQ=255)
         tttCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACg                                                             <  1:2944437/70‑1 (MQ=255)
               aaCGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTgg                                                       <  1:561689/70‑1 (MQ=255)
               aaCGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTgg                                                       <  1:789766/70‑1 (MQ=255)
                  gTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATg                                                   <  1:3854304/71‑1 (MQ=255)
                   tCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTgg                                                       <  1:3384534/66‑1 (MQ=255)
                   tCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTgg                                                       <  1:2035175/66‑1 (MQ=255)
                   tCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTgg                                                       <  1:3092732/66‑1 (MQ=255)
                   tCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTgg                                                       <  1:139713/66‑1 (MQ=255)
                     cTGAATTTCACGTTGCAGTTTGTTCTTACGACTGAACAGACGAACAAGTTCCAGTACGTCCTGGaa                                                     >  1:2913176/1‑66 (MQ=255)
                     cTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGaa                                                     >  1:2607045/1‑66 (MQ=255)
                     cTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGaa                                                     >  1:2873234/1‑66 (MQ=255)
                     cTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAg                                                >  1:343482/1‑71 (MQ=255)
                       gAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGc                                              >  1:1924645/1‑71 (MQ=255)
                       gAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGc                                              >  1:3468426/1‑71 (MQ=255)
                       gAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGc                                              >  1:2637278/1‑71 (MQ=255)
                       gAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGc                                              >  1:1551449/1‑71 (MQ=255)
                        aaTTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGc                                              <  1:2620452/70‑1 (MQ=255)
                         aTTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGaa                                                 <  1:1996615/66‑1 (MQ=255)
                          tttCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATACCAGTACGTCCTGGAATGAAGGCTTg                                           >  1:1313302/1‑71 (MQ=255)
                           ttCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAgg                                               >  1:421947/1‑66 (MQ=255)
                            tCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTgg                                          >  1:2905752/1‑70 (MQ=255)
                             cACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTgg                                          >  1:958428/1‑69 (MQ=255)
                             cACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGt                                         <  1:2673395/70‑1 (MQ=255)
                             cACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGt                                         <  1:797846/70‑1 (MQ=255)
                             cACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGt                                         <  1:1600622/70‑1 (MQ=255)
                              aCGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTa                                        <  1:1291019/70‑1 (MQ=255)
                               cGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGt                                      <  1:1167809/71‑1 (MQ=255)
                               cGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGt                                      <  1:2564978/71‑1 (MQ=255)
                               cGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGt                                      <  1:1692229/71‑1 (MQ=255)
                               cGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGt                                      <  1:81412/71‑1 (MQ=255)
                               cGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGt                                      <  1:913385/71‑1 (MQ=255)
                                      gTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtt                               <  1:1412459/71‑1 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             >  1:1072032/1‑70 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             >  1:709424/1‑70 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             >  1:179184/1‑70 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             >  1:1453916/1‑70 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAATTCCAGTa                                                               >  1:1146383/1‑36 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             >  1:3732136/1‑70 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             >  1:989196/1‑70 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             >  1:3129350/1‑70 (MQ=255)
                                         tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             <  1:2797259/70‑1 (MQ=255)
                                          gTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                           >  1:1546811/1‑71 (MQ=255)
                                          gTTCTTACGACGGAACAGACGAACAAATTACAGTACGTCCTGGAATGAAGGCTTgg                                          <  1:1907242/56‑1 (MQ=255)
                                           ttCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtt                               <  1:3722542/66‑1 (MQ=255)
                                           ttCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTtat                          <  1:2682228/71‑1 (MQ=255)
                                           ttCTTACGACGAAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                           <  1:965175/70‑1 (MQ=255)
                                             cTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATAc                        >  1:1627781/1‑71 (MQ=255)
                                                   acgGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             <  1:1417679/60‑1 (MQ=255)
                                                      gAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             <  1:1848622/57‑1 (MQ=255)
                                                      gAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATaa               >  1:1851749/1‑71 (MQ=255)
                                                      gAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATaa               >  1:316456/1‑71 (MQ=255)
                                                      gAACAGACGAACAAATTCCAGTAAGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATaa               >  1:3938316/1‑71 (MQ=255)
                                                         cAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                             <  1:3577297/54‑1 (MQ=255)
                                                          aGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAAtgt            >  1:2753974/1‑70 (MQ=255)
                                                           gACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAAtgt            >  1:722243/1‑69 (MQ=255)
                                                              gAACAAATTGCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAATGTGagag       <  1:845507/71‑1 (MQ=255)
                                                                 cAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAATGTGAGAGTc     >  1:1319360/1‑70 (MQ=255)
                                                                 cAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                           >  1:1041201/1‑48 (MQ=255)
                                                                 cAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAAt              >  1:27136/1‑61 (MQ=255)
                                                                  aaaTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAAtgtg           <  1:774934/63‑1 (MQ=255)
                                                                   aaTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                           <  1:2422114/46‑1 (MQ=255)
                                                                   aaTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                           <  1:1412823/46‑1 (MQ=255)
                                                                   aaTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAATGTGAGAGTCGGa  >  1:3109978/1‑71 (MQ=255)
                                                                    |                                                                     
ACGGATCTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAATGTGAGAGTCGGA  >  W3110S.gb/2081391‑2081528

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: