Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 798,591 A→T 11.1% I88N (ATT→AAT)  ybhA ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb798,5910AT11.1% 113.6 / 7.3 45I88N (ATT→AAT) ybhApredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (17/23);  new base T (2/3);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

ATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAATAGGTGCCATTACAGCA  >  W3110S.gb/798519‑798661
                                                                        |                                                                      
aTTGCATCATCGACATACATCAGACCGTGAATGTG‑‑GTTCATTCAGCATCTCAATGAGTTGCAGGGCTTtaa                                                                        >  1:358010/1‑71 (MQ=255)
    catcatCGACATACATCAGACCGTGAATGTGGTGTTCATTCGGCATCTCAATGAGTTGCAGGGCTTtaat                                                                       <  1:1070624/70‑1 (MQ=255)
          cGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACggg                                                                 <  1:3365046/70‑1 (MQ=255)
          cGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACgg                                                                  >  1:2186194/1‑69 (MQ=255)
            acatacatCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTtaat                                                                       <  1:3390749/62‑1 (MQ=255)
                 catCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGt                                                         >  1:2652800/1‑71 (MQ=255)
                 catCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGt                                                         <  1:1928466/71‑1 (MQ=255)
                 catCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGt                                                         >  1:1922168/1‑71 (MQ=255)
                 catCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGt                                                         >  1:1545086/1‑71 (MQ=255)
                   tCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTcc                                                       >  1:2432910/1‑71 (MQ=255)
                   tCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTcc                                                       <  1:328674/71‑1 (MQ=255)
                    cAGACCGTGAATGTGGTGTTCATACAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGt                                                         >  1:3673737/1‑68 (MQ=255)
                     aGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGc                                                     >  1:1328599/1‑71 (MQ=255)
                      gACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCt                                                    <  1:2450329/71‑1 (MQ=255)
                            gAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTATTAACGGGCATTGGGTCCGCTTCCAg                                               >  1:364282/1‑70 (MQ=255)
                            gAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTATTAACGGGCATTGGGTCCGCTTCCAg                                               >  1:2404688/1‑70 (MQ=255)
                             aaTGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAg                                               <  1:533983/69‑1 (MQ=255)
                             aaTGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGc                                              <  1:1818462/70‑1 (MQ=255)
                             aaTGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGc                                              <  1:3844521/70‑1 (MQ=255)
                               tgtgGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACg                                           >  1:2327319/1‑71 (MQ=255)
                                 tggtgTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGt                                         >  1:570318/1‑71 (MQ=255)
                                    tgtTCATTCAGCATCTCAATGAGTTGCAGGGCTTtaataa                                                                     <  1:2195737/40‑1 (MQ=255)
                                      ttcattcaGCATCTCAATGAGTTGCAGGGCTTtaa                                                                        <  1:2853626/35‑1 (MQ=255)
                                      ttcattcaGCATCTCAATGAGTTGCAGGGCTTTATTAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTg                                    <  1:754510/71‑1 (MQ=255)
                                      ttcattcaGCATCTCAATGAGTTGCAGGGCTTTATTAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTg                                    <  1:2438173/71‑1 (MQ=255)
                                      ttaattcaGCATCTCAATGAGTTGCAGGGCTTTATTAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTg                                    <  1:1217155/68‑1 (MQ=255)
                                         attcaGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCa                                                <  1:1705136/56‑1 (MQ=255)
                                            caGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGa                               >  1:3307370/1‑70 (MQ=255)
                                                   tcAATGAGTTGCAGGGCTTTAATAACGGGCATTgg                                                           >  1:2730436/1‑35 (MQ=255)
                                                   tcAATGAGTTGCAGGGCTTTAATAACGGGCATTgg                                                           >  1:2066335/1‑35 (MQ=255)
                                                   tcAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATa                       >  1:217803/1‑71 (MQ=255)
                                                   tcAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATa                       >  1:3409973/1‑71 (MQ=255)
                                                     aaTGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATTATAATCATACa                     <  1:1742625/71‑1 (MQ=255)
                                                     aaTGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCa                         <  1:3687988/67‑1 (MQ=255)
                                                     aaTGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACa                     <  1:3024065/71‑1 (MQ=255)
                                                     aaTGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACa                     <  1:3788964/71‑1 (MQ=255)
                                                       tGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACaaa                   <  1:1003801/71‑1 (MQ=255)
                                                       tGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACaaa                   <  1:3292195/71‑1 (MQ=255)
                                                         aGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCa                                             >  1:3263668/1‑43 (MQ=255)
                                                         aGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAAt                  <  1:465818/70‑1 (MQ=255)
                                                         aGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAAt                  <  1:713363/70‑1 (MQ=255)
                                                         aGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAAt                  <  1:2607744/70‑1 (MQ=255)
                                                          gTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAATAg                <  1:2136577/71‑1 (MQ=255)
                                                          gTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAATAg                <  1:2447467/71‑1 (MQ=255)
                                                               aGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAATAGGTGc            >  1:3294465/1‑70 (MQ=255)
                                                               aGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAATAGGTGc            >  1:544089/1‑70 (MQ=255)
                                                                 ggCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACa                     <  1:3424989/59‑1 (MQ=255)
                                                                 ggCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACa                     <  1:3806592/59‑1 (MQ=255)
                                                                 ggCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACa                     <  1:526080/59‑1 (MQ=255)
                                                                   cTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAATAGGTGCCATTa       <  1:1297050/71‑1 (MQ=255)
                                                                       aataaCGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATAc                      >  1:3733740/1‑52 (MQ=255)
                                                                        ataaCGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAATAGGTGCCATTACAGCa  <  1:3565948/71‑1 (MQ=255)
                                                                        |                                                                      
ATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCATACAAATAGGTGCCATTACAGCA  >  W3110S.gb/798519‑798661

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: