Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,649 C→T 11.8% 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,6490CT11.8% 71.0 / 4.2 34intergenic (+256/+459)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base T (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

TACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGAC  >  NC_000913/2304504‑2304789
                                                                                                                                                 |                                                                                                                                            
tACACCGCTGTGAAGAGCGGCACCCCGTAGGTAAGACAAAGCGGTCAAGCCGCATCCGACATCAAACGCCCCAGCCGGTTCCCTCATGATACGATGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTGAACCCCGGAGTTCGGa                                                                                                                                           >  6:356693/1‑149 (MQ=14)
   acCGCTGGGAAGTGCGGCACCCCGTAGGTCTGCCAAGGCGGCCACTCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                        <  7:356111/149‑1 (MQ=14)
     cgccgggAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCTAGCCTGTTGCCTGAAGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                                                                                                                                      <  2:247303/143‑1 (MQ=11)
            aaGTGCGGCCCCCCGTAGGTCAGACACTGCGGTCACGCCGCATCCGACATCCCACGCCCGAGCCGCTTGCCTGAGGCGACGCTGGCGCGTCTTATCATGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTa                                                                                                                                <  1:258115/148‑1 (MQ=11)
             agcgcgcCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcg                                                                                                                              <  7:200106/146‑1 (MQ=255)
                gCGGCCCCCCGTAGTTCATACAAGCAGCTCACGCCTCAACCGACATCCTACGCCCGAGGCGGTTGCCACATGCGACGATGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccg                                                                                                                           <  7:15594/149‑1 (MQ=14)
                       cccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGa                                                                                                                    >  7:257997/1‑149 (MQ=255)
                          gTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACAt                                                                                                                 <  8:68881/149‑1 (MQ=255)
                                 aGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTACCGc                                                                                                          >  7:219285/1‑149 (MQ=38)
                                   aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                                                                                                           <  6:272241/114‑1 (MQ=38)
                                   aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                                                                                                           >  5:272241/1‑114 (MQ=38)
                                   aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCCTCATATCAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCCGTGACGCCGCATCCGACATCTCACGccc                                                                                                        >  7:175387/1‑149 (MQ=11)
                                     aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                       >  8:374123/1‑148 (MQ=34)
                                            tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                               <  8:257997/149‑1 (MQ=18)
                                               cgacgcAGCCCACATCCATAGCTCGAGCCGCGTGCCGGATGCTCGGCTGGCGATTCTAATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGcc                                                                                            <  3:201311/146‑1 (MQ=11)
                                                gccgcATCCGACATCCAACGCCCTAGCCTGTTGCCTGATGCGACGCTGTCGCGTATTATCAGGCCTACACATCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                           <  8:219285/149‑1 (MQ=11)
                                                      tCCGACATCCAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                                                                                                                                              >  7:342784/1‑92 (MQ=18)
                                                      tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                                                                                                         >  1:72557/1‑97 (MQ=33)
                                                      tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                                                                                                                                         <  2:72557/97‑1 (MQ=33)
                                                       ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGc                                                                                                          <  4:157421/127‑1 (MQ=32)
                                                       ccGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGc                                                                                                          >  3:157421/1‑127 (MQ=32)
                                                         gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCg                                                                                                                                    >  4:306298/1‑99 (MQ=18)
                                                         gACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCg                                                                                                                                    <  3:306298/99‑1 (MQ=18)
                                                             tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTaa                                                                                                             >  7:145435/1‑118 (MQ=32)
                                                             tCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTaa                                                                                                             <  8:145435/118‑1 (MQ=32)
                                                              ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcgccgc                                                                                                                          >  3:111653/1‑104 (MQ=14)
                                                              ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTAcgccgc                                                                                                                          <  4:111653/104‑1 (MQ=14)
                                                              ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                             >  8:47870/1‑149 (MQ=32)
                                                              ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                             >  5:278355/1‑149 (MQ=32)
                                                               cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                        <  7:361626/89‑1 (MQ=33)
                                                               cAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                        >  8:361626/1‑89 (MQ=33)
                                                                   gCCCGAGGCGTTGGCCTGATGCTAAGCTGGCGAGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCTGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCTAGCCGGTTGCCTGATGCGACGCTGGCACGTc                                                                        <  6:338549/149‑1 (MQ=11)
                                                                   gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTACCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                        >  1:279052/1‑149 (MQ=14)
                                                                                                                                        cccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGGAGGTTGGATAAGGCGGTTACGCCGCATCCGa   >  5:65042/1‑149 (MQ=14)
                                                                                                                                          cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGGTACGCCGCAGCCGAc  >  1:98825/1‑148 (MQ=17)
                                                                                                                                                 |                                                                                                                                            
TACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGAC  >  NC_000913/2304504‑2304789

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