Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,193,779 G→T L143I (CTT→ATT)  yeiB ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,193,7790GT96.8% 83.0 / NA 31L143I (CTT→ATT) yeiBconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (10/20);  minor base A (1/0);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GCCGATGCATCCGGCGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAAGCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCATCGCGCACCAGCC  >  REL606/2193706‑2193847
                                                                         |                                                                    
gCCGATGCATCCGGCGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATcc                                                                    <  2:1162901/76‑1 (MQ=255)
    aTGCATCCGGCGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCccac                                                                >  1:320274/1‑76 (MQ=255)
     tGCATCCGGCGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAAAGCCAATCCccacc                                                               >  1:243209/1‑76 (MQ=255)
       cATCCGGCGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAg                                                             >  2:748840/1‑76 (MQ=255)
        aTCCGGCGTCCAAGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCGCCAGa                                                            >  2:51641/1‑76 (MQ=255)
            ggCGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATaaa                                                        <  2:984414/76‑1 (MQ=255)
             gCGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAg                                                       >  1:263551/1‑76 (MQ=255)
              cGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGc                                                      >  2:411262/1‑76 (MQ=255)
              cGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGc                                                      <  1:732673/76‑1 (MQ=255)
              cGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGc                                                      <  2:66741/76‑1 (MQ=255)
                  cAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGa                                                  <  2:410211/76‑1 (MQ=255)
                          gCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTa                                          <  1:821318/76‑1 (MQ=255)
                           cTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCTTGTAt                                         >  2:55643/1‑76 (MQ=255)
                            tGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTAtt                                        >  2:1171874/1‑76 (MQ=255)
                                 tGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACa                                   <  2:1162391/76‑1 (MQ=255)
                                   gCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAgg                                 <  1:595127/76‑1 (MQ=255)
                                      gTCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCtt                              <  2:102839/76‑1 (MQ=255)
                                       tCGGAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCttt                             >  1:450904/1‑76 (MQ=255)
                                         ggAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCtnttt                           <  1:1206572/76‑1 (MQ=255)
                                          gAAATCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTa                          <  1:430647/76‑1 (MQ=255)
                                            aaTCAGCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTAcc                        <  1:659104/76‑1 (MQ=255)
                                                 gCCCCAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATgg                   <  2:400412/76‑1 (MQ=255)
                                                   cccAATAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGgcg                 <  2:766094/76‑1 (MQ=255)
                                                      aaTAACANCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCAt              >  2:1075296/1‑76 (MQ=255)
                                                      aaTAACAGNAAATCGCCAATCCCCGCCAGATAAAGGATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCAt              <  2:231403/76‑1 (MQ=255)
                                                      aaTAACAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCAt              >  1:187503/1‑76 (MQ=255)
                                                          aCAGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCATcgcg          <  2:750200/76‑1 (MQ=255)
                                                           cAGCAAAACGCCAAACCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCATCgcgc         >  1:56813/1‑76 (MQ=255)
                                                            aGCAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCATCGCGCa        <  1:475379/76‑1 (MQ=255)
                                                              cAAAACGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCATCGCGCAcc      <  2:419129/76‑1 (MQ=255)
                                                               aaaaCGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCATCGCGCACCa     <  1:802136/76‑1 (MQ=255)
                                                                  aCGCCAATCCCCACCAGATAAAGCATGACGCCTGTATTNAAGAGGCTTTTTACCGATGGCGCATCGCGCACCAgcc  <  1:1307478/76‑1 (MQ=255)
                                                                         |                                                                    
GCCGATGCATCCGGCGTCCAGGCGCGGCTGGTCTGGCTGTCGGAAATCAGCCCCAATAACAGCAAAACGCCAAGCCCCACCAGATAAAGCATGACGCCTGTATTAAACAGGCTTTTTACCGATGGCGCATCGCGCACCAGCC  >  REL606/2193706‑2193847

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: