Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,850,204 C→G 14.3% intergenic (+76/‑111) hemN → / → fnrA oxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,850,2040CG14.3% 150.3 / 17.5 56intergenic (+76/‑111)hemN/fnrAoxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base G (3/5);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

GCGCTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGCC  >  NC_002947/4850056‑4850353
                                                                                                                                                    |                                                                                                                                                     
gcgcTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCa                                                                                                                                                      >  2:671291/1‑150 (MQ=255)
  gcTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCAcc                                                                                                                                                    >  1:534460/1‑150 (MQ=255)
   cTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCt                                                                                                                                                   >  2:805087/1‑150 (MQ=255)
          tGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            <  2:580444/150‑2 (MQ=255)
          tGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            <  1:816195/150‑2 (MQ=255)
          tGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            <  1:1158751/150‑2 (MQ=255)
                      gACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAgcgc                                                                                                                                >  1:898028/1‑150 (MQ=255)
                       aCGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAgtgtg                                                                                                                               >  1:472717/1‑146 (MQ=255)
                              tCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            <  2:399344/130‑2 (MQ=255)
                              tCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            >  1:399344/1‑129 (MQ=255)
                                          caACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACt                                                                                                            >  1:87456/1‑150 (MQ=255)
                                           aaCCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTg                                                                                                           >  1:856421/1‑150 (MQ=255)
                                            aCCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTgg                                                                                                          <  1:671291/150‑1 (MQ=255)
                                               gCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCt                                                                                                       <  2:534460/150‑1 (MQ=255)
                                                 cAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTac                                                                                                     <  1:39791/150‑1 (MQ=255)
                                                         cgcgcgAGTGATCTGAAAAACAAGCCCTTCAACCGCATGGACAAGCTTCCTTGTCATGCACACTATGAGAACACAAGGCCTCCAGGCAGGCCACCTGCCGGGGTTGCGCCAGCGCGCACCATCACCGTGAACACTGGCCTACACCCtgtg                                                                                             <  1:527865/148‑1 (MQ=255)
                                                          tcGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGcc                                                                                                                                                       <  1:3784/91‑1 (MQ=255)
                                                          tcGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGcc                                                                                                                                                       >  2:3784/1‑91 (MQ=255)
                                                          tcGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTgtgt                                                                                            >  2:1196705/1‑150 (MQ=255)
                                                           cgcgAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGATCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            <  1:446805/101‑2 (MQ=255)
                                                           cgcgAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGATCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                            >  2:446805/1‑100 (MQ=255)
                                                               aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGcc                                                                                                                                                       >  1:36608/1‑86 (MQ=255)
                                                               aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGcc                                                                                                                                                       <  2:36608/86‑1 (MQ=255)
                                                               aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                                       <  2:856421/150‑1 (MQ=255)
                                                                 tGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGt                                                                                     <  1:1196705/150‑1 (MQ=255)
                                                                      tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGca                                                                                                                >  2:400278/1‑118 (MQ=255)
                                                                      tGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGca                                                                                                                <  1:400278/118‑1 (MQ=255)
                                                                       gAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCt                                                                               >  1:671188/1‑150 (MQ=255)
                                                                              aaGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACt                                                                        >  2:711741/1‑150 (MQ=255)
                                                                                ggCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTa                                                                      <  1:805087/150‑1 (MQ=255)
                                                                                    ttCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAtgtg                                                                  >  1:289220/1‑150 (MQ=255)
                                                                                            cATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTa                                                                                                      >  1:720464/1‑106 (MQ=255)
                                                                                            cATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTa                                                                                                      <  2:720464/106‑1 (MQ=255)
                                                                                                 aCAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGCCTTATGTGTGCTTTCCCACaa                                                     <  2:898028/150‑1 (MQ=255)
                                                                                                  cAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAg                                                    <  2:87456/150‑1 (MQ=255)
                                                                                                   aaGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAgg                                                   >  1:841916/1‑150 (MQ=255)
                                                                                                          cTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGGCGGCGTTCCCACAAGGATCGAtt                                            >  2:280308/1‑150 (MQ=255)
                                                                                                              tCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAgcgc                                                                                                                                <  1:152740/62‑1 (MQ=255)
                                                                                                              tCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAgcgc                                                                                                                                >  2:152740/1‑62 (MQ=255)
                                                                                                               cAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAAt                                       <  1:621999/150‑1 (MQ=255)
                                                                                                               cAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTCGCCTACACCCTTTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAAt                                       <  1:333155/150‑1 (MQ=255)
                                                                                                                        tATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCAAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGcc                              >  2:462145/1‑150 (MQ=255)
                                                                                                                         aTGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGAt                                                 <  2:60426/130‑1 (MQ=255)
                                                                                                                         aTGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGAt                                                 >  1:60426/1‑130 (MQ=255)
                                                                                                                         aTGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGGTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTAACCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTAGAATGTCCGAGCCa                             >  1:1154450/1‑150 (MQ=255)
                                                                                                                                 ccaAGGTATCCCGACAGGACCACTCTCGCGGTTTCCCCAGCGCGCACCATCTGCTAGCACAGCGGCCTACACCCTGTTACGTGAGTAACCCCTACAACGTATGAGTGGTTTCCCACTAGCCTTGTTTCAAATGTCCGAGCCAGTCAAACt                     <  2:175194/149‑1 (MQ=255)
                                                                                                                                  caAGGACTCCCGGCCGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACGCTTACGACTTATGTGTGCTTTCCCACAAAGATCGATTCAAATGTCCGCGCCAGTCAAACTg                    >  2:81435/1‑150 (MQ=255)
                                                                                                                                   aaGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTgc                   >  2:724473/1‑150 (MQ=255)
                                                                                                                                     gaaaccccGGCAGGTGGCCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTgcgc                 >  2:184724/6‑150 (MQ=255)
                                                                                                                                         tCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGCGCCAATCAAGATGCGCCcac             >  2:974909/1‑150 (MQ=255)
                                                                                                                                           ccGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCcacac           >  1:212000/1‑150 (MQ=255)
                                                                                                                                            cGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGAt                                             <  2:441019/115‑1 (MQ=255)
                                                                                                                                            cGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGAt                                             >  1:441019/1‑115 (MQ=255)
                                                                                                                                            cGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                                   >  2:1106393/1‑125 (MQ=255)
                                                                                                                                            cGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                                   <  1:1106393/125‑1 (MQ=255)
                                                                                                                                             ggCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTcc                                                          >  2:376966/1‑101 (MQ=255)
                                                                                                                                             ggCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTcc                                                          <  1:376966/101‑1 (MQ=255)
                                                                                                                                              gCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAAc        <  2:212000/150‑1 (MQ=255)
                                                                                                                                                aGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCa      <  1:724473/150‑1 (MQ=255)
                                                                                                                                                aGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCa      <  1:711741/150‑1 (MQ=255)
                                                                                                                                                    cACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGcc  <  2:289220/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCTCGGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGCACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGCCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACACAACCAGGCC  >  NC_002947/4850056‑4850353

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