Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,470,735 G→A 10.0% H188H (CAC→CAT)  alg44 ← alginate biosynthesis protein Alg44

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,470,7350GA10.0% 106.8 / 2.5 40H188H (CAC→CAT) alg44alginate biosynthesis protein Alg44
Reads supporting (aligned to +/- strand):  ref base G (15/21);  new base A (2/2);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

CCTTGAGCATGTCCAGCATGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCTTGAAGCGGCCGAAGGCGCTCAGGCCCGACCCACCGTCCTTCTGC  >  NC_002947/1470587‑1470873
                                                                                                                                                    |                                                                                                                                          
ccTTGAGCATGTCCAGCATGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCg                                                                                                                                            >  4:55967/1‑149 (MQ=255)
 cTTGAGCATGTCCAGCATGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGt                                                                                                                                           <  6:64319/149‑1 (MQ=255)
 cTTGAGCATGTCCAGCATGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGt                                                                                                                                           <  7:335894/149‑1 (MQ=255)
     aGCATGTCCAGCATGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGc                                                                                                                                        <  3:55967/148‑1 (MQ=255)
       cATGTCCATCATGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGa                                                                                                                                     <  2:243710/149‑1 (MQ=255)
              aGCATGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATa                                                                                                                              <  3:335853/149‑1 (MQ=255)
                  tGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATg                                                                                                                          <  1:178574/149‑1 (MQ=255)
                     tggtggTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCa                                                                                                                       >  4:129227/1‑149 (MQ=255)
                                   gCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGCCCACACCGGCTTCGGCGTGGCTGACGAAATCCATGCCATACAACGATTTGGc                                                                                                         >  8:230243/1‑149 (MQ=255)
                                     cAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGa                                                                                                                                     <  3:105651/119‑1 (MQ=255)
                                     cAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGa                                                                                                                                     >  4:105651/1‑119 (MQ=255)
                                      aaTGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGcgac                                                                                                      <  4:95697/149‑1 (MQ=255)
                                           cgcgcCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGa                                                                                                                                     <  7:339681/113‑1 (MQ=255)
                                                            gCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCATGGc                                                                                                                                        <  8:217611/93‑1 (MQ=255)
                                                            gCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCATGGc                                                                                                                                        >  7:217611/1‑93 (MQ=255)
                                                            gCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCATGGc                                                                                                                                        <  6:198341/93‑1 (MQ=255)
                                                            gCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCATGGc                                                                                                                                        >  5:198341/1‑93 (MQ=255)
                                                               cACCGCTTTCGACCAGGCTGCTTACGGTCCCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCg                                                                                        <  8:258319/138‑1 (MQ=255)
                                                               cACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCg                                                                                        >  7:258319/1‑138 (MQ=255)
                                                                aCCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCa                                                                            >  1:68105/1‑149 (MQ=255)
                                                                     tttCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGaa                                                                        >  5:286982/1‑148 (MQ=255)
                                                                      ttCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAaca                                                                      >  7:23696/1‑149 (MQ=255)
                                                                              ggCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCaga                                                              <  3:129227/149‑1 (MQ=255)
                                                                              ggCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCaga                                                              <  2:68105/149‑1 (MQ=255)
                                                                                         gTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGGGCCGGTgcc                                                   >  4:89409/1‑147 (MQ=255)
                                                                                          tACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACt                                                  >  6:80942/1‑149 (MQ=255)
                                                                                           aCCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCAt                                                                                                                      >  5:69601/1‑80 (MQ=255)
                                                                                           aCCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCAt                                                                                                                      <  6:69601/80‑1 (MQ=255)
                                                                                              gTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCt                                              >  2:58125/1‑149 (MQ=255)
                                                                                                      gCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCTTGAAGCgg                                      <  5:80942/149‑1 (MQ=255)
                                                                                                            tGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCTTGAAGCGGCCGAAg                                <  3:89409/149‑1 (MQ=255)
                                                                                                                gTTGGTGGTCGGCAATGCGACCAAACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCTTGAAGCGGCCGAAGgcgc                            >  4:66961/1‑149 (MQ=255)
                                                                                                                 ttGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACt                                                  <  1:334112/126‑1 (MQ=255)
                                                                                                                 ttGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACt                                                  >  2:334112/1‑126 (MQ=255)
                                                                                                                     tggtCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGt                                                      >  7:320214/1‑118 (MQ=255)
                                                                                                                     tggtCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGt                                                      <  8:320214/118‑1 (MQ=255)
                                                                                                                                 cGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACTAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGAc                                                         <  3:50001/103‑1 (MQ=255)
                                                                                                                                 cGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCTTGAAGCGGCCGAAGGCGCTCAGGCCCGACCCAcc            <  5:87409/148‑1 (MQ=255)
                                                                                                                                 cGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACAAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGAc                                                         >  4:50001/1‑103 (MQ=255)
                                                                                                                                        acCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCTTGAAGCGGCCGAAGGCGCTCAGGCCCGACCCACCGTCCTTCt    <  3:66961/149‑1 (MQ=255)
                                                                                                                                          cGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCTTGAAGCGGCCGAAGGCGCTCAGGCCCGACCCACCGTCCTTCTGc  <  1:58125/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
CCTTGAGCATGTCCAGCATGCTGGTGGTGAAGCTGGCCAATGGCGCGCCCTTGGCGATCTGCCCACCGCTTTCGACCAGGCTGCTTACGGTACCGTCGCGCGGCATGGTGACGTTGGTGGTCGGCACTGCGACCACACCGGCTTCGGCGTGGCTGACGAAATACATGCCATACAACGATTTGGCGACGAAGCCGAAGGCCGCAACGCCGACCACGAACACGCCCAGAGTGACGGTGACTGCCTTGAAGCGGCCGAAGGCGCTCAGGCCCGACCCACCGTCCTTCTGC  >  NC_002947/1470587‑1470873

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