Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,649 C→T 28.6% intergenic (+256/+459) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,6490CT28.6% 27.3 / 10.9 21intergenic (+256/+459)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base T (3/3);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

ACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAC  >  NC_000913/2304505‑2304796
                                                                                                                                                |                                                                                                                                                   
acacCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAAGGCCCGAGCCGGTAGCCTGATGCGAAGCTGCCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGATAAACCCCGTCGGTCGGAt                                                                                                                                                 >  1:210780/1‑149 (MQ=14)
             gTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGCGCCGGTTGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAACCCGTAGGTCGGATAAGGCGGTTAcgc                                                                                                                                    >  7:224824/1‑149 (MQ=255)
                              cccgACAAGTCGGTCACGCCGCATCCGACCTCCAACGCCCGAGCCGCATGCCACAGGCGACCCGCGCGCGGCTCATCCGCCCTACACCGCCGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                                                                                                   <  3:182148/146‑1 (MQ=11)
                                     aGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                                                                                                                                                     <  3:90257/108‑1 (MQ=35)
                                     aGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                                                                                                                                                     >  4:90257/1‑108 (MQ=35)
                                      ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                               <  1:134997/113‑1 (MQ=35)
                                      ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                               >  2:134997/1‑113 (MQ=35)
                                      ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTCCGCCGCCTCCGCCATCTAACGCCCAAGc                                                                                                           >  3:208551/1‑149 (MQ=255)
                                      ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACTCCGCATCCGACATCTAACGCCCAAGc                                                                                                           >  1:38964/1‑149 (MQ=21)
                                      ggCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGc                                                                                                           <  7:17943/149‑1 (MQ=34)
                                                    aTCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGGTGGCGCGTCTTAGCAGGCGTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                               <  1:42430/147‑1 (MQ=255)
                                                     tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                                                                                                                                               >  1:157692/1‑98 (MQ=18)
                                                     tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                                                                                                                                               <  6:236766/98‑1 (MQ=18)
                                                     tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                                                                                                                                               <  6:233231/98‑1 (MQ=18)
                                                     tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                                                                                                                                               >  5:236766/1‑98 (MQ=18)
                                                     tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                                                                                                                                               >  5:233231/1‑98 (MQ=18)
                                                     tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                                                                                                                                               <  2:157692/98‑1 (MQ=18)
                                                           aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCg                                                                                                                                           >  2:206020/1‑96 (MQ=33)
                                                           aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCg                                                                                                                                           <  1:206020/96‑1 (MQ=33)
                                                           aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAc                                                                                                                          <  7:47305/113‑1 (MQ=32)
                                                           aTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAc                                                                                                                          >  8:47305/1‑113 (MQ=32)
                                                                                                                                       cccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGa          >  4:218929/1‑149 (MQ=31)
                                                                                                                                               tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAc  >  7:79689/1‑149 (MQ=34)
                                                                                                                                                |                                                                                                                                                   
ACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAAC  >  NC_000913/2304505‑2304796

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