Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,342,193 A→G T95A (ACG→GCG)  melA → alpha‑galactosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,342,1930AG97.4% 118.1 / ‑6.5 39T95A (ACG→GCG) melAalpha‑galactosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (16/22);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGT  >  NC_000913/4342048‑4342329
                                                                                                                                                 |                                                                                                                                         
ggAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCgg                                                                                                                                        <  2:1610/149‑1 (MQ=255)
ggAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTg                                                                                                                                      <  2:287836/151‑1 (MQ=255)
  aaGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGa                                                                                                                                     <  1:495912/150‑1 (MQ=255)
           atatTGTGGTGCGTAAGCTGATGGATTCAGAAGGGGCCAGGGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCg                                                                                                                            <  1:497980/150‑1 (MQ=255)
                     gCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAg                                                                                                                 <  2:617316/151‑1 (MQ=255)
                            cTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTCGGGGTCTGTAAGCgg                                                                                                              >  1:56369/1‑147 (MQ=255)
                              gATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATgg                                                                                                         <  1:143091/150‑1 (MQ=255)
                                               ccAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAAACAtt                                                                                         >  1:24555/1‑149 (MQ=255)
                                                   cGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTTCGCGGGGACTGTTTTCGAGGGCTGTAAGCGGCATGGTCTGGAACAAACCaattcc                                                                                      >  1:446310/1‑143 (MQ=255)
                                                     gCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATAc                                                                                 <  2:410204/151‑1 (MQ=255)
                                                      cAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGGACCTTGCGCGGGGACTGGTTTTGAGGGCTTTAAGCGGCATTGTCTGGGACAAACCATTGCCGATa                                                                                  >  1:424316/1‑149 (MQ=255)
                                                       aaaaTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTTCGCGGGGGCTGGTTTCGAGGGCTGTAAGCGGCATGGGTTGGGACAAACCATTGCCGGTACGt                                                                              >  1:337380/1‑151 (MQ=255)
                                                       aaaaTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGt                                                                              <  2:354919/151‑1 (MQ=255)
                                                        aaaTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGGGGCTGGTTTTGAGGGCTGTAAGCGGCATGGGCTGGAACAAACCATTGCCGCTACGtt                                                                             >  2:198172/1‑151 (MQ=255)
                                                                            cAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgc                                                         <  2:556394/151‑1 (MQ=255)
                                                                            cAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgc                                                         <  2:158397/151‑1 (MQ=255)
                                                                             agaaagaaGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcg                                                        <  1:148752/151‑1 (MQ=255)
                                                                                aagaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGGGACTGATTTCGAGGGCTGTAAGCGGCATGGTCTGGAAAAAACCATTGCCGATACGGTGGGGGCGGGCGGGATTATgcg                                                        >  1:230649/1‑148 (MQ=255)
                                                                                 agaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGGTATGAAACTTGGGGGGGGAATGATTTTGAGGGCTGTAAAGGGCATGGGCTGGAACAAAACAATGCCGGTACGTTGGGGCCGGGCGGTATTATgcgcg                                                      >  1:467164/1‑149 (MQ=255)
                                                                                    agCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTTCACGGTGACTTATTTTGAGGGCTTTAAGCGGCATGGTCTGGAAAAAAACCTTGCCGATACGTTGGGGCCGGGCGGGATTATgcgcgc                                                     >  1:625985/1‑147 (MQ=255)
                                                                                    agCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTa                                                 <  2:432610/151‑1 (MQ=255)
                                                                                         tAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTa                                             <  1:2098/150‑1 (MQ=255)
                                                                                            agGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCCTGCGCGGTGACCGATTTTGAGGGCTGTAAAAGGGAAGGGCTGGGACAAAACATTGCCGATAAGGTTGGGGCGGGGGGGATTATGCGCGCGCTACGTACCa                                          >  2:394818/1‑150 (MQ=255)
                                                                                                     aTTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCTTTCGCGGGGACTGATTTCGAGGTCTGTAAGCGGCAAGGGCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGGATTATTCGCGCGCTACGTACCATTccccctc                                 >  1:483527/1‑145 (MQ=255)
                                                                                                     aTTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCt                                <  2:519039/151‑1 (MQ=255)
                                                                                                       tttGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCTTGCGCGGGGGCTGATTTCGAGGGCTGTAAGGGGCATGGTCTGGAACAAACCATTGCCGGTACGTTTGGGGCCGGGGGGATTAAGCGCGCGCCACCTACCATTCCGcct                                  >  2:244157/1‑145 (MQ=255)
                                                                                                         tGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCTGGCGCGGGGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGGATTATGCGCGCGCTACGTACCATTCCGCATCtgt                              >  2:519816/1‑149 (MQ=255)
                                                                                                          gTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTgg                            <  1:290294/150‑1 (MQ=255)
                                                                                                             gtggtgGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCaaa                        <  2:163455/151‑1 (MQ=255)
                                                                                                                  ggCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTg                    <  2:412358/150‑1 (MQ=255)
                                                                                                                  ggCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTg                    <  1:468421/150‑1 (MQ=255)
                                                                                                                       ttCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAgg               <  2:515752/150‑1 (MQ=255)
                                                                                                                         cAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACa            <  2:24555/151‑1 (MQ=255)
                                                                                                                           gATTGGCGGTTATGAACCTTGCGCGGGGACTGGTTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGGTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAgg               >  2:379159/1‑146 (MQ=255)
                                                                                                                               ggCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACgg      >  2:248782/1‑151 (MQ=255)
                                                                                                                               ggCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTAACATTCCGCATCTGTGGCAAATTTGCGAGGACATGACgg      >  2:518533/1‑151 (MQ=255)
                                                                                                                                gCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGa         <  2:553017/147‑1 (MQ=255)
                                                                                                                                 cGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGa     >  2:123576/1‑150 (MQ=255)
                                                                                                                                   gTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGt  <  1:112947/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
GGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGT  >  NC_000913/4342048‑4342329

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