Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 141,137 T→A 10.2% Q237L (CAG→CTG)  algB ← alginate biosynthesis transcriptional regulator AlgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947141,1370TA10.2% 114.1 / 4.6 39Q237L (CAG→CTG) algBalginate biosynthesis transcriptional regulator AlgB
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base A (2/2);  total (20/19)
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.49e-01

TGGCAGCCAGGATACGTACATCGGCACGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCAGGCTTTGCGTGCACGCTTGCTCCAACCATGGAT  >  NC_002947/140988‑141285
                                                                                                                                                     |                                                                                                                                                    
ttgCAGCCAGGATACGTACATCGGCACGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCt                                                                                                                                                      <  1:441417/148‑1 (MQ=255)
       cAGGATACGTACATCGGCACGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCAGGCTGAc                                                                                                                                               >  1:592930/1‑150 (MQ=255)
          gATACGTACATCGGCACGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCAGGCTGAcggg                                                                                                                                            <  2:592930/150‑4 (MQ=255)
           aTACGTACATCGGCACGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGAcccgc                                                                                                                                           <  2:158270/150‑1 (MQ=255)
               gTACATCGGCACGCCGGGTAATGGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCAGGCTGAcggg                                                                                                                                            <  1:429899/145‑4 (MQ=255)
               gTACATCGGCACGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACATTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGa                                                                                                                                       >  2:348780/1‑150 (MQ=255)
               gTACATCGGCACGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCAGGCTGAcggg                                                                                                                                            >  2:429899/1‑142 (MQ=255)
                         aCGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATATAGCGGAAGAATGTCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTc                                                                                                                             >  2:663237/1‑150 (MQ=255)
                                     cGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTTCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGt                                                                                                                 >  2:355643/1‑150 (MQ=255)
                                     cGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGt                                                                                                                 <  1:710895/150‑1 (MQ=255)
                                       ggTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGa                                                                                                               >  1:305525/1‑150 (MQ=255)
                                             cGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGgcgc                                                                                                         <  1:513919/150‑1 (MQ=255)
                                             cGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGCACTCTGGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGgcgc                                                                                                         <  2:292725/150‑1 (MQ=255)
                                              gACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGcgcc                                                                                                        >  2:530010/1‑150 (MQ=255)
                                                   gTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCggg                                                                                                   >  2:678577/1‑150 (MQ=255)
                                                    ttCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCgggg                                                                                                  >  2:836911/1‑149 (MQ=255)
                                                          ttCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGAcc                                                                                            <  1:111216/150‑1 (MQ=255)
                                                               tGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATa                                                                                       <  2:485923/150‑1 (MQ=255)
                                                                gTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAg                                                                                      <  1:678577/150‑1 (MQ=255)
                                                                    tGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCTTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCg                                                                                  <  1:530010/150‑1 (MQ=255)
                                                                           aaCGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCa                                                                           >  2:308108/1‑150 (MQ=255)
                                                                            aCGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGa                                                                                                                                                <  1:680903/80‑1 (MQ=255)
                                                                            aCGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGa                                                                                                                                                >  2:680903/1‑80 (MQ=255)
                                                                              gCAACATTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCa                                                                        <  1:348780/150‑1 (MQ=255)
                                                                                aaCAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAAt                                                                                        <  2:304823/132‑1 (MQ=255)
                                                                                aaCAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAAt                                                                                        >  1:304823/1‑132 (MQ=255)
                                                                                    gTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCgg                                                                  >  1:441976/1‑150 (MQ=255)
                                                                                            tGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTACATCAGCTCGGCGTTAAGc                                                          >  1:44124/1‑150 (MQ=255)
                                                                                                       cGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGtt                                               >  1:250150/1‑150 (MQ=255)
                                                                                                                   gATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCa                                   >  1:399524/1‑150 (MQ=255)
                                                                                                                          tcgAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCAGGCTTTg                            >  2:560468/1‑150 (MQ=255)
                                                                                                                           cgAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCAGGCTTTGc                           <  2:305525/150‑1 (MQ=255)
                                                                                                                                      gTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCa                                       >  2:145190/1‑127 (MQ=255)
                                                                                                                                      gTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCa                                       <  1:145190/127‑1 (MQ=255)
                                                                                                                                       tGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCAc                                      >  2:514943/1‑127 (MQ=255)
                                                                                                                                       tGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCAc                                      <  1:514943/127‑1 (MQ=255)
                                                                                                                                          ccccGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCAGGCTTTGCGTGCACGCTTGCTcc            >  2:810519/1‑150 (MQ=255)
                                                                                                                                          ccccGTCAGCCTGGCTGACCCGACCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCAGGCTTTGCGTGCACGCTTGCTcc            <  2:218373/150‑1 (MQ=255)
                                                                                                                                                   ccTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCAGGCTTTGCGTGCACGCTTGCTCCAACCatgga   <  2:441976/150‑1 (MQ=255)
                                                                                                                                                    cTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCAGGCTTTGCGTGCACGCTTGCTCCAACCatggat  <  1:308108/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCAGCCAGGATACGTACATCGGCACGCCGGGTAACCGGGTCGCCGACGCGTTCGTATTCCTTGTCCTGGATGAAACGCAACAGTTTGGGCTGCAAGGTCAGCGGGAAATCGCCGATCTCGTCGAGGAATAGCGTGCCCCCGTCAGCCTGGCTGACCCGCCCGAGGGTGCTTTCGCTGGCGCCGGTGAAGGCGCCGCGGGTATGACCGAATAGCTCGCTTTCCATCAGCTCGGCGTTAAGCGACGGGCAGTTGATGGTCACGCAGGCTTTGCGTGCACGCTTGCTCCAACCATGGAT  >  NC_002947/140988‑141285

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