Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,191,074 A→C 55.8% intergenic (+415/+485) yeaJ → / ← rnd predicted diguanylate cyclase/ribonuclease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,191,0740AC55.8% 40.8 / 50.9 52intergenic (+415/+485)yeaJ/rndpredicted diguanylate cyclase/ribonuclease D
Reads supporting (aligned to +/- strand):  ref base A (9/14);  new base C (24/5);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 1.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

TGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCG  >  minE/1191010‑1191140
                                                                |                                                                  
tGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAgcgc                                                               >  1:1002260/1‑70 (MQ=255)
 gTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCCCAAgcg                                                                >  1:2719077/1‑68 (MQ=255)
  tCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTag                                                            <  1:2939628/71‑1 (MQ=255)
   ccATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTag                                                            >  1:2830873/1‑70 (MQ=255)
     aTCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTagcag                                                         >  1:2606585/1‑71 (MQ=255)
        aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGt                                                      <  1:229493/71‑1 (MQ=255)
        aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCg                                                       <  1:58912/70‑1 (MQ=255)
          gTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAgc                                                                 >  1:2351800/1‑58 (MQ=255)
          gTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAgc                                                                 <  1:1707042/58‑1 (MQ=255)
          gTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAgc                                                                 <  1:1702516/58‑1 (MQ=255)
          gTTAAGCATCGGATGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGttt                                                    <  1:3093434/71‑1 (MQ=255)
            tAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTa                                                  >  1:335517/1‑71 (MQ=255)
            tAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTa                                                  >  1:3210262/1‑71 (MQ=255)
            tAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTa                                                  >  1:2894871/1‑71 (MQ=255)
            tAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTa                                                  >  1:76937/1‑71 (MQ=255)
            tAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTa                                                  >  1:2505196/1‑71 (MQ=255)
             aaGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTAc                                                 >  1:376754/1‑71 (MQ=255)
              aGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTAc                                                 <  1:136199/70‑1 (MQ=255)
                cATCGGCTGCGGTAAGCCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGt                                              <  1:3406513/71‑1 (MQ=255)
                 aTCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGt                                              >  1:1568542/1‑70 (MQ=255)
                    ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGa                                            >  1:1608846/1‑69 (MQ=255)
                    ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGaaa                                          <  1:776269/71‑1 (MQ=255)
                      cTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAg                                         >  1:3541574/1‑70 (MQ=255)
                           gTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGAtt                                   <  1:1012947/71‑1 (MQ=255)
                             aaGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTc                                 >  1:628412/1‑71 (MQ=255)
                              aGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTc                                 >  1:1089398/1‑70 (MQ=255)
                                 cATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTa                             >  1:2649160/1‑71 (MQ=255)
                                  aTCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTAc                            >  1:395485/1‑71 (MQ=255)
                                  aTCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTAc                            >  1:1133451/1‑71 (MQ=255)
                                    cTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTAcc                           >  1:262735/1‑70 (MQ=255)
                                       ccGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCgg                         <  1:2454615/69‑1 (MQ=255)
                                           cAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAAc                   >  1:1203028/1‑71 (MQ=255)
                                            aaTGTCTGCGCTTTTTCCACCAGCGCTAGGAGCGTTTTACCGTGAAAGCGGATTTCGc                               >  1:2451093/1‑58 (MQ=255)
                                             aTGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAAcc                  >  1:1198065/1‑70 (MQ=255)
                                             aTGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAAcc                  >  1:1443649/1‑70 (MQ=255)
                                             aTGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGAt                                    <  1:653230/52‑1 (MQ=255)
                                               gTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCg                                       >  1:1042581/1‑47 (MQ=255)
                                               gTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACACAg               >  1:10988/1‑71 (MQ=255)
                                                tCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTAc                            >  1:1857481/1‑57 (MQ=255)
                                                    cgcTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGc             <  1:3134546/68‑1 (MQ=255)
                                                    cgcTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGc             <  1:989355/68‑1 (MQ=255)
                                                    cgcTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGc             <  1:3279188/68‑1 (MQ=255)
                                                     gcTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCAAGGCTGTc         >  1:3083148/1‑71 (MQ=255)
                                                     gcTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTc         >  1:2622607/1‑71 (MQ=255)
                                                       tttttCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTAc                            <  1:793299/50‑1 (MQ=255)
                                                         tttCCACCAGCGCTAGCAGCGTTTTACCGTGaaat                                         >  1:3849953/1‑34 (MQ=37)
                                                         tttCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTAc                            >  1:2770867/1‑48 (MQ=255)
                                                         tttCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAg      >  1:81108/1‑70 (MQ=255)
                                                         tttCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGt     >  1:118816/1‑71 (MQ=255)
                                                          ttCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCt                              <  1:93106/45‑1 (MQ=255)
                                                          ttCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGt     <  1:1991591/70‑1 (MQ=255)
                                                            ccACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTc   <  1:2026585/70‑1 (MQ=255)
                                                             cacaAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCg  <  1:1524034/70‑1 (MQ=255)
                                                                |                                                                  
TGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCG  >  minE/1191010‑1191140

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: