Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,547,986 C→T *695* (TAG→TAA)  malQ ← 4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,547,9860CT97.0% 100.3 / NA 33*695* (TAG→TAA) malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/14);  minor base A (1/0);  total (19/14)
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.98e-01

TGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGC  >  NC_000913/3547853‑3548119
                                                                                                                                     |                                                                                                                                     
tgctgctgAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActt                                                                                                                                  >  1:208241/1‑139 (MQ=255)
      tgAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtcttctt                                                                                                                            <  1:346919/139‑1 (MQ=255)
      tgAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtcttctt                                                                                                                            <  2:59572/139‑1 (MQ=255)
        aaTATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCg                                                                                                                          >  1:595388/1‑139 (MQ=255)
               tGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGc                                                                                                                   <  2:208241/139‑1 (MQ=255)
               tGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGc                                                                                                                   >  1:77630/1‑139 (MQ=255)
                      aTAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCTACCGCACTTTACTTCTTCTTCGCTGCAGCTCTgcgc                                                                                                            >  2:502232/1‑139 (MQ=255)
                      aTAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTgcgc                                                                                                            >  2:288357/1‑139 (MQ=255)
                         acacTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGt                                                                                                         <  1:288357/139‑1 (MQ=255)
                         acacTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGt                                                                                                         <  2:591156/139‑1 (MQ=255)
                                  gCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCaaa                                                                                                <  1:361009/139‑1 (MQ=255)
                                         gcggcgTGACCGCCTTATCCGGCCAACGATTCGAGATGAATTAGTAGGCAGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCACTGCGCCGTCTCTCCAAATCCTTCa                                                                                         >  1:88086/1‑139 (MQ=255)
                                                       ttATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACg                                                                           >  2:749593/1‑139 (MQ=255)
                                                              gCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCcatcat                                                                    <  2:464496/139‑1 (MQ=255)
                                                                         gggATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTg                                                                                                               >  1:499804/1‑85 (MQ=255)
                                                                         gggATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTg                                                                                                               <  2:499804/85‑1 (MQ=255)
                                                                         gggATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATc                                                         >  1:270533/1‑139 (MQ=255)
                                                                                     tagGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGtt                                             >  1:791137/1‑139 (MQ=255)
                                                                                         ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                         >  2:582652/1‑139 (MQ=255)
                                                                                         ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                         >  1:143771/1‑139 (MQ=255)
                                                                                         ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                         >  2:87530/1‑139 (MQ=255)
                                                                                                       tCAAGCATCGCATCCGGTAACAACCGCACTATACTTCTACTTCGCTGCAGCCCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCAATCAACCGCCAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGc                           >  1:478988/1‑139 (MQ=255)
                                                                                                        cAAGCATCGCGTCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCATCACGGCGCCTTCAGTCCAAATCCTTCAGCAACTTTTTCACGCCATCATCGGCAAACATCGACTCTCGCGTTGCGGAAAGCTGGCGTCGcc                          <  2:88086/139‑1 (MQ=255)
                                                                                                        cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                          <  2:595388/139‑1 (MQ=255)
                                                                                                            catcgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGtt                      >  1:17019/1‑139 (MQ=255)
                                                                                                              tcgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGtttt                    >  2:562872/1‑139 (MQ=255)
                                                                                                                gcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTta                  >  2:599111/1‑139 (MQ=255)
                                                                                                                gcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTta                  <  1:115428/139‑1 (MQ=255)
                                                                                                                 catcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTtat                 >  1:62881/1‑139 (MQ=255)
                                                                                                                         aTCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTaa         <  2:549391/139‑1 (MQ=255)
                                                                                                                          tCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAAc        <  1:708357/139‑1 (MQ=255)
                                                                                                                            aaCCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAG‑GTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTgg     >  1:360222/1‑139 (MQ=255)
                                                                                                                                gCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGc  <  1:628049/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
TGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGC  >  NC_000913/3547853‑3548119

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