Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,159,046 T→G L186R (CTG→CGG)  oxyR → oxidative and nitrosative stress transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,159,0460TG95.9% 158.5 / ‑2.3 50L186R (CTG→CGG) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base G (25/22);  minor base A (0/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGC  >  NC_000913_3_bme_pgi/4158903‑4159184
                                                                                                                                               |                                                                                                                                          
cGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAgg                                                                                                                                      <  1:99987/150‑1 (MQ=255)
cGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAgg                                                                                                                                      >  2:7380/1‑150 (MQ=255)
cGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAgg                                                                                                                                      <  1:541295/150‑1 (MQ=255)
 ggCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAggg                                                                                                                                     >  1:63964/1‑150 (MQ=255)
  gCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAgggg                                                                                                                                    >  2:294650/1‑150 (MQ=255)
  gCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATCCGTACCGACGGCCGATCGGCCAgggg                                                                                                                                    >  2:271437/1‑150 (MQ=255)
       cTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGaaaaa                                                                                                                               >  1:275252/1‑150 (MQ=255)
                gTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATg                                                                                                                      >  2:210461/1‑150 (MQ=255)
                     ccTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTCGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGa                                                                                                                 <  2:53108/150‑1 (MQ=255)
                          cgcTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATg                                                                                                            <  2:252539/150‑1 (MQ=255)
                                 aaaaGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGACAAGGTTGCTGGCTATCTATGAAGATCACCCGTGAGCGAACCGCGAATGCGTACCGATGGCCGAAAAAAAAAAAAAAAAACTGCTGATGCTGGAAGATGGTCACTg                                                                                                     <  4:316511/149‑1 (MQ=255)
                                  aaaGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGt                                                                                                    >  1:99948/1‑150 (MQ=255)
                                  aaaGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGt                                                                                                    <  2:197926/150‑1 (MQ=255)
                                    agagagCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGttt                                                                                                  <  3:42748/150‑1 (MQ=255)
                                          cGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGa                                                                                            >  1:290946/1‑150 (MQ=255)
                                                    attGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGAt                                                                                           >  1:310572/1‑141 (MQ=255)
                                                    attGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGAt                                                                                           <  2:310572/141‑1 (MQ=255)
                                                    attGAAGTGCCGGTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATg                                                                                  >  4:142695/1‑150 (MQ=255)
                                                         aGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGttt                                                                             <  2:63964/150‑1 (MQ=255)
                                                         aGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCAGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGttt                                                                             <  1:29391/150‑1 (MQ=255)
                                                                   tttGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGaa                                                                   <  2:275252/150‑1 (MQ=255)
                                                                     tgatgaGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCt                                                                           >  4:101706/1‑140 (MQ=255)
                                                                     tgatgaGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCt                                                                           <  3:101706/140‑1 (MQ=255)
                                                                        tgaGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCgg                                                              <  2:99948/150‑1 (MQ=255)
                                                                         gaGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCggg                                                             <  1:7380/150‑1 (MQ=255)
                                                                             cAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCgg                                                         <  1:294650/150‑1 (MQ=255)
                                                                                  ttGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGaa                                                    >  1:252834/1‑150 (MQ=255)
                                                                                     cTGGCTATCTATGAAGATCGCCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGGTCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGAt                                                 >  4:51287/1‑150 (MQ=255)
                                                                                       ggCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCgg                                                         >  2:3046/1‑140 (MQ=255)
                                                                                       ggCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCgg                                                         <  1:3046/140‑1 (MQ=255)
                                                                                          tatcttgGAAGATGACACGTGGGCGAACCGCGAATGTGTACCGATGTTCGATCTGTCAGTGGAAAAACTGATGATGCTGGACGAAGGTAACTGTTTGCGCGATCAGGCAATGGGTTTCATTTTTTAAGCCCGTGCGGATGAAGATacaca                                            >  4:153070/1‑150 (MQ=255)
                                                                                          tatctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATacaca                                            <  2:271767/150‑1 (MQ=255)
                                                                                           atctatGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATacacac                                           <  1:210461/150‑1 (MQ=255)
                                                                                               atGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAActg                                                                                                                            <  1:15585/65‑1 (MQ=255)
                                                                                               atGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAActg                                                                                                                            >  2:15585/1‑65 (MQ=255)
                                                                                               atGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTCTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTcc                                       >  1:520676/1‑150 (MQ=255)
                                                                                                tGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCcg                                      >  4:295109/1‑150 (MQ=255)
                                                                                                 gAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCcgc                                     >  1:37071/1‑150 (MQ=255)
                                                                                                             tGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTg                         <  3:153070/150‑1 (MQ=255)
                                                                                                             tGGGCGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTg                         <  2:252834/150‑1 (MQ=255)
                                                                                                                 cGAACCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGaaa                     >  1:433675/1‑150 (MQ=255)
                                                                                                                    aCCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAActgct                                                                                                                          >  2:208635/1‑46 (MQ=255)
                                                                                                                    aCCGCGAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAActgct                                                                                                                          <  1:208635/46‑1 (MQ=255)
                                                                                                                      cgcgAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTg                >  4:206011/1‑150 (MQ=255)
                                                                                                                       gcgAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGcgcg                                                                                             >  1:72514/1‑72 (MQ=255)
                                                                                                                       gcgAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGcgcg                                                                                             <  2:72514/72‑1 (MQ=255)
                                                                                                                       gcgAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGAc                                >  3:41075/1‑133 (MQ=255)
                                                                                                                       gcgAATGCGTACCGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGAc                                <  4:41075/133‑1 (MQ=255)
                                                                                                                                   cGATGGCCGATCGGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACAtggtgg   >  3:356892/1‑150 (MQ=255)
                                                                                                                                    gATGGCCGATCGGGCAGGGGAAAAACTGCTGATGATGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTggc  >  4:227213/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                          
CGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGC  >  NC_000913_3_bme_pgi/4158903‑4159184

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