Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 133,201 T→C 7.2% A208A (GCT→GCC)  dsbA → thiol:disulfide interchange protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947133,2010TC 7.2% 154.2 / 2.3 55A208A (GCT→GCC) dsbAthiol:disulfide interchange protein
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base C (2/2);  total (26/29)
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.89e-01

CAAGGAATTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGT  >  NC_002947/133060‑133347
                                                                                                                                             |                                                                                                                                                  
cAAGGAATTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAg                                                                                                                                              >  2:178800/1‑148 (MQ=255)
 aaGGAATTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAg                                                                                                                                              >  5:93912/1‑147 (MQ=255)
  aGGAATTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAg                                                                                                                                           >  2:124963/1‑149 (MQ=255)
      aTTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAggcgg                                                                                                                                       <  1:124963/149‑1 (MQ=255)
            gaagaaGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTGCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCAc                                                                                                                                 <  5:189985/149‑1 (MQ=255)
              agaagTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCAccc                                                                                                                               <  6:268830/149‑1 (MQ=255)
               gaagTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCa                                                                                                                              >  3:111536/1‑149 (MQ=255)
                  gTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGc                                                                                                                           <  1:30863/149‑1 (MQ=255)
                  gTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGc                                                                                                                           <  8:104240/149‑1 (MQ=255)
                  gTATCAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGc                                                                                                                           <  2:245221/149‑1 (MQ=255)
                       agaTCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCg                                                                                                                       <  2:175623/148‑1 (MQ=255)
                        gaTCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCt                                                                                                                     >  7:245441/1‑149 (MQ=255)
                         aTCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCtt                                                                                                                    <  6:93912/149‑1 (MQ=255)
                            aCCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCtt                                                                                                                    >  2:284128/1‑146 (MQ=255)
                             ccGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCaga                                                                                                                <  8:128726/149‑1 (MQ=255)
                                    aCCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCCGCCTACTTGgcagcgc                                                                                                                                    >  3:134841/1‑121 (MQ=255)
                                    aCCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCCGCCTACTTGgcagcgc                                                                                                                                    <  4:140022/122‑2 (MQ=255)
                                    aCCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCCGCCTACTTGgcagcgc                                                                                                                                    <  4:134841/122‑2 (MQ=255)
                                    aCCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCCGCCTACTTGgcagcgc                                                                                                                                    >  3:140022/1‑121 (MQ=255)
                                      cGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGt                                                                                                       >  3:254238/1‑149 (MQ=255)
                                       gAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTg                                                                                                      <  1:178800/149‑1 (MQ=255)
                                           cTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCAt                                                                                                  >  8:253163/1‑149 (MQ=255)
                                              gtcgtcAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTgg                                                                                               <  4:254238/149‑1 (MQ=255)
                                                   cAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTg                                                                                           >  8:123050/1‑148 (MQ=255)
                                                     aCGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCAc                                                                                        <  7:253163/149‑1 (MQ=255)
                                                          aaGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGcc                                                                                   >  7:66576/1‑149 (MQ=255)
                                                          aaGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGcc                                                                                   <  8:66576/149‑1 (MQ=255)
                                                           aGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCg                                                                                  >  8:157104/1‑149 (MQ=255)
                                                              aCCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCaga                                                                                                                >  5:358217/1‑116 (MQ=255)
                                                               ccGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGTCCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAgg                                                                              <  1:284128/149‑1 (MQ=255)
                                                               ccGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAgg                                                                              <  2:333747/149‑1 (MQ=255)
                                                                   ttCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCaa                                                                          >  1:144367/1‑149 (MQ=255)
                                                                    tCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAgacggc                                                                          >  5:33224/1‑143 (MQ=255)
                                                                     cGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCaa                                                                          <  7:289850/147‑1 (MQ=255)
                                                                        ccTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACAAGc                                                                     <  7:157104/149‑1 (MQ=255)
                                                                                        ggACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGcccc                                                     >  6:305339/1‑149 (MQ=255)
                                                                                          aCCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCgg                                                   <  4:111536/149‑1 (MQ=255)
                                                                                              gAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCt                                               >  2:227630/1‑149 (MQ=255)
                                                                                               aaGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTg                                              >  7:335351/1‑149 (MQ=255)
                                                                                                  gCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGccc                                           >  5:178849/1‑149 (MQ=255)
                                                                                                          aaTGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCAc                                                                >  2:68827/1‑120 (MQ=255)
                                                                                                          aaTGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCAc                                                                <  1:68827/120‑1 (MQ=255)
                                                                                                                gCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCAcc                                                               <  6:196501/115‑1 (MQ=255)
                                                                                                                gCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCAcc                                                               >  5:196501/1‑115 (MQ=255)
                                                                                                                    aCCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCgg                         <  2:144367/149‑1 (MQ=255)
                                                                                                                      cAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGa                                         <  6:165440/131‑1 (MQ=255)
                                                                                                                      cAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGa                                         >  5:165440/1‑131 (MQ=255)
                                                                                                                            aTCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCCTCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGc                           <  7:43867/139‑1 (MQ=255)
                                                                                                                            aTCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGc                           >  8:43867/1‑139 (MQ=255)
                                                                                                                                aCAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCa                  <  3:32360/144‑1 (MQ=255)
                                                                                                                                aCAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCa                  >  4:32360/1‑144 (MQ=255)
                                                                                                                                  aaGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCt                                               >  7:176870/1‑113 (MQ=255)
                                                                                                                                  aaGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCt                                               <  8:176870/113‑1 (MQ=255)
                                                                                                                                  aaGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCaa           <  8:335351/149‑1 (MQ=255)
                                                                                                                                           gcTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCAACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGt  <  5:305339/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CAAGGAATTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGT  >  NC_002947/133060‑133347

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