Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,304,575 G→A intergenic (+182/+533) eco → / ← mqo serine protease inhibitor ecotin/malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,5750GA91.7% 14.8 / ‑3.3 12intergenic (+182/+533)eco/mqoserine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (3/8);  total (3/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 = 9.41e-01

GCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAC  >  NC_000913/2304285‑2304847
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                 
gcGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGtt                                                                                                                                                                                                                                                                           >  1:31633/1‑299 (MQ=255)
                     gATCGTGGTGTATACGCCAGACAAGGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                                                                     <  2:712976/300‑1 (MQ=17)
                             tgtATACGCCAGACAATGTAGCTGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCTTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGTCCTACACCGCTTTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAg                                                                                                                                                                                                                                            <  2:407899/300‑1 (MQ=17)
                              gtATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCa                                                                                                                                                                                                                                             >  1:108577/1‑298 (MQ=255)
                                   cGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCt                                                                                                                                                                                                                                        >  1:270754/1‑298 (MQ=255)
                                       aGACGATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACACCGCGGTAGTTCGTTAAACTGCCGTGAAGTGTGGCAGCTAGTGGGTTAGATAAGGCGGTCATGTGGCATCTGATATCCAACGCCCTAGCTGGTTGCCTGATGCTACGCTGGCGTGTCTTATCGGGCCTACACCGCTGTGAAGTGCGTCACCCCGTAGGTCGGGTAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGTGT‑TTTATCAggtttacatt                                                                                                                                                                                                                                   <  2:237580/299‑9 (MQ=12)
                                          aatggTAGATGTGAAGNTACCGCGTCTGGAAGGCGNGAAGAGAGAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTNCTTATCAGGCCTacac                                                                                                                                                                                                                                    <  1:700276/295‑1 (MQ=14)
                                           acgctAGATGTGAAGNTACCGCGTCTGGAAGGCGNGAAGAGAGAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTNCTTATCAGGCCTacac                                                                                                                                                                                                                                    <  1:700277/294‑1 (MQ=14)
                                                   atgttaaggaccgcggtgggAAGGGG‑GAAGGTA‑AAATTGACAGCTCGTTAGTTCGCTAAACTGCCGTTAAATGCGGCACCCCTTAGGTCAGACACGTCGGTCACGCCGCATCCGACATCCAACGCCCGAGCTGGTTGCCGGATTCGACACTGTCGCTTCTGATCATGCCTACACCGCTTTTAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCAAGCCGGTTGCCTGAGGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGc                                                                                                                                                                                                                       <  2:463229/282‑1 (MQ=2)
                                                                 gTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGGAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGt                                                                                                                                                                                                         <  1:415940/300‑1 (MQ=11)
                                                                     ggAAGGCG‑GAAGAGA‑AAATCGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGTTCAGACAAGGCTGTCACGCCGCATCCGACAGCCAACGCCCGAACCGGTTGCCTGATGTGACACTGGCGCGTCGGATCATGCCTACACTGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAg                                                                                                                                                                                                                                            <  2:409433/261‑1 (MQ=11)
                                                                                                                                                                                                                                                                          tCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac  <  1:204731/300‑1 (MQ=39)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                 
GCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAG‑TACCGCGTCTGGAAGGCG‑GAAGAGA‑AAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGT‑CTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAC  >  NC_000913/2304285‑2304847

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