Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,393,339 A→G 8.9% N30S (AAT→AGT)  yfbO → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,393,3390AG 8.9% 178.8 / 9.1 67N30S (AAT→AGT) yfbOhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base G (2/4);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

ATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATT  >  W3110S.gb/2393269‑2393406
                                                                      |                                                                   
aTTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                     >  1:3892944/1‑71 (MQ=255)
aTTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                     >  1:1993212/1‑71 (MQ=255)
    gTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                 <  1:1674894/71‑1 (MQ=255)
    gTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                 <  1:1996960/71‑1 (MQ=255)
    gTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                 <  1:1921101/71‑1 (MQ=255)
    gTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATc                                                                 <  1:1896383/71‑1 (MQ=255)
    gTTTGATGAATATATTATTTCAGGCAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATc                                                                 <  1:3239474/71‑1 (MQ=255)
     tttGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCa                                                                <  1:3574328/71‑1 (MQ=255)
       tgatgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCaa                                                               <  1:564185/70‑1 (MQ=255)
          tgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGt                                                           <  1:2812004/71‑1 (MQ=255)
          tgaATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGt                                                           <  1:2083429/71‑1 (MQ=255)
              tatatTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATg                                                       >  1:1091907/1‑71 (MQ=255)
              tatatTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATg                                                       >  1:2133635/1‑71 (MQ=255)
              tatatTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATg                                                       >  1:772702/1‑71 (MQ=255)
                tattatTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGtt                                                     <  1:3311490/71‑1 (MQ=255)
                tattatTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGtt                                                     <  1:855469/71‑1 (MQ=255)
                tattatTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGtt                                                     <  1:4390058/71‑1 (MQ=255)
                  ttattTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                     >  1:885566/1‑53 (MQ=255)
                       tCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTaataat                                               <  1:1870757/70‑1 (MQ=255)
                        cAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATcc                                             >  1:4156752/1‑71 (MQ=255)
                        cAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATcc                                             >  1:495476/1‑71 (MQ=255)
                         aGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCAATGGTGATGTTAATAATCCt                                            <  1:1317383/71‑1 (MQ=255)
                         aGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTta                                                    >  1:3193051/1‑63 (MQ=255)
                         aGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATcc                                             <  1:3677212/70‑1 (MQ=255)
                           ggAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                          <  1:4437440/71‑1 (MQ=255)
                           ggAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                          >  1:696889/1‑71 (MQ=255)
                              ataatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                            <  1:237409/66‑1 (MQ=255)
                               taatTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                     >  1:4003513/1‑40 (MQ=255)
                               taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaa                                                   <  1:2900195/58‑1 (MQ=255)
                               taatTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATAct                                      >  1:2664662/1‑71 (MQ=255)
                               taatTATGACACCTTTAGAACGCATTACACAAATTGTTAATATCAATGGTGATGTtaa                                                   <  1:2059528/58‑1 (MQ=255)
                                aatTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                     >  1:2284310/1‑39 (MQ=255)
                                aatTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                     >  1:2146356/1‑39 (MQ=255)
                                aatTATGACACCTTTAGAACGCATTACCCAACTTGTTaa                                                                     >  1:3255152/1‑39 (MQ=255)
                                 atTATGACACCTTTAGAACGCATTACCCAACTTGTTAGTATCAATGGTGATGTTAATAATCCTGATActcc                                    >  1:2434898/1‑71 (MQ=255)
                                 atTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGa                                                         >  1:577752/1‑50 (MQ=255)
                                 atTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATc                                              >  1:4398820/1‑61 (MQ=255)
                                 atTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATActcc                                    >  1:5476/1‑71 (MQ=255)
                                  ttATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATActcct                                   >  1:3096181/1‑71 (MQ=255)
                                    aTGACACCTTTAGAACGCATTACCCAACTTGTTAATATCaa                                                               >  1:2752767/1‑41 (MQ=255)
                                        cacCTTTAGAACGCATTACCCAACTTGTTAGTATCAATGGTGATGTTAATAATCCTGATACTCCTCGTcc                              <  1:2515740/70‑1 (MQ=255)
                                          ccTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTc                               <  1:1888510/67‑1 (MQ=255)
                                          ccTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCAtt                           >  1:3993436/1‑71 (MQ=255)
                                           atttAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTg                          >  1:893350/2‑71 (MQ=255)
                                            tttAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGt                         >  1:1295632/1‑71 (MQ=255)
                                             ttAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTaataatgc                                             >  1:2969177/1‑48 (MQ=255)
                                             ttAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCt                                            <  1:1125452/51‑1 (MQ=255)
                                              tAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaata                                                 >  1:3366658/1‑45 (MQ=255)
                                              tAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaata                                                 >  1:3958802/1‑45 (MQ=255)
                                              tAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGtt                        <  1:1143506/70‑1 (MQ=255)
                                               aGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTtaata                                                 >  1:4609852/1‑44 (MQ=255)
                                                       ttACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAg              <  1:4464223/71‑1 (MQ=255)
                                                        tACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAg              >  1:1821709/1‑70 (MQ=255)
                                                         aCCCAACTTGTTAGTATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGa            >  1:2497653/1‑71 (MQ=255)
                                                         aCCCAACTTGTTAATATCAATGGTGATGTtaataa                                                >  1:2258649/1‑35 (MQ=255)
                                                         aCCCAACTTGTTAATATCAATGGTGATGTTAATAATcc                                             >  1:789142/1‑38 (MQ=255)
                                                         aCCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                          >  1:1257085/1‑41 (MQ=255)
                                                         aCCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                          >  1:1444622/1‑41 (MQ=255)
                                                         aCCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTcc                              >  1:3138159/1‑53 (MQ=255)
                                                         aCCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGa            >  1:3193115/1‑71 (MQ=255)
                                                          cccAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGa               <  1:4460341/67‑1 (MQ=255)
                                                            cAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAttt         <  1:2263889/71‑1 (MQ=255)
                                                            cAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAttt         <  1:2861289/71‑1 (MQ=255)
                                                             aaCTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                          >  1:1803728/1‑37 (MQ=255)
                                                              aCTTGTTAATATCAATGGTGATGTTAATAATCCTGa                                          <  1:4182542/36‑1 (MQ=255)
                                                               cTTGTTAGTATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt      <  1:2529583/71‑1 (MQ=255)
                                                               cTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt      >  1:2477969/1‑71 (MQ=255)
                                                                ttGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCAtt                           <  1:3116098/49‑1 (MQ=255)
                                                                ttGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt      >  1:2769686/1‑70 (MQ=255)
                                                                ttGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt      <  1:4049992/70‑1 (MQ=255)
                                                                 tGTTAATATCAATGGTGATGTTAATAATCCTGATATTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTa    >  1:3108839/1‑71 (MQ=255)
                                                                 tGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTa    >  1:3199583/1‑71 (MQ=255)
                                                                   ttAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTAt                      >  1:3524891/1‑51 (MQ=255)
                                                                   ttAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTAtt  <  1:2608767/71‑1 (MQ=255)
                                                                      atatCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGAtttttt      <  1:768769/64‑1 (MQ=255)
                                                                      atatCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTAt   >  1:2632851/1‑67 (MQ=255)
                                                                      |                                                                   
ATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATT  >  W3110S.gb/2393269‑2393406

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: ‑