Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,463,031 A→T 7.9% intergenic (‑100/+42) 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,0310AT 7.9% 142.6 / 2.1 48intergenic (‑100/+42)algI/algLalginate O‑acetylase AlgI/alginate lyase
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base T (2/2);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCC  >  NC_002947/1462884‑1463177
                                                                                                                                                   |                                                                                                                                                  
gAAGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAg                                                                                                                                                  <  1:91710/150‑1 (MQ=255)
                        gaacaCGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGccccc                                                                                                                          >  2:207123/1‑150 (MQ=255)
                          acacGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGc                                                                                                                        >  2:462708/1‑150 (MQ=255)
                          acacGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGATGGAAGGCAACTTAGATCCCTGGGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGc                                                                                                                        >  1:674803/1‑150 (MQ=255)
                             cGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAAt                                                                                                                     <  1:516309/150‑1 (MQ=255)
                              gTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATg                                                                                                                    >  2:372873/1‑150 (MQ=255)
                                tGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGTATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATggg                                                                                                                  <  2:973833/150‑1 (MQ=255)
                                 ggACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATgggg                                                                                                                 <  1:462708/150‑1 (MQ=255)
                                       gaagaCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCTCCGCAATGGGGGTCAAt                                                                                                           >  2:33747/1‑150 (MQ=255)
                                        aagaCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGccccc                                                                                                                          >  1:1055201/1‑134 (MQ=255)
                                        aagaCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGccccc                                                                                                                          <  2:1055201/134‑1 (MQ=255)
                                           aCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATgggg                                                                                                                 <  2:1120729/140‑1 (MQ=255)
                                           aCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATgggg                                                                                                                 >  1:1120729/1‑140 (MQ=255)
                                           aCCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCAt                                                                                                       >  2:433573/1‑150 (MQ=255)
                                             cATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATgg                                                                                                                   >  1:1113479/1‑136 (MQ=255)
                                             cATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATgg                                                                                                                   <  2:1113479/136‑1 (MQ=255)
                                                       ctcCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGctctct                                                                                           <  2:857109/150‑1 (MQ=255)
                                                       ctcCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGctctct                                                                                           <  2:276311/150‑1 (MQ=255)
                                                           tAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATggggg                                                                                                                >  2:649788/1‑125 (MQ=255)
                                                           tAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATggggg                                                                                                                <  1:649788/125‑1 (MQ=255)
                                                           tAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCAc                                                                                       >  2:1416028/1‑150 (MQ=255)
                                                              tttGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCg                                                                                    <  1:372873/150‑1 (MQ=255)
                                                               ttGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCgg                                                                                   <  1:1111359/150‑1 (MQ=255)
                                                                   tCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTc                                                                                                                                           >  2:913645/1‑90 (MQ=255)
                                                                   tCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCTTGCCCGCTc                                                                                                                                           <  1:913645/90‑1 (MQ=255)
                                                                     acaGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTa                                                                             >  2:891000/1‑150 (MQ=255)
                                                                     acaGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTa                                                                             >  1:1145023/1‑150 (MQ=255)
                                                                            ccAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCcgc                                                                      <  2:1145023/150‑1 (MQ=255)
                                                                             gaGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCcgcg                                                                     >  2:979420/2‑150 (MQ=255)
                                                                              aGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGt                                                                    >  2:829172/1‑150 (MQ=255)
                                                                                            cccAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCa                                                      >  2:753946/1‑150 (MQ=255)
                                                                                              cAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTAATCTCTGCACGAGGGCCGTAGGCCCGCGTCACATCTCCGCCCAGc                                                    >  2:308523/1‑150 (MQ=255)
                                                                                               aacccccccAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGcc                                                   <  2:674803/147‑1 (MQ=255)
                                                                                                        aGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGt                                          <  1:1416028/150‑1 (MQ=255)
                                                                                                            cAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGaa                                      <  2:976845/150‑1 (MQ=255)
                                                                                                             aaGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAACCGGAAGCTGTTGAAc                                     >  2:970396/1‑150 (MQ=255)
                                                                                                                  tGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTa                                                                             >  2:1174003/1‑105 (MQ=255)
                                                                                                                  tGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTa                                                                             <  1:1174003/105‑1 (MQ=255)
                                                                                                                        gCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTc                          <  1:979420/150‑1 (MQ=255)
                                                                                                                             ttAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCt                     >  1:1048272/1‑150 (MQ=255)
                                                                                                                             ttAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCACCCGGAAGCTGTTGAACGGCTGCATGTCGTGCt                     >  1:243020/1‑150 (MQ=255)
                                                                                                                                 aTCCCTGTGGGAGCGGGCCAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTgcg                 <  1:433573/150‑1 (MQ=255)
                                                                                                                                        tgGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGc          >  1:491352/1‑150 (MQ=255)
                                                                                                                                         gggAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCg         >  1:225697/1‑150 (MQ=255)
                                                                                                                                           gAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTc                                                                                <  1:1169213/77‑1 (MQ=255)
                                                                                                                                           gAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTc                                                                                >  2:1169213/1‑77 (MQ=38)
                                                                                                                                               gggCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAACCGGAAGCTGTTGAACGGCTGCAGGTCGTGCTTGCGCTCAAGCGTGTCGc   <  1:970396/150‑1 (MQ=255)
                                                                                                                                                ggCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCGTTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGcc  >  2:790704/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAAGATCGGCAAGAACAGGAACAGGAACACGTTGGACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCCCAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTCGCC  >  NC_002947/1462884‑1463177

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