Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,903,168 G→A H178H (CAC→CAT)  bglB ← 6‑phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,903,1680GA92.3% 69.2 / ‑2.4 26H178H (CAC→CAT) bglB6‑phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (9/15);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCA  >  NC_000913/3903034‑3903304
                                                                                                                                      |                                                                                                                                        
ccGATTTTGGCTTCGGGCAGCAGGCTATGACAGCCGTTATCTGCCCGCACACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTaa                                                                                                                                     <  2:309847/140‑1 (MQ=255)
                        cTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGt                                                                                                             <  2:55246/140‑1 (MQ=255)
                         tATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTc                                                                                                            <  2:551743/140‑1 (MQ=255)
                                 gCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATaa                                                                                                    >  1:609236/1‑140 (MQ=255)
                                     tAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGcc                                                                                                <  1:152001/140‑1 (MQ=255)
                                           cccGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTa                                                                                          <  1:12413/140‑1 (MQ=255)
                                           cccGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTa                                                                                          >  2:12413/1‑140 (MQ=255)
                                           cccGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGATGATTTCATTAACCGTCGGCCATAACGCCACTTTa                                                                                          >  2:223130/1‑140 (MQ=255)
                                           cccGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGACACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTa                                                                                          >  2:415799/1‑140 (MQ=255)
                                                  aCTGGCAACCAGTTGATAGTGGTTAGCCTGATAAATTTCTGCCTCGCCTCTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGt                                                                                   <  2:249038/140‑1 (MQ=255)
                                                    tGGCAACCAGTTGATGGTGGATCGCCTGATATACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAg                                                                                 <  2:191100/140‑1 (MQ=255)
                                                    tGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAg                                                                                 <  1:96522/140‑1 (MQ=255)
                                                    tGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAg                                                                                 <  2:493671/140‑1 (MQ=255)
                                                    tGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAg                                                                                 <  2:411333/140‑1 (MQ=255)
                                                    tGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAg                                                                                 <  1:117773/140‑1 (MQ=255)
                                                    tGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAg                                                                                 <  1:640783/140‑1 (MQ=255)
                                                          ccAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTa                                                                           <  1:415799/140‑1 (MQ=255)
                                                              ttGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGa                                                                       <  1:37219/140‑1 (MQ=255)
                                                              ttGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGa                                                                       <  2:587113/140‑1 (MQ=255)
                                                                gATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGAcc                                                                     >  2:70382/1‑140 (MQ=255)
                                                                     tggATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCg                                                                >  1:19452/1‑140 (MQ=255)
                                                                                     cTTCCGCCCCGCCACTCCCTTCAGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAg                                                <  2:297469/140‑1 (MQ=255)
                                                                                           ccTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGGCCGTGCCGGCGTAGTGCTCGAAGTGATCg                                          >  1:297469/1‑140 (MQ=255)
                                                                                                     ctctTCTGCCAGCCCCACGCCCGTGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTc                                >  2:98244/1‑140 (MQ=255)
                                                                                                                            tGAATGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGttt         >  2:655707/1‑140 (MQ=255)
                                                                                                                              aaTGGCGCATGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATACTCGAAGTAATCGATGGCCGCTCGAGTAGCCCAACCGCCGTAGttttt       >  2:198051/1‑140 (MQ=255)
                                                                                                                                   cgcaTGTAACGAGATGTTGTTTTCATCACAGGTCACCCATAGCCCCTTTTTATGTTGGGAGGGAGTAAAGCCCGGGGGGGTTTAAGGCTCGAAGTGCGCGTTGTCCCCTGGATTCCCCCAACCGCCGTAGTTTTTCACCa  <  2:609236/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                        
CCGATTTTCGCTTCGGGGAGCAGGCTATGACAGGCTTTAACTGCCCGCGCACTGGCAACCAGTTGATGGTGGATAGCCTGATAAACTTCCGCCTCGCCACTCTCTTCTGCCAGCCCCACGCCCGTGAATGGCGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCA  >  NC_000913/3903034‑3903304

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