Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,463,025 A→T 17.3% 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,0250AT17.3% 62.5 / 9.3 35intergenic (‑101/+41)algI/algLalginate O‑acetylase AlgI/alginate lyase
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base T (4/2);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTC  >  NC_002947/1462911‑1463167
                                                                                                                   |                                                                                                                                              
gACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAGGCGCCCCCCAACCCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCAGTGGGGGGGGGCACGCCCGCGCCCACCGCTTCCGCCCCGGCAATgggg                                                                                                              >  8:71733/1‑149 (MQ=255)
 aCGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATgggg                                                                                                              <  7:225330/148‑1 (MQ=255)
  cGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGt                                                                                                            <  5:225305/149‑1 (MQ=255)
                           aTTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCCACTTAGATCCCTGTGGGGGCGGGCAAGCCCGCTCCCACACGTTCCGCCCCCGCAATGGGGGGCAATTCAATTCTTGCTCTATGccca                                                                                   >  8:18184/1‑146 (MQ=255)
                                ttCACAGCATCCAGGGCAAACCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCAACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCggg                                                                               <  2:244409/148‑1 (MQ=255)
                                  cacaGCATCCAGGGCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTc                                                                             <  5:78655/148‑1 (MQ=255)
                                    caGCATCCAGGGCAAAGCGCCCCCCAAACCCCCATAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTccc                                                                                                                                      <  1:74161/90‑1 (MQ=39)
                                            aGGTCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTc                                                                                  <  3:109003/133‑1 (MQ=255)
                                            aGGGCAAAGCGCCCCCCCAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCGCCCACAGGTTGCGCCCCCGCAATGGGGGTCAATTCATTTCTTGGGCTCTGCCCACGGGTGGTAGACCcgcg                                                                  >  6:108582/1‑149 (MQ=255)
                                            aGGGCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTa                                                                          <  5:228400/141‑1 (MQ=255)
                                            aGGGCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTAATTTCTTGCTCTCTGCACTCggggagta                                                                          >  6:228400/1‑136 (MQ=255)
                                               gCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTca                                                               <  7:18184/149‑1 (MQ=255)
                                               gCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACGTAGATACCTGTGGGGGCGGGCAAGCCCGCGCCCACAGGTTCCGCCCCCGCAATGGGGGTGCAGTGATTTGTTGGGGTCTGCACTCGGGTCGTAGACCCGCGTca                                                               >  4:123415/1‑149 (MQ=255)
                                                  aaGCGCCCCCCCAACCCCC‑CAGGGCACCGCGGGACGGCAACCAAGCTCCCTGTGCGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAAATCATTTCTTGCTCTCTGCACTCGGGTCTACGACCCGAGtgcaga                                                            <  8:89669/149‑6 (MQ=255)
                                                          cccAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCgccgcc                                                     >  2:255028/1‑148 (MQ=255)
                                                          cccAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCgccgcc                                                     >  1:261524/1‑148 (MQ=255)
                                                           ccAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCa                                                   <  1:248296/149‑1 (MQ=255)
                                                                  ccc‑cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCgg                                              >  3:15029/1‑147 (MQ=255)
                                                                      cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATcgcc                                                        <  1:169788/134‑1 (MQ=255)
                                                                      cAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATcgcc                                                        >  2:169788/1‑134 (MQ=255)
                                                                            aaGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACACGGGTCGTAGACCCGCGTCACAAAGCCGCCCACCCGGAAGCGGTTGa                                    >  2:126853/1‑148 (MQ=255)
                                                                               gCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACgg                                <  5:108582/149‑1 (MQ=255)
                                                                                          aCTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGt                     <  2:261524/149‑1 (MQ=255)
                                                                                            ttAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGc                   >  1:140729/1‑149 (MQ=255)
                                                                                                 tCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTgc               <  1:255028/148‑1 (MQ=255)
                                                                                                    cTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACAt                                                            >  6:170503/1‑100 (MQ=255)
                                                                                                     tgtgGGAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGt                     <  1:58410/138‑1 (MQ=255)
                                                                                                     tgtgGGAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGt                     >  2:58410/1‑138 (MQ=255)
                                                                                                      gtgGGAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAAc                                  >  7:6287/1‑124 (MQ=255)
                                                                                                      gtgGGAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAAc                                  <  8:6287/124‑1 (MQ=255)
                                                                                                         ggAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCg      >  8:60696/1‑149 (MQ=255)
                                                                                                         ggAGCGGGCTTGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCg      >  4:236414/1‑149 (MQ=255)
                                                                                                         ggAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGc       <  2:140729/148‑1 (MQ=255)
                                                                                                           aGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAAGTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCgtg    >  1:199676/1‑149 (MQ=255)
                                                                                                             cGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTc  >  6:257091/1‑149 (MQ=255)
                                                                                                                   |                                                                                                                                              
GACGAGAAGACCATCCCGGTTCTCCTAATTTGTTCACAGCATCCAGGGCAAAGCGCCCCCCAAACCCCC‑CAGTGCAAGGCTGGAAGGCAACTTAGATCCCTGTGGGAGCGGGCAAGCCCGCTCCCACAGGTTCCGCCCCCGCAATGGGGGTCAATTCATTTCTTGCTCTCTGCACTCGGGTCGTAGACCCGCGTCACATCGCCGCCCAGCCGGAAGCTGTTGAACGGCTGCATGTCGTGCTTGCGCTCAAGCGTGTC  >  NC_002947/1462911‑1463167

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: ‑