Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,842,414 G→A *475* (TAG→TAA ascB → cryptic 6‑phospho‑beta‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,842,4140GA97.4% 113.6 / ‑6.5 39*475* (TAG→TAAascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/18);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACAT  >  NC_000913/2842280‑2842546
                                                                                                                                      |                                                                                                                                    
gAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAgtaagtaa                                                                                                                                  <  2:459928/139‑1 (MQ=255)
        cAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCg                                                                                                                          >  2:266277/1‑139 (MQ=255)
                    ctttgtctttgtTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGa                                                                                                              <  2:219282/139‑1 (MQ=255)
                                   ccGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                               >  1:1538779/1‑139 (MQ=255)
                                     gTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                             >  2:1027059/1‑139 (MQ=255)
                                       gacgacGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTATTGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCGTTATCCGGCCTACa                                                                                           >  2:986900/1‑139 (MQ=255)
                                         cgacgCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATg                                                                                         <  1:266277/139‑1 (MQ=255)
                                                  cAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGaa                                                                                >  2:497726/1‑139 (MQ=255)
                                                            cTGACGCGCATGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGc                                                                      <  2:1252081/139‑1 (MQ=255)
                                                                 gcgcACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTAATTGCCAGTAATGGGGAAGGTTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACACCTGATCGGGCCTACATGTCGCACGAACGATGAAGGCCCGAt                                                                 >  1:35345/1‑139 (MQ=255)
                                                                                  cATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCggg                                                <  2:1572919/139‑1 (MQ=255)
                                                                                     tCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATg                                                                                                                      >  2:831792/1‑66 (MQ=255)
                                                                                     tCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATg                                                                                                                      <  1:831792/66‑1 (MQ=255)
                                                                                     tCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCaa                                             >  1:1455296/1‑139 (MQ=255)
                                                                                       tggtggTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAAtt                                           >  1:39647/1‑139 (MQ=255)
                                                                                        ggtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                              <  1:1514785/87‑1 (MQ=255)
                                                                                        ggtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                              >  2:1514785/1‑87 (MQ=255)
                                                                                         gtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGaa                                                                                                             >  2:1177237/1‑71 (MQ=255)
                                                                                         gtggtATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGaa                                                                                                             <  1:1177237/71‑1 (MQ=255)
                                                                                               taAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCAttt                                   >  1:322780/1‑139 (MQ=255)
                                                                                                    aaGTGATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCAt                              <  2:322780/139‑1 (MQ=255)
                                                                                                       tGATTGCCAGTACTGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCag                           >  1:1486319/1‑139 (MQ=255)
                                                                                                       tGATTGCCAGTAATGGGGAAGATTTAGAgtaagtaa                                                                                                                                  >  1:765436/1‑36 (MQ=255)
                                                                                                       tGATTGCCAGTAATGGGGAAGATTTAGAgtaagtaa                                                                                                                                  >  2:1026931/1‑36 (MQ=255)
                                                                                                       tGATTGCCAGTAATGGGGAAGATTTAGAgtaagtaa                                                                                                                                  <  2:765436/36‑1 (MQ=255)
                                                                                                       tGATTGCCAGTAATGGGGAAGATTTAGAgtaagtaa                                                                                                                                  <  1:1026931/36‑1 (MQ=255)
                                                                                                       tGATTGCCAGTAATGGGGAAGATTTAGAgtaagtaa                                                                                                                                  >  2:1346288/1‑36 (MQ=255)
                                                                                                       tGATTGCCAGTAATGGGGAAGATTTAGAgtaagtaa                                                                                                                                  <  1:1346288/36‑1 (MQ=255)
                                                                                                        gATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCa                                                                                    <  2:1277283/81‑1 (MQ=255)
                                                                                                        gATTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCa                                                                                    >  1:1277283/1‑81 (MQ=255)
                                                                                                         aTTGCCAGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCagag                         >  2:520139/1‑139 (MQ=255)
                                                                                                               aGTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGtt                   <  1:497726/139‑1 (MQ=255)
                                                                                                                gTAATGGGGAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTc                  <  1:520139/139‑1 (MQ=255)
                                                                                                                       ggAAGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCAt           <  1:71610/139‑1 (MQ=255)
                                                                                                                         aaGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                   <  2:672495/49‑1 (MQ=255)
                                                                                                                         aaGATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                   >  1:672495/1‑49 (MQ=255)
                                                                                                                           gATTTAGAGTAAGTCACAGTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                     >  2:984013/1‑45 (MQ=18)
                                                                                                                           gATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTAt                                                                                                     <  1:984013/45‑1 (MQ=255)
                                                                                                                           gATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCa                               >  1:1587792/1‑115 (MQ=255)
                                                                                                                           gATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCa                               <  2:1587792/115‑1 (MQ=255)
                                                                                                                           gATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAcc               >  2:124261/1‑131 (MQ=255)
                                                                                                                           gATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCAcc               <  1:124261/131‑1 (MQ=255)
                                                                                                                           gATTTAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAg       <  2:35345/139‑1 (MQ=255)
                                                                                                                             tttAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGaa     <  2:1486319/139‑1 (MQ=255)
                                                                                                                             tttAGAGTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGaa     <  2:1538779/139‑1 (MQ=255)
                                                                                                                                agagTAAGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACTTGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACAt  >  2:851700/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGACGCGCACGCGTAAGAAATCATTCTGGTGGTATAAAAAAGTGATTGCCAGTAATGGGGAAGATTTAGAGTAGGTAACAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACAT  >  NC_000913/2842280‑2842546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: