Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,635,354 C→A 7.0% M86I (ATG→ATT)  gpr ← spore germination protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,635,3540CA 7.0% 184.3 / 3.8 57M86I (ATG→ATT) gprspore germination protease
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base A (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

AAGCGCATCAGGTGTCACGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGTGATCACGCTCTTTTTC  >  BS168_reference_genome/2635228‑2635491
                                                                                                                              |                                                                                                                                         
aaGCGCATCAGGTGTCACGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCa                                                                                                                               >  2:248810/1‑139 (MQ=255)
  gcgcATCAGGTGTCACGTTCCAGTTGCCGAGCCCAACGATCACACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGTATAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCAtt                                                                                                                             <  4:348011/139‑1 (MQ=255)
  gcgcATCACGGGGCACGTTCCAGTTGCCGAGCCCAACGATCCAACAGCTCGCATCTTATGAGATATTAAGGTCTTCAAGAACAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCAtt                                                                                                                             <  2:349352/139‑1 (MQ=255)
    gcATCAGGTGTCACGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATtct                                                                                                                           <  4:150085/139‑1 (MQ=255)
       tCAGGTGTCACGTTCCAGTTGCCGAGCCCAACGCTCAAACAGCTCGCAGCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAAGCATTCTccc                                                                                                                        >  1:224359/1‑139 (MQ=255)
       tCAGGTGTCACGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTccc                                                                                                                        >  4:214650/1‑139 (MQ=255)
       tCAGGTGTCACGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTccc                                                                                                                        >  6:46924/1‑139 (MQ=255)
         aGGTGTCGCGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCt                                                                                                                       <  7:263820/138‑1 (MQ=255)
         aGGTGTCACGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTg                                                                                                                      <  8:323994/139‑1 (MQ=255)
          ggTGGCACCTTCCAGATGCCGAGCCCATCGATCAACCAGCTCGCAACTTTTGAGATATTAAGGCTTTCAAGAAAAGGTGAGAATTCCTCCGCAAAGACCGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGa                                                                                                                     <  4:344229/139‑1 (MQ=255)
            tgtCACGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAACGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGAtt                                                                                                                   <  2:334253/139‑1 (MQ=255)
                aCGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCGTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCt                                                                                                               >  5:318077/1‑139 (MQ=255)
                aCGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTccc                                                                                                                        <  4:329532/130‑1 (MQ=255)
                aCGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTccc                                                                                                                        >  3:329532/1‑130 (MQ=255)
                aCGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTccc                                                                                                                <  3:214650/138‑1 (MQ=255)
                aCGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCt                                                                                                               >  5:83029/1‑139 (MQ=255)
                aCGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCt                                                                                                               >  1:287363/1‑139 (MQ=255)
                aCGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCt                                                                                                               >  7:107338/1‑139 (MQ=255)
                       aGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTc                                                                                                        <  5:46924/139‑1 (MQ=255)
                       aGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTc                                                                                                        <  1:265097/139‑1 (MQ=255)
                        gTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCAtt                                                                                                                             >  2:230673/1‑117 (MQ=255)
                        gTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCAtt                                                                                                                             <  1:230673/117‑1 (MQ=255)
                        gTTGCCGAGACCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTTCTTCa                                                                                                       <  7:298485/139‑1 (MQ=255)
                                 cccAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTaaa                                                                                              >  1:144165/1‑139 (MQ=255)
                                      cGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAAt                                                                                                                                 >  8:232829/1‑99 (MQ=255)
                                      cGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAAt                                                                                                                                 <  7:232829/99‑1 (MQ=255)
                                                 cTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTc                                                                                                        <  6:26039/113‑1 (MQ=255)
                                                 cTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTc                                                                                                        >  5:26039/1‑113 (MQ=255)
                                                           ttGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGa                                                                       <  6:318077/136‑1 (MQ=255)
                                                           ttGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGa                                                                       <  2:287363/136‑1 (MQ=255)
                                                           ttGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGa                                                                       <  6:83029/136‑1 (MQ=255)
                                                              agaTATTCAGGTTTTCAAGAAAAGCTGAGAATTGCTCCGCAAATACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTACCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGct                                                                  <  6:57717/138‑1 (MQ=255)
                                                                           ttCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCttt                                                    >  5:20419/1‑139 (MQ=255)
                                                                           ttCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCttt                                                    >  8:20152/1‑139 (MQ=255)
                                                                           ttCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGAATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCttt                                                    >  7:47132/1‑139 (MQ=255)
                                                                           ttCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGAATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGATCTGCGCCTTCCttt                                                    >  7:52741/1‑139 (MQ=255)
                                                                              aaGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGt                                                  <  8:107338/138‑1 (MQ=255)
                                                                                  aaaGCTGAGAATTCCTCCGCAAAGACAGCAGAGAGCTTCTCCTGCATTTCAGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGAGACCTTCGTTTTTTACCGAAAGGTCTGCGCCTTCCTGTGTaaaa                                              >  1:174057/1‑136 (MQ=255)
                                                                                                ctccGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTg                                                                                                                      <  2:86354/52‑1 (MQ=255)
                                                                                                ctccGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTg                                                                                                                      >  1:86354/1‑52 (MQ=255)
                                                                                                ctccGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGcct                                                          >  4:71104/1‑112 (MQ=255)
                                                                                                ctccGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGcct                                                          <  3:71104/112‑1 (MQ=255)
                                                                                                ctccGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGAtt                               >  3:232027/1‑139 (MQ=255)
                                                                                                     cAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACAt                                             <  2:106908/120‑1 (MQ=255)
                                                                                                     cAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACAt                                             >  1:106908/1‑120 (MQ=255)
                                                                                                            atcagAGACCTTCTCCTGCATTTACGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCACGAAAGCTCTGCGCCTTCCTTTTAAACATCAACTGGGGTTATTTTAAAGccgccg                   <  2:174057/137‑1 (MQ=255)
                                                                                                             gcagAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGccgccgg                  >  1:179222/1‑138 (MQ=255)
                                                                                                             gcagAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGt                  >  2:141062/1‑139 (MQ=255)
                                                                                                                  gaCCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGTGATCa             <  4:337915/139‑1 (MQ=255)
                                                                                                                      ttCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCttt                                                    >  1:246180/1‑96 (MQ=255)
                                                                                                                      ttCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCttt                                                    <  2:246180/96‑1 (MQ=255)
                                                                                                                        ctcCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGaaa                                                                     >  7:15924/1‑77 (MQ=255)
                                                                                                                        ctcCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGaaa                                                                     <  8:15924/77‑1 (MQ=255)
                                                                                                                        ctcCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGTGATCACGctct       >  3:118713/1‑139 (MQ=255)
                                                                                                                        ctcCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGGGATCAGGctcg       >  6:19694/1‑138 (MQ=255)
                                                                                                                           cTGAATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGTGATCACGCTCttt     <  8:47132/138‑1 (MQ=255)
                                                                                                                           cTGAATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGTGATCACGCTCttt     <  8:52741/138‑1 (MQ=255)
                                                                                                                              cATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGTGATCACGCTCTTTTTc  >  4:366654/1‑138 (MQ=255)
                                                                                                                              |                                                                                                                                         
AAGCGCATCAGGTGTCACGTTCCAGTTGCCGAGCCCAACGATCAAACAGCTCGCATCTTTTGAGATATTAAGGTTTTCAAGAAAAGCTGAGAATTCCTCCGCAAAGACAGCAGAGACCTTCTCCTGCATTTCCGAATCATTCTCCCTGATTCCCTGTGCTTCAAGCGTTAAATAGCGCCCTTCTTTTTTTCCCGAAAGCTCTGCGCCTTCCTTTGTAACATCAACTGTGCGGATTTTAATGCCGCCGTGATCACGCTCTTTTTC  >  BS168_reference_genome/2635228‑2635491

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑