Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,191,074 A→C 9.6% 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,0740AC 9.6% 200.6 / 15.7 83intergenic (+415/+485)yeaJ/rndpredicted diguanylate cyclase/ribonuclease D
Reads supporting (aligned to +/- strand):  ref base A (46/29);  new base C (3/5);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTA  >  minE/1191009‑1191144
                                                                 |                                                                      
aTGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAgcgc                                                                   <  1:385864/71‑1 (MQ=255)
   tCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTag                                                                >  1:2439360/1‑71 (MQ=255)
   tCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTag                                                                >  1:3836259/1‑71 (MQ=255)
         aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGt                                                          <  1:3838548/71‑1 (MQ=255)
         aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGt                                                          >  1:3144890/1‑71 (MQ=255)
         aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGt                                                          >  1:1128827/1‑71 (MQ=255)
         aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGt                                                          >  1:2823985/1‑71 (MQ=255)
         aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGt                                                          >  1:1188212/1‑71 (MQ=255)
         aGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGt                                                          >  1:2374440/1‑71 (MQ=255)
             tAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAgcg                                                                    >  1:4502783/1‑57 (MQ=255)
             tAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAgcg                                                                    >  1:4607299/1‑57 (MQ=255)
              aaGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTAc                                                     >  1:1468829/1‑71 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACaa                                                                       >  1:2030400/1‑46 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTa                                                                 <  1:1466175/52‑1 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTa                                                                 <  1:3648065/52‑1 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTa                                                                 <  1:1200429/52‑1 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTa                                                                 <  1:4803076/52‑1 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTa                                                                 <  1:2743847/52‑1 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTa                                                                 <  1:390046/52‑1 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGt                                                  >  1:2649070/1‑67 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGa                                                >  1:1627004/1‑69 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGa                                                >  1:3282461/1‑69 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGa                                                >  1:3774901/1‑69 (MQ=255)
                     ggCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGa                                                >  1:3705979/1‑69 (MQ=255)
                       cTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGc                                            <  1:3967781/71‑1 (MQ=255)
                       cTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGc                                            <  1:3952399/71‑1 (MQ=255)
                       cTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGc                                            <  1:1277437/71‑1 (MQ=255)
                       cTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGc                                            <  1:1429474/71‑1 (MQ=255)
                       cTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGc                                            <  1:3707596/71‑1 (MQ=255)
                       cTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGc                                            <  1:2613300/71‑1 (MQ=255)
                       cTGCGGTAAGGCATCATCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGc                                            <  1:2607085/71‑1 (MQ=255)
                            gTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAg                                             <  1:1159833/65‑1 (MQ=255)
                              aaGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACaa                                                                       >  1:2766268/1‑37 (MQ=255)
                              aaGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACaa                                                                       >  1:3016226/1‑37 (MQ=255)
                              aaGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACaa                                                                       >  1:4752435/1‑37 (MQ=255)
                              aaGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACaa                                                                       >  1:4000622/1‑37 (MQ=255)
                                 gCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCt                                  <  1:2908969/71‑1 (MQ=255)
                                    tCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTc                                     >  1:1259917/1‑65 (MQ=255)
                                    tCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTc                                     >  1:2206013/1‑65 (MQ=255)
                                    tCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTc                                     >  1:1266074/1‑65 (MQ=255)
                                    tCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTc                                     >  1:3186529/1‑65 (MQ=255)
                                    tCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTc                                     >  1:1051114/1‑65 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTag                                                                >  1:5132314/1‑36 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCg                              >  1:981970/1‑70 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCg                              >  1:3257066/1‑70 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCg                              >  1:1102630/1‑70 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCg                              >  1:2718790/1‑70 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCg                              >  1:1458252/1‑70 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCg                              >  1:1254692/1‑70 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCg                              >  1:2336976/1‑70 (MQ=255)
                                      ttCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCg                              >  1:2415069/1‑70 (MQ=255)
                                      ttCCGGAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGa                            >  1:4001034/1‑71 (MQ=255)
                                        ccGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGa                                                >  1:4171470/1‑50 (MQ=255)
                                        ccGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGa                            >  1:2861723/1‑70 (MQ=255)
                                        ccGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGa                            >  1:268247/1‑70 (MQ=255)
                                        ccGGCAATGTCTGCGCTTTTTCCACAAGCGCTagc                                                               >  1:4031149/1‑35 (MQ=255)
                                        ccGGCAATGTCTGCGCTTTTTCCACAAGCGCTagc                                                               >  1:1473988/1‑35 (MQ=255)
                                        ccGGCAATGTCTGCGCTTTTTCCACAAGCGCTagc                                                               >  1:3364821/1‑35 (MQ=255)
                                        ccGGCAATGTCTGCGCTTTTTCCACAAGCGCTagc                                                               >  1:3189356/1‑35 (MQ=255)
                                         cGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATa                          <  1:589503/71‑1 (MQ=255)
                                         cGGCAATGTCTGCGCTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATa                          <  1:1582787/71‑1 (MQ=255)
                                          ggCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATa                          >  1:3361730/1‑70 (MQ=255)
                                            cAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAAc                       >  1:1946850/1‑71 (MQ=255)
                                            cAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAAc                       >  1:88019/1‑71 (MQ=255)
                                             aaTGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTg                                                 <  1:3835397/44‑1 (MQ=255)
                                             aaTGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTg                                                 <  1:4981982/44‑1 (MQ=255)
                                             aaTGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTg                                                 <  1:3110853/44‑1 (MQ=255)
                                             aaTGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATaaa                        >  1:1187337/1‑69 (MQ=255)
                                             aaTGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATaaa                        >  1:3333972/1‑69 (MQ=255)
                                             aaTGTCTGCGCTGTTTCCACAAGCGCTAGCAGCGTTTTACCGTg                                                 <  1:3407316/44‑1 (MQ=255)
                                                gTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAg                   >  1:3171586/1‑71 (MQ=255)
                                                gTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAg                   >  1:2164819/1‑71 (MQ=255)
                                                 tCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGc                                   >  1:4199398/1‑54 (MQ=255)
                                                 tCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGc                                   >  1:1611879/1‑54 (MQ=255)
                                                 tCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGc                                   >  1:2532299/1‑54 (MQ=255)
                                                     cgcTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCt                <  1:4176423/69‑1 (MQ=255)
                                                     cgcTTTTTCCACCAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCt                <  1:1849054/69‑1 (MQ=255)
                                                     cgcTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTg               <  1:31917/70‑1 (MQ=255)
                                                     cgcTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTg               <  1:542936/70‑1 (MQ=255)
                                                     cgcTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTCAAAGCGGATTTCGCTACCGGATAAACCCAGGCTg               <  1:4693089/70‑1 (MQ=255)
                                                       cTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCg                                           <  1:2983510/40‑1 (MQ=255)
                                                       cTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCg                                           <  1:774652/40‑1 (MQ=255)
                                                       cTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCg                                           <  1:739165/40‑1 (MQ=255)
                                                       cTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCg                                           <  1:3924470/40‑1 (MQ=255)
                                                       cTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCg                                           <  1:1825618/40‑1 (MQ=255)
                                                       cTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCg                                           <  1:4205098/40‑1 (MQ=255)
                                                          tttCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGt         >  1:5066483/1‑71 (MQ=255)
                                                          tttCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGt         >  1:2343796/1‑71 (MQ=255)
                                                              cacaAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCg      >  1:4803888/1‑70 (MQ=255)
                                                              cacaAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGc     >  1:4401485/1‑71 (MQ=255)
                                                                 aaGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTa  <  1:2212229/71‑1 (MQ=255)
                                                                 |                                                                      
ATGTCCATCAGGTTAAGCATCGGCTGCGGTAAGGCATCTTCCGGCAATGTCTGCGCTTTTTCCACAAGCGCTAGCAGCGTTTTACCGTGAAAGCGGATTTCGCTACCGGATAAACCCAGGCTGTCCAGTTCGCCTA  >  minE/1191009‑1191144

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: