Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,519,318 G→C R300G (CGG→GGG)  yfeR ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,519,3180GC86.8% 160.0 / 19.6 53R300G (CGG→GGG) yfeRLysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base C (23/23);  total (26/27)
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.12e-01

CGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGTAATGCGCTTTACCACCAG  >  NZ_CP009273/2519183‑2519453
                                                                                                                                       |                                                                                                                                       
cGGATGACAAACTGGGCGAAGGTTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACATCGTATCCCCTTcccccc                                                                                                                                       >  2:635532/1‑138 (MQ=255)
cGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTcccccc                                                                                                                                       >  2:763277/1‑138 (MQ=255)
 ggATTACAAACTGGGCTAAGATTTTTCCCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCa                                                                                                                                     <  2:864711/139‑1 (MQ=255)
  gATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAg                                                                                                                                    >  1:974766/1‑139 (MQ=255)
     gACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCa                                                                                                                                 >  1:588728/1‑139 (MQ=255)
      aCAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCa                                                                                                                                 <  1:48384/138‑1 (MQ=255)
       cAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCAtt                                                                                                                               <  1:98681/139‑1 (MQ=255)
             gggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCa                                                                                                                         >  1:1003200/1‑139 (MQ=255)
              ggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCa                                                                                                                         <  2:703743/138‑1 (MQ=255)
              ggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCa                                                                                                                         <  2:588728/138‑1 (MQ=255)
              ggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCa                                                                                                                         <  1:895704/138‑1 (MQ=255)
               gCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCAtt                                                                                                                       >  2:570114/1‑139 (MQ=255)
                   aGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGccc                                                                                                                   <  1:734950/139‑1 (MQ=255)
                   aGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGccc                                                                                                                   <  1:73940/139‑1 (MQ=255)
                   aGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGccc                                                                                                                   <  1:1257177/139‑1 (MQ=255)
                        tttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGt                                                                                                               <  2:884215/138‑1 (MQ=255)
                          nncnAGGNTGCGATAATCTTCTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTg                                                                                                                 <  2:412162/132‑1 (MQ=255)
                                  tGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACg                                                                                                    <  2:992879/139‑1 (MQ=255)
                                              cTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCa                                                                                                                         <  1:61277/106‑1 (MQ=255)
                                              cTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCa                                                                                                                         >  2:61277/1‑106 (MQ=255)
                                              cTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCt                                                                                        >  2:168238/1‑139 (MQ=255)
                                               tCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                                                                       >  1:911291/1‑139 (MQ=255)
                                               tCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                                                                       >  2:646667/1‑139 (MQ=255)
                                                 cATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                                                                       >  2:955312/1‑137 (MQ=255)
                                                   tCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgc                                                                                   >  1:883042/1‑139 (MQ=255)
                                                      aaTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcggc                                                                                <  1:168238/139‑1 (MQ=255)
                                                           aaCTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTg                                                                                                                 <  1:1189317/101‑1 (MQ=255)
                                                           aaCTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTg                                                                                                                 >  2:1189317/1‑101 (MQ=255)
                                                                   attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGaa                                                                   >  2:355593/1‑139 (MQ=255)
                                                                   attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                                                                    <  1:763277/138‑1 (MQ=255)
                                                                   attgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                                                                    <  1:955312/138‑1 (MQ=255)
                                                                    ttgattACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAAc                                                                  <  1:1048927/139‑1 (MQ=255)
                                                                    ttgattACTACCGGGGTTCGCGGCAGTATGAATTCACTCTATATCTGATACAACGGCTCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCTTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAAc                                                                  <  2:930744/139‑1 (MQ=255)
                                                                                 gggCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg                                                     >  1:933374/1‑139 (MQ=255)
                                                                                 gggCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcg                                                     >  2:1032452/1‑139 (MQ=255)
                                                                                    cTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCAc                                                  <  2:314625/139‑1 (MQ=255)
                                                                                                        ctctATATCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTca                               >  2:1209641/1‑138 (MQ=255)
                                                                                                               tCTGATACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCg                        <  2:858467/138‑1 (MQ=255)
                                                                                                                   aTACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGt                                                               >  1:607715/1‑95 (MQ=255)
                                                                                                                   aTACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGt                                                               <  2:607715/95‑1 (MQ=255)
                                                                                                                   aTACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGt                    >  1:1214939/1‑138 (MQ=255)
                                                                                                                   aTACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGt                    >  2:820840/1‑138 (MQ=255)
                                                                                                                   aTACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGt                    >  1:540338/1‑138 (MQ=255)
                                                                                                                   aTACAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACCGGCGTTCAACAACCGGCGt                    >  2:1070741/1‑138 (MQ=255)
                                                                                                                     aCAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGTAAt                 >  2:204914/1‑139 (MQ=255)
                                                                                                                       aaCGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGTAATgc               <  1:1242137/139‑1 (MQ=255)
                                                                                                                       aaCGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGTAATgc               >  2:330454/1‑139 (MQ=255)
                                                                                                                         cGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCCGCATCAACTGGCGTTCAACAACCGGCGTAATgcgc             >  2:762572/1‑139 (MQ=255)
                                                                                                                         cGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcgg                                                                                 <  1:302781/71‑1 (MQ=255)
                                                                                                                         cGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcgg                                                                                 >  2:302781/1‑71 (MQ=255)
                                                                                                                         cGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCa                                                   <  1:219601/101‑1 (MQ=255)
                                                                                                                         cGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCa                                                   >  2:219601/1‑101 (MQ=255)
                                                                                                                               cccTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGTAATGCGCTTTac        <  1:330454/138‑1 (MQ=255)
                                                                                                                                    ccccccCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGTAATGCGCTTTACCACCAg  >  1:946780/1‑139 (MQ=255)
                                                                                                                                    cccGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCa                                                   <  2:263404/90‑1 (MQ=255)
                                                                                                                                    cccGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCa                                                   >  1:263404/1‑90 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAACTGGCGTTCAACAACCGGCGTAATGCGCTTTACCACCAG  >  NZ_CP009273/2519183‑2519453

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: