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,8160TC96.8% 194.8 / ‑1.1 63E117G (GAA→GGA) ygbAnitrous oxide‑stimulated promoter family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (30/31);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTT  >  NZ_CP009273/2849730‑2849893
                                                                                      |                                                                             
gTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCAt                                                                            <  2:497437/90‑1 (MQ=255)
  cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                          >  2:794018/1‑90 (MQ=255)
  cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                          >  1:1210493/1‑90 (MQ=255)
       aTGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGttt                                                                     >  1:291698/1‑90 (MQ=255)
          aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                  >  2:702748/1‑90 (MQ=255)
          aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                  >  1:778434/1‑90 (MQ=255)
          aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                  >  2:1475104/1‑90 (MQ=255)
          aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                  >  1:1422583/1‑90 (MQ=255)
           aGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTc                                                                 <  2:205369/90‑1 (MQ=255)
              ggCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGa                                                              >  2:542323/1‑90 (MQ=255)
              ggCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGa                                                              >  2:324815/1‑90 (MQ=255)
              ggCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGa                                                              <  1:917292/90‑1 (MQ=255)
                   ggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGata                                                         <  1:1221436/90‑1 (MQ=255)
                         gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTc                                                   >  2:714473/1‑90 (MQ=255)
                         gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTc                                                   >  1:1168854/1‑90 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              <  1:1390681/90‑1 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              <  1:1364311/90‑1 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              >  2:829156/1‑90 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              <  1:1012615/90‑1 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              <  2:970163/90‑1 (MQ=255)
                                      tGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGc                                      <  2:1104519/90‑1 (MQ=255)
                                        aTTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAc                                    <  2:657999/90‑1 (MQ=255)
                                           tAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCgg                                 >  1:847592/1‑90 (MQ=255)
                                            aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                >  1:421588/1‑90 (MQ=255)
                                            aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                >  2:492500/1‑90 (MQ=255)
                                            aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                >  1:54148/1‑90 (MQ=255)
                                            aaCTGCAAATTGCCAGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                <  1:1512811/90‑1 (MQ=255)
                                             aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAc                               <  1:829927/90‑1 (MQ=255)
                                             aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAc                               <  2:778434/90‑1 (MQ=255)
                                             aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAc                               <  2:181819/90‑1 (MQ=255)
                                                gCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCt                            >  1:819992/1‑90 (MQ=255)
                                                gCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCAATGAGGTTTTTTCGGACGATATGTTTCCGGCAGTTCTGGCACCGGACGCt                            >  2:1190817/1‑90 (MQ=255)
                                                 cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                           <  1:790282/90‑1 (MQ=255)
                                                 cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                           <  2:1210493/90‑1 (MQ=255)
                                                       gCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATc                     >  1:1375415/1‑90 (MQ=255)
                                                            aCAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatga                <  1:542323/90‑1 (MQ=255)
                                                            aCAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatga                <  1:324815/90‑1 (MQ=255)
                                                               gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat             >  2:648483/1‑90 (MQ=255)
                                                               gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat             >  1:677922/1‑90 (MQ=255)
                                                                    gCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc        <  1:702748/90‑1 (MQ=255)
                                                                    gCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc        >  1:745703/1‑90 (MQ=255)
                                                                    gCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc        <  1:601318/90‑1 (MQ=255)
                                                                     ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa       <  2:1449673/90‑1 (MQ=255)
                                                                     ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa       <  1:273211/90‑1 (MQ=255)
                                                                     ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa       <  2:452736/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  2:568372/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  1:1119859/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  2:847592/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  1:894295/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  1:714473/90‑1 (MQ=255)
                                                                        gTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACgg    >  1:574768/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                 <  2:1163428/76‑1 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                 >  1:1163428/1‑76 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  2:209778/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  2:1472290/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  2:1281995/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  2:1276391/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  1:143117/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  1:1189130/1‑90 (MQ=255)
                                                                          ccGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGtt  <  2:1442250/90‑1 (MQ=255)
                                                                          ccGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGtt  <  1:794018/90‑1 (MQ=255)
                                                                                  cTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  2:466306/78‑1 (MQ=255)
                                                                                  cTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      >  1:466306/1‑78 (MQ=255)
                                                                                      |                                                                             
GTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTT  >  NZ_CP009273/2849730‑2849893

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