Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,849,816 T→C E117G (GAA→GGA)  ygbA ← nitrous oxide‑stimulated promoter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,849,8160TC98.6% 257.6 / ‑4.6 71E117G (GAA→GGA) ygbAnitrous oxide‑stimulated promoter family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (38/32);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGG  >  NZ_CP009273/2849732‑2849901
                                                                                    |                                                                                     
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  2:1686442/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  1:1208730/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  1:1243530/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  2:818293/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  2:1233933/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  1:851443/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  1:2061397/1‑90 (MQ=255)
   ggATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGt                                                                               >  1:342511/1‑90 (MQ=255)
   ggATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGt                                                                               >  1:437950/1‑90 (MQ=255)
   ggATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGt                                                                               >  1:2511437/1‑90 (MQ=255)
   ggATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGt                                                                               >  2:290492/1‑90 (MQ=255)
     aTGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGttt                                                                             >  1:1595934/1‑90 (MQ=255)
        aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                          >  1:1450666/1‑90 (MQ=255)
            ggCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGa                                                                      >  2:440331/1‑90 (MQ=255)
             gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGAc                                                                     <  1:1754805/90‑1 (MQ=255)
             gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGAc                                                                     <  1:987450/90‑1 (MQ=255)
             gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGAc                                                                     <  2:2128446/90‑1 (MQ=255)
             gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGAc                                                                     <  2:2014293/90‑1 (MQ=255)
             gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGAc                                                                     <  1:3013857/90‑1 (MQ=255)
                       gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTc                                                           >  2:88271/1‑90 (MQ=255)
                       gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTc                                                           >  1:1888183/1‑90 (MQ=255)
                       gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTc                                                           >  2:315224/1‑90 (MQ=255)
                           aaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                       >  2:2515446/1‑90 (MQ=255)
                            aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                      <  2:2052108/90‑1 (MQ=255)
                            aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                      <  2:1208730/90‑1 (MQ=255)
                            aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                      <  2:341490/90‑1 (MQ=255)
                            aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                      <  2:85811/90‑1 (MQ=255)
                            aTCGGCATTGATTTAACTGCAAATTGCCGGACAAATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                      <  2:1205127/90‑1 (MQ=255)
                               ggCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGtt                                                   >  1:2569984/1‑90 (MQ=255)
                                     gATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCa                                             >  2:2414702/1‑90 (MQ=255)
                                      aTTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAc                                            <  2:1861983/90‑1 (MQ=255)
                                          aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                        >  2:1323757/1‑90 (MQ=255)
                                           aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAc                                       <  2:1494526/90‑1 (MQ=255)
                                             tGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGc                                     <  1:1778804/90‑1 (MQ=255)
                                              gCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGata                                                                 <  2:2176024/61‑1 (MQ=255)
                                              gCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGata                                                                 >  1:2176024/1‑61 (MQ=255)
                                              gCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCt                                    >  1:248238/1‑90 (MQ=255)
                                               cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                                   <  1:101918/90‑1 (MQ=255)
                                               cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                                   >  1:125112/1‑90 (MQ=255)
                                                 aaTTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGt                                 <  1:440331/90‑1 (MQ=255)
                                                  aTTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                <  2:2511437/90‑1 (MQ=255)
                                                             gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat                     >  1:1157960/1‑90 (MQ=255)
                                                             gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat                     >  1:56250/1‑90 (MQ=255)
                                                             gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat                     >  1:1304807/1‑90 (MQ=255)
                                                              aTCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagatg                    <  2:1304807/90‑1 (MQ=255)
                                                                   ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa               <  2:851443/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc              <  2:2609254/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc              <  1:1243946/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc              <  1:2211008/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc              <  2:1450666/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc              <  1:2515446/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc              <  1:2714614/90‑1 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGGTGAGATGACGCACGGt           >  1:283714/1‑90 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa               >  2:2832093/1‑86 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa               <  1:2832093/86‑1 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc              <  2:624031/87‑1 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc              >  1:624031/1‑87 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt           >  1:1353167/1‑90 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt           >  2:708066/1‑90 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt           >  2:832145/1‑90 (MQ=255)
                                                                        ccGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGtt          <  1:2203944/90‑1 (MQ=255)
                                                                        ccGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGtt          <  2:1353167/90‑1 (MQ=255)
                                                                          ggCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTaa        >  2:2082330/1‑90 (MQ=255)
                                                                           gCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAg       >  2:2325921/1‑90 (MQ=255)
                                                                              tACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATc    <  2:2171951/90‑1 (MQ=255)
                                                                              tACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATc    <  2:1217712/90‑1 (MQ=255)
                                                                              tACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATc    <  1:1185355/90‑1 (MQ=255)
                                                                                cTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTaa        >  2:2309076/1‑84 (MQ=255)
                                                                                cTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTaa        <  1:2309076/84‑1 (MQ=255)
                                                                                cTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCgg  >  1:1983425/1‑90 (MQ=255)
                                                                                cTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCgg  >  2:1730666/1‑90 (MQ=255)
                                                                                    |                                                                                     
CACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGG  >  NZ_CP009273/2849732‑2849901

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