Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,546,448 C→A 7.7% intergenic (+200/+207) fdnI → / ← yddM formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,546,4480CA 7.7% 2970.4 / 34.7 964intergenic (+200/+207)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (334/527);  new base A (92/10);  total (427/537)
Fisher's exact test for biased strand distribution p-value = 1.59e-24
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 965 total aligned reads displayed.

AGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAT  >  CP009273/1546380‑1546512
                                                                    |                                                                
   ttGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCaa                                                                <  1:770351/68‑1 (MQ=37)
       aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATg                                                             >  1:33490511/1‑67 (MQ=40)
       aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATg                                                             >  1:49533666/1‑67 (MQ=40)
       aGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAAtt                                                             >  1:38119959/1‑66 (MQ=255)
          gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                                           <  1:6923347/66‑1 (MQ=40)
          gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                                           <  1:44355086/66‑1 (MQ=40)
            aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                                           >  1:24115389/1‑64 (MQ=40)
            aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCt                                                        <  1:14910580/67‑1 (MQ=40)
            aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTg                                                       <  1:56623043/68‑1 (MQ=255)
            aCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGgc                                                           >  1:41705205/1‑64 (MQ=25)
             cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                      <  1:26859915/68‑1 (MQ=255)
             cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                      <  1:34911107/68‑1 (MQ=255)
             cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                      <  1:17039687/68‑1 (MQ=255)
             cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                      >  1:32218813/1‑68 (MQ=255)
             cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                      <  1:33243787/68‑1 (MQ=255)
             cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                      <  1:55683400/68‑1 (MQ=255)
             cGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                      <  1:50043980/68‑1 (MQ=255)
                   ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAg                                                 >  1:1344023/1‑67 (MQ=255)
                    gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                                               >  1:5628492/1‑68 (MQ=255)
                    gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGc                                               >  1:10028461/1‑68 (MQ=255)
                    gCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGc                                               >  1:53218741/1‑68 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGCg                                              >  1:39678796/1‑68 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGc                                               <  1:43284476/67‑1 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                              >  1:55112262/1‑68 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                              >  1:57128687/1‑68 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                              >  1:23582995/1‑68 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                              >  1:45024119/1‑68 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                              >  1:30539472/1‑68 (MQ=255)
                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                              >  1:49173074/1‑68 (MQ=255)
                        gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCa                                           <  1:2236408/68‑1 (MQ=255)
                        gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCa                                           <  1:31791369/68‑1 (MQ=255)
                          tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtt                                         <  1:405315/68‑1 (MQ=255)
                          tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtt                                         <  1:8840755/68‑1 (MQ=255)
                          tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtt                                         <  1:53757490/68‑1 (MQ=255)
                          tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtt                                         <  1:24369853/68‑1 (MQ=255)
                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                       <  1:38618695/68‑1 (MQ=255)
                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                       <  1:3392316/68‑1 (MQ=255)
                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                       <  1:42897150/68‑1 (MQ=255)
                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                       <  1:6079354/68‑1 (MQ=255)
                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAttt                                        >  1:54044619/1‑67 (MQ=255)
                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCAttt                                        >  1:47780424/1‑67 (MQ=255)
                              cGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCa                                     >  1:11764305/1‑68 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                      <  1:8429635/66‑1 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                      <  1:44709521/66‑1 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                    >  1:29521266/1‑68 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                    >  1:41182306/1‑68 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                    >  1:46058514/1‑68 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                    <  1:25608980/68‑1 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                    >  1:21286058/1‑68 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCAATTTCat                                    >  1:26259370/1‑68 (MQ=255)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGACATTTTAat                                    >  1:45582631/1‑68 (MQ=18)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCGTTAAAAGGACACAttttatt                                    >  1:57632421/1‑65 (MQ=14)
                               gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGATTAAAATGGCGCTGTAAAAGGCGTCATTTTCat                                    >  1:7241693/1‑68 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:5801656/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:46483950/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:51184167/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:50562689/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:1063743/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:18559020/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:24849993/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:34312475/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:35764307/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:19452272/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:36220388/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:18990550/68‑1 (MQ=255)
                                 aGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCatat                                  <  1:1118410/68‑1 (MQ=255)
                                  gcgcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGAGTCATTTTAATAtt                                 >  1:15397327/1‑68 (MQ=17)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                               >  1:17323239/1‑68 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                               <  1:53016259/68‑1 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                               >  1:31540213/1‑68 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                               >  1:39032035/1‑68 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                               <  1:47226789/68‑1 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                               >  1:51570560/1‑68 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                               <  1:27644157/68‑1 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCGGTAAAAGGCGTCATTTTCATATTGt                               <  1:964887/68‑1 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                               >  1:27341169/1‑68 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGt                               >  1:36862107/1‑68 (MQ=255)
                                    gcATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATCTCTCTGTAAAAGGCGTCAAATTGATATTGt                               >  1:28146618/1‑68 (MQ=17)
                                     cATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGt                                                      <  1:12486317/44‑1 (MQ=255)
                                           gCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                        >  1:20010702/1‑68 (MQ=255)
                                            cAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACaa                         <  1:1794967/66‑1 (MQ=255)
                                            cAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                        <  1:22698774/67‑1 (MQ=255)
                                            cAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAAc                        <  1:16568046/67‑1 (MQ=255)
                                            cAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACa                          <  1:31034800/65‑1 (MQ=255)
                                            aaaTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACaa                         <  1:43818479/65‑1 (MQ=255)
                                              aTGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTa                     >  1:32731553/1‑68 (MQ=255)
                                                gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAg                    >  1:28893684/1‑67 (MQ=255)
                                                  tGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                  <  1:15966905/67‑1 (MQ=255)
                                                  tGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                  <  1:5240790/67‑1 (MQ=255)
                                                  tGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                  <  1:54355956/67‑1 (MQ=255)
                                                  tGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                  <  1:29967990/67‑1 (MQ=255)
                                                  tGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                  <  1:35306166/67‑1 (MQ=255)
                                                  tGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGc                  <  1:37328487/67‑1 (MQ=255)
                                                     gTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTg              >  1:20824019/1‑68 (MQ=255)
                                                     gTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTg              >  1:14245984/1‑68 (MQ=255)
                                                     gTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTg              >  1:42488408/1‑68 (MQ=255)
                                                           catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATg        <  1:8990436/68‑1 (MQ=255)
                                                            atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGc       <  1:21784755/68‑1 (MQ=255)
                                                                gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGa   >  1:13729927/1‑68 (MQ=255)
                                                                  ttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGa   <  1:48245552/66‑1 (MQ=40)
                                                                    |                                                                
AGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAT  >  CP009273/1546380‑1546512

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