Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,676,832 G→T intergenic (‑25/+158) dctA ← / ← yhjK C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,676,8320GT98.2% 169.5 / NA 57intergenic (‑25/+158)dctA/yhjKC4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (28/27);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATC  >  CP009273/3676686‑3676977
                                                                                                                                                  |                                                                                                                                                 
aaGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGt                                                                                                                                                 <  1:383980/149‑1 (MQ=255)
   ccGTCGCCAAGCGGTTTCATTTGCTCGCCTATTACAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTCTTAAACAGAGAGGTTTTCATAGGGTGTCCTGTAGTAAAACCACATGTCtt                                                                                                                                              >  1:324596/1‑149 (MQ=255)
        gCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAg                                                                                                                                         <  2:264663/149‑1 (MQ=255)
          cAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGACAGCAAAGGCATTTGAGAAAACAGGCTTTCATGGCGTGTCCTTTAGTAAAACCACATGCCTTGTAAGGt                                                                                                                                       >  1:442185/1‑149 (MQ=255)
            aGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCCTTAGTAAAACCACATATCTTGTAAGGtca                                                                                                                                     >  2:37050/1‑147 (MQ=255)
            aGCGGTTTCATTTGCTCGCCCATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTa                                                                                                                                     <  2:603873/149‑1 (MQ=255)
              cGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATg                                                                                                                                   <  2:479764‑M1/149‑2 (MQ=255)
              cGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATg                                                                                                                                   <  1:323740‑M1/149‑2 (MQ=255)
              cGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCCATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATg                                                                                                                                   <  1:510422‑M1/149‑2 (MQ=255)
                gTTTCATTTGCTCGCCGATTTGAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggg                                                                                                                                 <  1:754141‑M1/149‑4 (MQ=255)
                    cATTTGCTCGCCTATTTCAGGATAGGAATGGGCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATTTCTTGTAAGGTTATggggtac                                                                                                                             <  2:953978‑M1/149‑8 (MQ=255)
                     aTTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGGTTTCATAGGGTGTCCGTTAGTAAAACGACAGGTCTTGTAAGGTTATggggtacc                                                                                                                            >  2:985057‑M1/1‑141 (MQ=255)
                     aTTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATAGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACACAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACTACATGTCTTGTAAGGTTATggggtacc                                                                                                                            >  2:985034‑M1/1‑141 (MQ=255)
                          cTCGCCTATTTCAGTATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcgg                                                                                                                       <  1:769971‑M1/149‑14 (MQ=255)
                          cTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcgg                                                                                                                       <  1:100987‑M1/149‑14 (MQ=255)
                           tCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggg                                                                                                                      <  2:314372‑M1/149‑15 (MQ=255)
                            cGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAGCCACATGTCTTGTAAGGTTATggggtacctgcgggg                                                                                                                     <  2:422534‑M1/149‑16 (MQ=255)
                                    tCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcg                                                                                                             <  2:379832‑M1/149‑24 (MQ=255)
                                     cAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGGTATggggtacctgcgggggacctgggg                                                                                                            >  1:157840‑M1/1‑125 (MQ=255)
                                      cggaTAGAGATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCAAATGTTTTAATTTTTTAGttggtacctgcggggtacctgcggc                                                                                                           <  1:243983‑M1/148‑26 (MQ=255)
                                      aGGATAGAAAGGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggc                                                                                                           <  1:130405‑M1/149‑26 (MQ=255)
                                          tAGAAATGGCCAAGGAGAATACCAATGGCTATCGGTGGCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGGTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggcctta                                                                                                       >  2:92532‑M1/1‑120 (MQ=255)
                                             aaaTGGCCAAGGAGAATACCAATGGATATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGAGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaa                                                                                                    <  1:616249‑M1/149‑33 (MQ=255)
                                                tGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataa                                                                                                 <  1:454908‑M1/149‑36 (MQ=255)
                                                tGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataa                                                                                                 >  1:402891‑M1/1‑114 (MQ=255)
                                                    cAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCCTAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacc                                                                                             >  1:140203‑M1/1‑110 (MQ=255)
                                                     aaGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggagaacatgcgagttaccttgggcctacccaacacactc                                                                                            >  1:229861‑M1/1‑109 (MQ=255)
                                                     aaGGAGAATACCAATGGCTATCGCTGGCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtaccggcggccttaaaataacaccc                                                                                            >  1:173264‑M1/1‑109 (MQ=255)
                                                                cAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatca                                                                                 <  1:414014‑M1/149‑52 (MQ=255)
                                                                  aTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcaca                                                                               <  1:713491‑M1/149‑54 (MQ=255)
                                                                       tATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGGGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaatgacagccagacaacatccaagaaat                                                                          >  1:419638‑M1/1‑91 (MQ=255)
                                                                       tATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGATTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcgggggacctgcggccttaaaatgacaccccgacaacatcacaggaat                                                                          >  1:576895‑M1/1‑91 (MQ=255)
                                                                         tCGCTGTAAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcccagaaatgt                                                                        >  2:263691‑M1/1‑89 (MQ=255)
                                                                                 aGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATgggatacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggat                                                                >  2:135687‑M1/1‑81 (MQ=255)
                                                                                   gACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTCTggggtacctgcggggtaccctcgtccctaaaataacacccagacaacatcacagaaatgtacctggatca                                                              >  1:935356‑M1/1‑79 (MQ=255)
                                                                                      ctggaggtaaaggCGGTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACTTGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataa                                                           >  2:739406‑M1/7‑76 (MQ=255)
                                                                                      cTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataa                                                           >  2:436889‑M1/1‑76 (MQ=255)
                                                                                       tGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTCGTAACACTACATGTCTTGTAAGGTTATgcggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataat                                                          <  1:741688‑M1/149‑75 (MQ=255)
                                                                                        gAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgcacctggatcataatt                                                         >  2:994630‑M1/1‑74 (MQ=255)
                                                                                            gTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaa                                                     <  1:869479‑M1/149‑80 (MQ=255)
                                                                                            gTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaa                                                     <  1:823850‑M1/149‑80 (MQ=255)
                                                                                                  gCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaa                                               >  1:419448‑M1/1‑64 (MQ=255)
                                                                                                        aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagcaatgtacccggatcataattgaaacaacaaggttaa                                         >  2:428409‑M1/1‑58 (MQ=255)
                                                                                                        aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaacgttaa                                         >  2:428322‑M1/1‑58 (MQ=255)
                                                                                                         aaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcatacttgaaacaaaaaggttaaa                                        >  2:626274‑M1/1‑57 (MQ=255)
                                                                                                           cAGAGAGGTTTTCATGGGGGGGCCTTTAGTAAAACCACATGTCTTGTAAGGTTATgggttacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataatttaaacaaaaaggttaaaaa                                      <  2:271810‑M1/149‑95 (MQ=255)
                                                                                                               gagGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagttt                                  <  1:436889‑M1/149‑99 (MQ=255)
                                                                                                                      tCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtcacaggtttatatttttaaccaaaacgttaaaaagtttgaactgg                           >  1:931006‑M1/1‑44 (MQ=255)
                                                                                                                        aTAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggct                         >  1:660760‑M1/1‑42 (MQ=255)
                                                                                                                         tAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctc                        <  1:37050‑M1/149‑109 (MQ=255)
                                                                                                                                  cTTTAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggctcataattgaaacaaaaaggttaaaaagtttgagctggctcggacaaacc               >  1:514362‑M1/1‑32 (MQ=255)
                                                                                                                                   tttAGTAAAACGACATGTCTTGTAACGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctcgcacaaacca              <  2:543955‑M1/149‑119 (MQ=255)
                                                                                                                                   tttAGTAAAACCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctcgcacaaacca              >  2:209825‑M1/1‑31 (MQ=255)
                                                                                                                                        tAAAACCAAATGGCTTGTAGGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctcgcacaaaccagcact         <  2:140203‑M1/149‑124 (MQ=255)
                                                                                                                                          aaaCCACATGTCTTGTAAGGTTATggggtacctgcggggtacctgcggccttaaaataacacccagacaacatcacagaaatgtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctcgcacaaaacagcacttt       >  1:176890‑M1/1‑24 (MQ=255)
                                                                                                                                             ccACACGTCTTTTAAGGTTATggggtacctgcggtgtacctccggccttaaaataacacccagacaacatcacataaaggtacctggatcataattgaaacaaaaaggttaaaaagtttgagctggctctcacaaaccagcacttttta    <  2:886742‑M1/149‑129 (MQ=255)
                                                                                                                                               acaTGTCTTGTAAGGTTATggggtacgtgcggggtacctgcggccttaaaataacacccagacaacagcacagaaatgtaccgggatcataattgaaacaaaaaggttaaaaagtttgggctggctcgcacaaaccagccctttttaaa  >  2:468635‑M1/1‑19 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATC  >  CP009273/3676686‑3676977

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