Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,467,305 C→T 7.9% G465D (GGT→GAT)  algE ← alginate production protein AlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,467,3050CT 7.9% 139.0 / 2.0 50G465D (GGT→GAT) algEalginate production protein AlgE
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base T (2/2);  total (23/27)
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.69e-01

ACTGGCCAACAGGCTGTGACGTAAGTGCGGGTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCAAGGTCCTTCTCGCCGTCGACCAGCTGACGG  >  NC_002947/1467156‑1467445
                                                                                                                                                     |                                                                                                                                            
aCTGGCCAACCGGCCGTGACGTAAGTGCGGGTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCCGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGACTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATacac                                                                                                                                              <  1:1151811/150‑1 (MQ=255)
aCTGGCCAACAGGCTGTGACGTAAGTGCGGGTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATacac                                                                                                                                              <  1:471614/150‑1 (MQ=255)
    gCCAACAGGCTGTGACGTAAGTGCGGGTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGc                                                                                                                                          <  1:124300/150‑1 (MQ=255)
          aGGCTGTGACGTAAGTGGGGGGGACGGTTCATCACGGTCAGCCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAg                                                                                                                                    <  1:496482/150‑1 (MQ=255)
           ggCTGTGACGTAAGTGCGGGTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAgg                                                                                                                                   >  2:541272/1‑150 (MQ=255)
                    gTAAGTGCGGGTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAgg                                                                                                                          >  2:643616/1‑150 (MQ=255)
                     tAAGTGCGGGTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAggg                                                                                                                         <  1:715737/150‑1 (MQ=255)
                     cccaacggggTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAggg                                                                                                                         <  2:610256/143‑1 (MQ=255)
                              gTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCt                                                                                                                >  1:2287/1‑150 (MQ=255)
                               tGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCTCACCAGGGCTGACGGCtc                                                                                                               <  2:710561/150‑1 (MQ=255)
                               tGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCtc                                                                                                               <  2:656896/150‑1 (MQ=255)
                                  aGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGCCGGCtcgtc                                                                                                            >  1:526075/1‑150 (MQ=255)
                                  aGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCtcgtc                                                                                                            <  1:705257/150‑1 (MQ=255)
                                  aGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCtcgtc                                                                                                            <  1:643616/150‑1 (MQ=255)
                                    gTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGa                                                                                                          >  2:570126/1‑150 (MQ=255)
                                    gTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGa                                                                                                          >  1:351336/1‑150 (MQ=255)
                                         tcaCGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGcc                                                                                                     >  1:360164/1‑150 (MQ=255)
                                               tCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTc                                                                                               >  1:653838/1‑150 (MQ=255)
                                                  gTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGACGCTTGCCTCGCTCagc                                                                                            <  2:351336/150‑3 (MQ=255)
                                                   tCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAgg                                                                                                                          <  1:816911/119‑1 (MQ=255)
                                                   tCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAgg                                                                                                                          >  2:816911/1‑119 (MQ=255)
                                                     cccTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGaa                                                                                         >  2:638481/1‑150 (MQ=255)
                                                      ccTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAg                                                                                        >  1:111377/1‑150 (MQ=255)
                                                        tACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACATCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGcc                                                                                      >  2:833747/1‑150 (MQ=255)
                                                        tACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACATCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGcc                                                                                      >  2:833790/1‑150 (MQ=255)
                                                          ccGGCCTCAAAATCCCCAGATAACGTCGACAGACGGGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCTTCGCCCGGCTTTAATACACCCGCACGCAGGCACACCATGACTGACGTCTCGTCGATTGCCTAGCTCATCGAAGCCgg                                                                                    <  2:397025/150‑1 (MQ=255)
                                                             gCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGgcag                                                                                 >  2:1098881/1‑150 (MQ=255)
                                                                   gAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGcc                                                                           <  2:360164/150‑1 (MQ=255)
                                                                       cgcCAGATGACGTCCACGAAGGCGCGTGGCATGTACTAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAAGACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCcttgc                                                                       <  2:546749/150‑1 (MQ=255)
                                                                         ccAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGa                                                                                                          <  2:227788/113‑1 (MQ=255)
                                                                         ccAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGa                                                                                                          >  1:227788/1‑113 (MQ=255)
                                                                            gATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGCCCCGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGaa                                                                  <  1:392937/150‑1 (MQ=255)
                                                                                aCGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTAc                                                              <  2:653838/150‑1 (MQ=255)
                                                                                   tCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTg                                                           <  2:1072567/150‑1 (MQ=255)
                                                                                     gacgaAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGt                                                         >  1:735915/1‑150 (MQ=255)
                                                                                     gacgaAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGt                                                         <  1:638481/150‑1 (MQ=255)
                                                                                       cgaAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCtcg                                                                                                              >  1:553181/1‑95 (MQ=255)
                                                                                       cgaAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCtcg                                                                                                              <  2:553181/95‑1 (MQ=255)
                                                                                          aGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATacac                                                                                                                                              >  2:466206/1‑60 (MQ=255)
                                                                                          aGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATacac                                                                                                                                              <  1:466206/60‑1 (MQ=255)
                                                                                              gcgGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCAc                                                <  2:949505/150‑1 (MQ=255)
                                                                                                           gAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGTCTTGCTTGAAGTACTTGGGCACGACCACGTCCATTTCCTgg                                   <  2:526075/150‑1 (MQ=255)
                                                                                                           gAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACg                                                     >  2:492724/1‑132 (MQ=255)
                                                                                                           gAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTgg                                   >  2:1052637/1‑150 (MQ=255)
                                                                                                           gAGTCCGTTTCCTTGCCGGAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACg                                                     <  1:492724/132‑1 (MQ=255)
                                                                                                              tCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCa                                >  2:566682/1‑150 (MQ=255)
                                                                                                                cGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCAAg                              <  1:541272/150‑1 (MQ=255)
                                                                                                                 gTTTCCTTGCCGTATGCGGCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCAAgg                             >  1:947806/1‑150 (MQ=255)
                                                                                                                      cTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGaa                                                                                         >  2:31352/1‑85 (MQ=255)
                                                                                                                      cTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGaa                                                                                         <  1:31352/85‑1 (MQ=255)
                                                                                                                                 cGTCGCCCGGCTTGAATACATCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCAAGGTCCTTCTCGCCGTCGa             <  1:833790/150‑1 (MQ=255)
                                                                                                                                 cGTCGCCCGGCTTGAATACATCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCAAGGTCCTTCTCGCCGTCGa             <  1:833747/150‑1 (MQ=255)
                                                                                                                                    cGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCAAGGTCCTTCTCGCCGTCGACCa          <  2:735915/150‑1 (MQ=255)
                                                                                                                                            ttGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCAAGGTCCTTCTCGCCGTCGACCAGCTGACgg  >  1:961357/1‑150 (MQ=255)
                                                                                                                                                     ccGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTg                                                                                                   <  1:913316/44‑1 (MQ=255)
                                                                                                                                                     ccGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTg                                                                                                   >  2:913316/1‑44 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
ACTGGCCAACAGGCTGTGACGTAAGTGCGGGTGAAGGTTCATCACGGTCAGTCCCCTACCGGCCTCAGAAGCGCCAGATGACGTCGACGAAGGCGCGGTGCATGTACGAGTCCGTTTCCTTGCCGTAGGCGTCGCCCGGCTTGAATACACCGGCACGCAGGCGCACCAGGGCTGACGGCTCGTCGATTGCCTGGCTCATCGAAGCCGGCAGCAGGCCTTGCTTGAAGTACTTGGTCACGACCACGTCCATTTCCTGGCCAAGGTCCTTCTCGCCGTCGACCAGCTGACGG  >  NC_002947/1467156‑1467445

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