Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,652,544 T→G 13.1% D85A (GAC→GCC)  nucB ← sporulation‑specific matrix degrading exported DNase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,652,5440TG13.1% 82.8 / 2.2 36D85A (GAC→GCC) nucBsporulation‑specific matrix degrading exported DNase
Reads supporting (aligned to +/- strand):  ref base T (15/16);  new base G (3/2);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGTACCGTCAGGATAGCTGCTCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATCGCATCCCTAATATGACTGCCGGTTTCC  >  BS168_reference_genome/2652408‑2652672
                                                                                                                                        |                                                                                                                                
tctGGTACCGTCAGGATAGCTGCTCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGGCCCATATCGGCCATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCCTCGGCCCCTCATCCCGGGCa                                                                                                                                >  6:194380/1‑139 (MQ=255)
         gTCAGGATAGCTGCTCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGc                                                                                                                       <  5:194380/139‑1 (MQ=255)
          tCAGGATAGCTGCTCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGcc                                                                                                                           <  8:218874/134‑1 (MQ=255)
          tCAGGATAGCTGCTCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGcc                                                                                                                           >  7:218874/1‑134 (MQ=255)
          tCAGGATAGCTGCTCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCt                                                                                                                      >  4:128526/1‑139 (MQ=255)
          tCAGGATAGCTGCTCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCt                                                                                                                      >  4:248186/1‑139 (MQ=255)
                   ctgctCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCCCTCATCCCGGTCATAGCCCGGCTTGGTCGGGa                                                                                                             >  7:68231/1‑139 (MQ=255)
                   ctgctCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGa                                                                                                             >  6:234557/1‑139 (MQ=255)
                      ctCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGAt                                                                                                            <  7:27793/137‑1 (MQ=255)
                      ctCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGGCATAGCCCGGCTTGGGCGGGAt                                                                                                            <  3:134228/137‑1 (MQ=255)
                      ctCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACAAATCTGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCCTCGGCCACACATCCCCGTCATAGCCCGCCGTGGTCGGGCTTc                                                                                                          >  3:62268/1‑139 (MQ=255)
                        cATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCCGAAGGCGTCACCTATCGGCCATCAGCCCCTGCCCCGCCGTCCTCGCCGACCGCCCTCGGCCACTCATCCCGGGCGTAGCCCGGCTTGGGCggggtggcc                                                                                                        >  6:173278/1‑133 (MQ=255)
                        cATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTccc                                                                                                        >  3:148474/1‑139 (MQ=255)
                        cATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCCGCATGGTCGGGAATcc                                                                                                         >  4:61292/1‑138 (MQ=255)
                        cATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCAGCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTccc                                                                                                        >  2:222620/1‑139 (MQ=255)
                                tCCCTGCCCCCGAGCCGGCGCCGCGATTTTCAGAAGGCGTCACATATCGGACATCAGCCCCTGCCCCGCCGTCCTCGCCGCCCGCCCTCGGCCACTCATCCCGGGCCTCTCCCGGCGGGGGCGGGATTCCCGTTaaa                                                                                                  >  2:74816/1‑137 (MQ=255)
                                     aCCCACGAGCCGGCGCCGCGATTATCAGAAGCCGTCACATATCCGACATCAGCCCCTGCACCGCCTTACTCGCAGACCGCCATCGGCCACTCAACCCCGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATtcct                                                                                           >  8:8662/1‑139 (MQ=255)
                                       ccACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATtcctcc                                                                                         <  5:234557/139‑1 (MQ=255)
                                           gAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGcc                                                                                     >  4:245895/1‑139 (MQ=255)
                                                ggCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTg                                                                                >  8:212404/1‑139 (MQ=255)
                                                 gcgcCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGt                                                                               <  4:309134/139‑1 (MQ=255)
                                                            aTCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCAtc                                                                     <  6:149436/138‑1 (MQ=255)
                                                            aTCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCAtc                                                                     <  7:114449/138‑1 (MQ=255)
                                                            aTCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCAtc                                                                     <  3:11722/138‑1 (MQ=255)
                                                             tCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCAtct                                                                    <  7:212404/138‑1 (MQ=255)
                                                             tCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCAtct                                                                    <  8:211437/138‑1 (MQ=255)
                                                               aGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATctct                                                                  <  4:148474/138‑1 (MQ=255)
                                                               aGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTg                                                                 <  7:12841/139‑1 (MQ=255)
                                                                                                   gCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATc                             >  6:242427/1‑139 (MQ=255)
                                                                                                    ccTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATCg                            <  3:128526/139‑1 (MQ=255)
                                                                                                                aCCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATCGCATCCCTAatat                >  3:122926/1‑139 (MQ=255)
                                                                                                                  cGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATCGCATCCCTAATATGa              <  3:248186/139‑1 (MQ=255)
                                                                                                                  cGCCATCGGCCACTCATCCCGGGCCTAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGGCCCTCTGCAATCGCAGCCCTACTCTGa              <  1:74816/139‑1 (MQ=255)
                                                                                                                       tCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATCGCATCCCTAATATGACTGcc         <  3:245895/139‑1 (MQ=255)
                                                                                                                             aCTCATCCCTGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATCGCATCCCTAATATGACTGCCGGTTTc   <  5:245881/139‑1 (MQ=255)
                                                                                                                              cTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATCGCATCCCTAATATGACTGCCGGTTTcc  >  4:126863/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                
TCTGGTACCGTCAGGATAGCTGCTCATTTGATTCCCTACCCACGAGCCGGCGCCGCGATTATCAGAAGGCGTCACATATCGGACATCAGCCCCTGCACCGCCTTCCTCGCAGACCGCCATCGGCCACTCATCCCGGTCATAGCCCGGCTTGGTCGGGATTCCCTTTAAAGATTCCTCCCGCCTTTTGTCTGCTCCATCTCTGTCAATGGTACAAATATCTGGATGTCCCTCTGCAATCGCATCCCTAATATGACTGCCGGTTTCC  >  BS168_reference_genome/2652408‑2652672

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: ‑