Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 3,698,796 C→T 9.0% P195S (CCC→TCC)  bcsG → DUF3260 family cellulose production inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA3,698,7960CT 9.0% 133.4 / 4.7 54P195S (CCC→TCC) bcsGDUF3260 family cellulose production inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base T (3/2);  total (31/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.14e-01

CTGCTATGGCTGAACGTACTTACCCTGGCGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTC  >  NC_000913_3_pae_tpiA/3698649‑3698944
                                                                                                                                                   |                                                                                                                                                    
cTGCTATGGCTGAACGTACTTACCCTGGCGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGcc                                                                                                                                                     <  2:720153/149‑1 (MQ=255)
 tGCTATGGCTGAACGTACTTACCCTGGCGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGccc                                                                                                                                                    <  1:315033/149‑1 (MQ=255)
            aaCGTACTTACCCTGGAGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGc                                                                                                                                         >  1:208337/1‑149 (MQ=255)
             aCGTACTTACCCTGGCGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCAGAACCTGCa                                                                                                                                        >  2:450851/1‑149 (MQ=255)
                   ttACCCTGGCGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGCCGCACCGGTAGTGGGTGATATGCCGCCACAAACTGCACCGCca                                                                                                                                  >  2:137480/1‑149 (MQ=255)
                         tGGCGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACg                                                                                                                            >  1:255709/1‑149 (MQ=255)
                         tGGCGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACAGACGACCACTGTAACAACGACGGGTGGGAACGCAGCGGCAACCGTTCCGGCGACGGGTGGCGCACCGGTAGTGGGGGATATGTGCGCCCAGACTCCACCGCcaacacag                                                                                                                            >  2:595738/1‑146 (MQ=255)
                               gACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGCAAACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAAc                                                                                                                      <  1:416252/149‑1 (MQ=255)
                               gACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAAc                                                                                                                      <  1:479365/149‑1 (MQ=255)
                                aCCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACCAACTGCACCGCCAACAACGGCGAAcc                                                                                                                     >  1:773232/1‑149 (MQ=255)
                                    aGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTaa                                                                                                                 >  2:358584/1‑149 (MQ=255)
                                     gTTTCTCCTTGGGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAAc                                                                                                                <  1:395961/149‑1 (MQ=255)
                                          tcCTTGTGGCCAGCCGGACACCCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCCGTAGTGGGTGATATGCCCGCACAAACTGCGCCGCCAACAACGGCCAACCTTAACGCCTg                                                                                                           <  2:913433/149‑1 (MQ=255)
                                             ttGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCGGGCt                                                                                                        >  1:847987/1‑149 (MQ=255)
                                                  gCCAGCCGGCCACCCCGCGTCCACGGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTAATATGCCCGCACAAAATGCACCGCCAACAACGGAGCACCTTAAAGCCGAGCTGaata                                                                                                   >  2:353235/1‑149 (MQ=255)
                                                  gCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACAGCACCGACAACAACGACGAACCTTAACGCCTGACTGacta                                                                                                   >  1:366326/1‑149 (MQ=255)
                                                       ccGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTc                                                                                              <  1:916460/149‑1 (MQ=255)
                                                       ccGCACCACCGACGACGGGTTGAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTc                                                                                              <  2:387754/149‑1 (MQ=255)
                                                             aaCCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCt                                                                                                                    >  2:379225/1‑121 (MQ=255)
                                                             aaCCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCt                                                                                                                    <  1:379225/121‑1 (MQ=255)
                                                                cgacgaCCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAAcgcg                                                                                     <  1:123355/149‑1 (MQ=255)
                                                                  atgacCACGGTAACAACGGCGAGTGGCAACGCAGCTGCAACCGTTGACTCCCCGCGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGa                                                                                   <  1:915102/147‑1 (MQ=255)
                                                                  acgacCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGa                                                                                   <  2:310693/149‑1 (MQ=255)
                                                                      ccACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGAGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCg                                                                               >  1:780425/1‑149 (MQ=255)
                                                                       cACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGCCGGGTGGCGCACCGGTAGTGGGCGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTCGCTGAATCACTTCTATAACACCGCGGCGa                                                                              >  2:533773/1‑149 (MQ=255)
                                                                       cACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGa                                                                              >  2:491814/1‑149 (MQ=255)
                                                                         cTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGaaa                                                                            <  1:439458/149‑1 (MQ=255)
                                                                            tAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGACAACAACGGCGAACCTTAACACCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGt                                                                         >  2:257976/1‑149 (MQ=255)
                                                                              acaacGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTaa                                                                       >  2:305773/1‑149 (MQ=255)
                                                                                aacGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAAt                                                                     >  1:395842/1‑149 (MQ=255)
                                                                                 acgacgGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATc                                                                    <  2:426257/149‑1 (MQ=255)
                                                                                  cgacgGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGTCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTgg                                                                                                          <  2:655180/110‑1 (MQ=255)
                                                                                  cgacgGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGTCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTgg                                                                                                          <  2:654835/110‑1 (MQ=255)
                                                                                  cgacgGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGTCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTgg                                                                                                          >  1:654835/1‑110 (MQ=255)
                                                                                  cgacgGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGTCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTgg                                                                                                          >  1:655180/1‑110 (MQ=255)
                                                                                          tgtcaCGCAGCGCGAACCGTTGCGTCGACGGGTGGTGCACCGGTAGCGGGTGATATGCCCGCACAAACTGCATCGCCAACAACGGCGAACTTTAACGCTTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTccc                                                           <  2:319233/145‑1 (MQ=255)
                                                                                           gtAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCg                                                          >  1:533115/1‑149 (MQ=255)
                                                                                            tAACGCAGCGGCAACCGTTGCGGCGACGGGTAGCGCACCGGTAGTGGGCGGTATGCCGGCACAAACTGAACCGCCAACAACGACGAACCTTAACGCCTGGCTGACTAATTTCTATAACGCGGAGGCGAAACGTAAATGGACCTTCCCGt                                                         >  2:157006/1‑149 (MQ=255)
                                                                                                gCAGCGGCAACCGTTTCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTc                                                     >  2:367728/1‑149 (MQ=255)
                                                                                                    cGGCAACCGTTGCGGCGACGGTTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTg                                                 <  2:319396/149‑1 (MQ=255)
                                                                                                          ccGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCt                                           >  1:474619/1‑149 (MQ=255)
                                                                                                                 ggcgACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCTAAACGTAAAGCGAACTTCACGTCTTCGCTGCCCGCTAATGCta                                    >  2:189349/1‑148 (MQ=255)
                                                                                                                   ctaCGGGTGGTGAACAGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAg                                  <  2:305825/147‑1 (MQ=255)
                                                                                                                    gACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGc                                 >  1:627915/1‑149 (MQ=255)
                                                                                                                        agtgacgcACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCAtt                             >  1:796143/6‑149 (MQ=255)
                                                                                                                         gTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCAttt                            <  2:208337/149‑1 (MQ=255)
                                                                                                                             cgcACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAac                        <  1:358584/149‑1 (MQ=255)
                                                                                                                              gcACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAaca                       >  1:119836/1‑148 (MQ=255)
                                                                                                                               cACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCCGCCCGCTGATGCTCAGCCATTTGAacta                      <  1:470214/149‑1 (MQ=255)
                                                                                                                                aCCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGCTCCTCAGCCAGTAGAactcc                     >  2:216651/1‑147 (MQ=255)
                                                                                                                                 ccGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGCCTGAATAATGTCTATAACGCGGAGGCGAAACGTAAATAGACCTTCCCGTCCTCGCTGCCCGCTCATCCTCAGACAATTGAactaat                    >  2:790306/1‑147 (MQ=255)
                                                                                                                                      agggggTTATATGCGCGGTCAAATTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGGAAAGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGa               <  1:31742/146‑1 (MQ=255)
                                                                                                                                        tGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGAtt             <  1:149643/149‑1 (MQ=255)
                                                                                                                                          ggTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTAATGATTaa           >  1:417128/1‑149 (MQ=255)
                                                                                                                                            tGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACa         >  1:830631/1‑149 (MQ=255)
                                                                                                                                             gATATGCCCCTACGAGCTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGTCGAAACGTAAATCGTCCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACAt        <  2:17853/149‑1 (MQ=255)
                                                                                                                                                   cccGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTc  <  2:255709/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGCTATGGCTGAACGTACTTACCCTGGCGGGACCAAGTTTCTCCTTGTGGCCAGCCGGACAACCGACGACCACTGTAACAACGACGGGTGGTAACGCAGCGGCAACCGTTGCGGCGACGGGTGGCGCACCGGTAGTGGGTGATATGCCCGCACAAACTGCACCGCCAACAACGGCGAACCTTAACGCCTGGCTGAATAATTTCTATAACGCGGAGGCGAAACGTAAATCGACCTTCCCGTCTTCGCTGCCCGCTGATGCTCAGCCATTTGAACTACTGGTGATTAACATCTGTTC  >  NC_000913_3_pae_tpiA/3698649‑3698944

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