Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,463,032 A→T 17.4% intergenic (‑101/+41) algI ← / ← algL alginate O‑acetylase AlgI/alginate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,463,0320AT17.4% 51.3 / 5.8 23intergenic (‑101/+41)algI/algLalginate O‑acetylase AlgI/alginate lyase
Reads supporting (aligned to +/- strand):  ref base A (11/8);  new base T (2/2);  total (13/10)
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.29e-01

AGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCC  >  NC_002947/1462886‑1463177
                                                                                                                                                  |                                                                                                                                                 
aGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGcc                                                                                                                                                <  1:336879/150‑1 (MQ=255)
          gaacaggaacagGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTcc                                                                                                                                          <  2:338571/146‑1 (MQ=255)
          gaacaggaacagGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTcc                                                                                                                                          >  1:338571/1‑146 (MQ=255)
             caggaacaggaaCACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGt                                                                                                                                   <  2:21658/150‑1 (MQ=255)
                                        gaCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCa                                                                                                        >  1:120883/1‑150 (MQ=255)
                                         aCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCAt                                                                                                       <  2:120883/150‑1 (MQ=255)
                                               ccGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCAAGTGCACGACTGGAAGGACAAGGAGAGCCCTGGGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTg                                                                                                 >  1:60057/1‑150 (MQ=255)
                                                     ctcCTAATTTGTTCACAGCATCCAGGACAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGctctct                                                                                           <  2:58857/150‑1 (MQ=255)
                                                                     aGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGTGGTCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGa                                                                           <  1:276174/150‑1 (MQ=255)
                                                                          ccAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCcgc                                                                      >  2:177217/1‑150 (MQ=255)
                                                                              ggCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCGACAGGGTCCGCCCCCGCAAGGGGGGTCAAGTCATTTCTTGCTCTCTGCACACGGTGCGTAGACCCGCGTca                                                                  >  1:271481/1‑150 (MQ=255)
                                                                               ccAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGTTCGTAGACCCGCGTcac                                                                 <  1:100180/149‑1 (MQ=255)
                                                                                           ccAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGGGGCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAg                                                     >  1:34563/1‑150 (MQ=255)
                                                                                                     cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTg                                           >  1:376337/1‑150 (MQ=255)
                                                                                                                    aGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCa                                                                                                        <  1:398159/74‑1 (MQ=255)
                                                                                                                    aGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCa                                                                                                        >  2:398159/1‑74 (MQ=255)
                                                                                                                                   cTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCCGTTGAACGGCTGCATGTCGTGCTTGCGCTCa             >  1:347764/1‑150 (MQ=255)
                                                                                                                                    ttttggAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGACCAGCTGCATGTCGTGCTTGCGCTCaa            >  1:17408/5‑150 (MQ=255)
                                                                                                                                    tgtgGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCaa            >  1:27472/1‑150 (MQ=255)
                                                                                                                                        ggAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCgt        <  2:186832/150‑1 (MQ=255)
                                                                                                                                        ggAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCgt        >  1:186832/1‑150 (MQ=255)
                                                                                                                                         gAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGAATGTCGTGCTTGCGCTAAAGCgtg       >  1:462006/1‑150 (MQ=255)
                                                                                                                                              ggCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGcc  <  2:27472/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCC  >  NC_002947/1462886‑1463177

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