Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,721,455 G→C 8.2% G264G (GGC→GGG)  fnrC ← Crp/Fnr family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,721,4550GC 8.2% 143.1 / 4.4 49G264G (GGC→GGG) fnrCCrp/Fnr family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base C (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

GTTGAATGCAGCGCTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGTTTCCAGGCGTAAACCCAG  >  NC_002947/3721306‑3721604
                                                                                                                                                     |                                                                                                                                                     
gTTGAATGCAGCGCTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGCGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAgggg                                                                                                                                                       <  1:133869/150‑1 (MQ=255)
 ttGAATGCAGCGCTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCAAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGc                                                                                                                                                      >  1:361371/1‑150 (MQ=255)
   gAATGCAGCGCTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGccc                                                                                                                                                    <  2:529460/150‑1 (MQ=255)
   gAATGCAGCGCTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGccc                                                                                                                                                    >  2:503452/1‑150 (MQ=255)
    aaTGCAGCGCTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCt                                                                                                                                                   >  2:1039999/1‑150 (MQ=255)
      tGCAGCGCTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGAGGACAACGATCCCGGCGTCGTGGAGGGGGCCTTGAGTGGCCCTTg                                                                                                                                                 >  2:850264/1‑150 (MQ=255)
        cAGCGCTGGAAACAGCCGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGACCTTGGc                                                                                                                                               <  2:566037/150‑1 (MQ=255)
         ggcgcTGGAACCAGCAGATGGCCCATGGCCGGACGTGAGCGAGAGTGTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCCACATTCACGGCGTCGTGCAGGGTGCCTTCAGGGGCACTTGGcc                                                                                                                                              <  2:231944/149‑1 (MQ=255)
          gcgcTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGcc                                                                                                                                                     >  2:711220/1‑142 (MQ=255)
          gcgcTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGcc                                                                                                                                                     <  1:711220/142‑1 (MQ=255)
          gcgcTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAAAATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGcag                                                                                                                                             <  2:879986/150‑3 (MQ=255)
             cTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTcac                                                                                                                                          <  1:1070414/150‑1 (MQ=255)
              tGGAAACAGCAGATGGCCCTTTGCCGGACGTGAGCGAGAGTTTTTCCAATCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTcaca                                                                                                                                         >  1:256018/1‑150 (MQ=255)
                        agATGGCCCATTGCCGGACGTGAGCGAGAACTTTTTCTAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCa                                                                                                                               <  2:88181/150‑1 (MQ=255)
                            ggCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCtt                                                                                                                           <  1:987965/150‑1 (MQ=255)
                               ccATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGtt                                                                                                                        <  2:419992/150‑1 (MQ=255)
                                cATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTc                                                                                                                       >  1:115526/1‑150 (MQ=255)
                                  ttGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGa                                                                                                                     <  1:915094/150‑1 (MQ=255)
                                      cGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAgcgc                                                                                                                 <  1:203200/150‑1 (MQ=255)
                                              aGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTg                                                                                                         <  2:813264/150‑1 (MQ=255)
                                                 gagaGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAg                                                                                                      <  2:184924/150‑1 (MQ=255)
                                                      tttttACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCata                                                                                                 >  2:217437/1‑150 (MQ=255)
                                                             aaGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGaa                                                                                          >  1:626301/1‑150 (MQ=255)
                                                                  ttttGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCt                                                                                     >  2:67089/1‑150 (MQ=255)
                                                                    ttGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTAc                                                                                   <  2:624460/150‑1 (MQ=255)
                                                                        cGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAAc                                                                               <  1:253812/150‑1 (MQ=255)
                                                                         gggTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAAcc                                                                              >  1:1103190/1‑150 (MQ=255)
                                                                            ttGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGc                                                                           >  1:408292/1‑150 (MQ=255)
                                                                                  tGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCa                                                                     <  2:115526/150‑1 (MQ=255)
                                                                                      gggAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTc                                                                 >  1:333801/1‑150 (MQ=255)
                                                                                         aGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTcgac                                                              <  1:713035/150‑1 (MQ=255)
                                                                                           tGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTcgacga                                                            >  1:474009/1‑150 (MQ=255)
                                                                                                  cTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTTCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAg                                                     >  2:310083/1‑150 (MQ=255)
                                                                                                     acaGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAgcgc                                                  <  1:503452/150‑1 (MQ=255)
                                                                                                      caGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAgcgcg                                                 >  1:391729/1‑150 (MQ=255)
                                                                                                       aGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCaa                                                                                             <  2:468283/105‑1 (MQ=255)
                                                                                                       aGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCaa                                                                                             >  1:468283/1‑105 (MQ=255)
                                                                                                             gcgGCCAACAATAACGACGTCTAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGCGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGAc                                          <  2:256018/150‑1 (MQ=255)
                                                                                                              cgGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACg                                         <  2:474009/150‑1 (MQ=255)
                                                                                                                   aaCATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGa                                    >  1:325651/1‑150 (MQ=255)
                                                                                                                          aCGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTcga                                                               >  2:62564/1‑116 (MQ=255)
                                                                                                                          aCGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTcga                                                               <  1:62564/116‑1 (MQ=255)
                                                                                                                           cGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCgg                            >  1:893354/1‑150 (MQ=255)
                                                                                                                                        ggTGCCTTCAGGGCCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGTTTCCa               >  2:857477/1‑150 (MQ=255)
                                                                                                                                        ggTGCCTTCAGGGCCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGTTTCCa               >  2:857689/1‑150 (MQ=255)
                                                                                                                                         gTGCCTTCAGGGGCCCTCGGCCTAACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGTTTCCAg              <  2:361371/150‑1 (MQ=255)
                                                                                                                                         gTGCCTTCAGGGCCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGTTTCCAg              <  1:857689/150‑1 (MQ=255)
                                                                                                                                         gTGCCTTCAGGGCCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGTTTCCAg              <  1:857477/150‑1 (MQ=255)
                                                                                                                                            ccTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGtt                   <  1:27268/142‑1 (MQ=255)
                                                                                                                                            ccTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGtt                   >  2:27268/1‑142 (MQ=255)
                                                                                                                                                     gCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCCGCGCGCGCAGACGGGCGACGGATCGGCACACGGTTTCCCGGCGTAAACCCAg  <  2:1066618/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGAATGCAGCGCTGGAAACAGCAGATGGCCCATTGCCGGACGTGAGCGAGAGTTTTTACAAGCGTTTTGGCGGGTTGGCTTGAGGGGAGTGCTCGCCTCACAGCGCCGCGGCCAACATTCACGGCGTCGAGCAGGGTGCCTTCAGGGGCCCTTGGCCTCACATTCACCGCCACCTTGTTCGAGCGCCATGAGTGCAGGCATATCAAGAATTTCTACCAACCTGCCATGCAGGTCGACGATACCCAGCGCGCGCAGACGGGCGACGGATCGGCACACGGTTTCCAGGCGTAAACCCAG  >  NC_002947/3721306‑3721604

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