Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,088,438 A→T 7.2% Y109F (TAC→TTC)  yhaI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,088,4380AT 7.2% 197.5 / 6.3 69Y109F (TAC→TTC) yhaIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAA  >  Exported/2088373‑2088507
                                                                 |                                                                     
gTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGa                                                                   >  1:1751347/1‑70 (MQ=255)
 tttATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGaaa                                                                 <  1:2800223/71‑1 (MQ=255)
   tATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAAc                                                                >  1:656025/1‑70 (MQ=255)
    aTTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACg                                                               <  1:1183070/70‑1 (MQ=255)
    aTTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACg                                                               <  1:908424/70‑1 (MQ=255)
     ttGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGa                                                              <  1:1269827/70‑1 (MQ=255)
      tggctggcTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTa                                                                       <  1:711194/60‑1 (MQ=255)
        gctggctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACg                                                                     >  1:2960586/1‑60 (MQ=255)
          tggctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCg                                                          >  1:1513797/1‑69 (MQ=255)
            gctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAg                                                       <  1:2444319/70‑1 (MQ=255)
            gctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGt                                                      >  1:940797/1‑71 (MQ=255)
                gTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGttt                                                    <  1:3322472/69‑1 (MQ=255)
                gTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTgg                                                  >  1:2094172/1‑71 (MQ=255)
                gTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTgg                                                  >  1:2032576/1‑71 (MQ=255)
                     tCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCa                                              >  1:948818/1‑70 (MQ=255)
                      cTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCaaa                                            >  1:998736/1‑71 (MQ=255)
                        ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGa                                                              <  1:2831140/51‑1 (MQ=255)
                        ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGa                                                              <  1:877461/51‑1 (MQ=255)
                        ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGa                                                              <  1:747196/51‑1 (MQ=255)
                        ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAAt                                           >  1:133507/1‑70 (MQ=255)
                         gTTTTTTTCTGCACAGAAGGTACTTCTTGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAAtt                                          <  1:991033/70‑1 (MQ=255)
                         gTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTa                                         <  1:3614027/71‑1 (MQ=255)
                         gTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTa                                         <  1:493242/71‑1 (MQ=255)
                          tttttttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTaa                                        <  1:1883243/71‑1 (MQ=255)
                            tttttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTAc                                                                      <  1:1476026/39‑1 (MQ=255)
                            tttttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAAt                                       <  1:1661962/70‑1 (MQ=255)
                            tttttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATc                                      >  1:3649611/1‑71 (MQ=255)
                                 cTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAg                                  >  1:1537249/1‑70 (MQ=255)
                                    cacaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAAt                               <  1:2860554/70‑1 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGaaa                                                                 >  1:2225127/1‑35 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGaaa                                                                 >  1:2387553/1‑35 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGaaa                                                                 >  1:1454653/1‑35 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTa                                   <  1:2417027/65‑1 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTa                                   <  1:2694569/65‑1 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAAt                               <  1:1860574/69‑1 (MQ=255)
                                       aGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTa                                   <  1:3518529/63‑1 (MQ=255)
                                       aGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAAt                               >  1:181967/1‑67 (MQ=255)
                                       aGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTg                             >  1:285557/1‑69 (MQ=255)
                                       aGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTAGGTTCAAATTAATCTTAGAATTg                             >  1:382726/1‑69 (MQ=255)
                                        gAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTggg                           >  1:943477/1‑70 (MQ=255)
                                         aaGGTACTTCTGGCAGCAATCGTTTCCGTAACGAt                                                             >  1:554743/1‑35 (MQ=21)
                                         aaGGTACTTCTGGCAGCAATCGTTTCCGTAACGAt                                                             >  1:3542184/1‑35 (MQ=21)
                                         aaGGTACTTCTGGCAGCAATCGTTTCCGTAACGAt                                                             >  1:3636099/1‑35 (MQ=21)
                                         aaGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGc                         >  1:3141962/1‑71 (MQ=255)
                                            gTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGc                         <  1:1259329/68‑1 (MQ=255)
                                               cTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCg                        <  1:3157854/66‑1 (MQ=255)
                                               cTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCg                        <  1:682525/66‑1 (MQ=255)
                                               cTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATAtt                   <  1:2993969/71‑1 (MQ=255)
                                               cTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATAtt                   <  1:1368469/71‑1 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTg                                                   >  1:2793462/1‑36 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                      >  1:3538681/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                      >  1:2808577/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                      >  1:2721553/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                      >  1:2703910/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                      >  1:1624977/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                      >  1:2306034/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTc                 <  1:1741141/70‑1 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCg                >  1:2929257/1‑71 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCg                >  1:695694/1‑71 (MQ=255)
                                                       agcaATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCt           <  1:2947200/71‑1 (MQ=255)
                                                       agcaATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCt           <  1:1493511/71‑1 (MQ=255)
                                                        cggaTCGTTTCGGTAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCt           <  1:2905270/67‑1 (MQ=255)
                                                          gttCGTTTCCGTAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCt           <  1:1861340/66‑1 (MQ=255)
                                                          aaTCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCtttt        >  1:207655/1‑71 (MQ=255)
                                                          aaTCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCtttt        <  1:2320696/71‑1 (MQ=255)
                                                             cGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAg                                  <  1:2015858/42‑1 (MQ=255)
                                                             cGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAg                                  <  1:2011637/42‑1 (MQ=255)
                                                              gTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTttatt    <  1:2513754/71‑1 (MQ=255)
                                                               ttACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAAtt                              <  1:2604305/44‑1 (MQ=255)
                                                               ttACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTttatt    >  1:1332242/1‑70 (MQ=255)
                                                                 aCGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTaa  >  1:710161/1‑70 (MQ=255)
                                                                 |                                                                     
GTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAA  >  Exported/2088373‑2088507

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: