Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,121,743 A→G 8.5% C155R (TGC→CGC)  PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,121,7430AG 8.5% 128.2 / 3.1 47C155R (TGC→CGC) PP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base G (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

TTGTTGACACTGAGCAGTGCCTCGGCGCGGCCCTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTAC  >  NC_002947/6121595‑6121884
                                                                                                                                                    |                                                                                                                                             
ttgttgACACTGAGCAGTGCCTCGGCGCGGCCCTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGCGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAg                                                                                                                                              >  2:828940/1‑150 (MQ=255)
    tgACACTGAGCAGTGCCTCGGCGCGGCCCTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCgg                                                                                                                                          <  1:1033612/150‑1 (MQ=255)
        acTGTGCTGTTCCTCGGCGCGGTCCCTGAGAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGtata                                                                                                                                      <  1:786233/150‑1 (MQ=255)
          tGAGCAGTGCCTCGGCGCGGCCCTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGa                                                                                                                                    >  1:635259/1‑150 (MQ=255)
             gCAGTGCCTCGGCGCGGCCCTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATTCGTTCGACTTTGCAGCCGGTATAGATAc                                                                                                                                 <  2:494928/150‑1 (MQ=255)
                  gCCTCGGCGCGGCCCTCGATAGGCAACCGCACGGGAGAAGCACCCCTCGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATg                                                                                                                            <  1:674074/150‑1 (MQ=255)
                   ccTCGGCGCGGCCCTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGc                                                                                                                           <  2:935653/150‑1 (MQ=255)
                            ggCCCTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCctct                                                                                                                  <  2:719451/150‑1 (MQ=255)
                              cccTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGt                                                                                                                >  2:847629/1‑150 (MQ=255)
                               ccTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTgat                                                                                                               >  2:416070/1‑148 (MQ=255)
                                    aTAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATc                                                                                                          >  1:181547/1‑150 (MQ=255)
                                          aaCCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAgc                                                                                                    >  2:988523/1‑150 (MQ=255)
                                             cGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAgcggc                                                                                                 <  2:996036/150‑1 (MQ=255)
                                              gCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATAc                                                                                                                                 >  1:15889/1‑117 (MQ=255)
                                              gCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATAc                                                                                                                                 <  2:15889/117‑1 (MQ=255)
                                                   ggAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCg                                                                                           >  1:595448/1‑150 (MQ=255)
                                                      gaAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTg                                                                                        >  1:138432/1‑150 (MQ=255)
                                                      gaAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTg                                                                                        <  1:459774/150‑1 (MQ=255)
                                                          caccacTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCa                                                                                    >  1:929759/1‑150 (MQ=255)
                                                            ccacTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGc                                                                                                                      >  2:723323/1‑114 (MQ=255)
                                                            ccacTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGc                                                                                                                      <  1:723323/114‑1 (MQ=255)
                                                                    tAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTgcggc                                                                          >  1:403039/1‑150 (MQ=255)
                                                                       cacCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGt                                                                       <  2:507825/150‑1 (MQ=255)
                                                                         ccAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTgc                                                                     <  2:897818/150‑1 (MQ=255)
                                                                                      gTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGTCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAg                                                        <  2:403039/150‑1 (MQ=255)
                                                                                      gTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAg                                                        <  2:138432/150‑1 (MQ=255)
                                                                                       tCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAgg                                                       <  1:847629/150‑1 (MQ=255)
                                                                                           tttCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTgc                                                                     >  1:726542/1‑132 (MQ=255)
                                                                                           tttCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTgc                                                                     <  2:726542/132‑1 (MQ=255)
                                                                                             tcgtcgAGGCTGGTGTACAAGCGCTGCGCCTGGGCATGGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCTCTGAGCCAGGGCGATg                                                 <  2:404912/150‑1 (MQ=255)
                                                                                                            tACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTg                                  <  1:926301/150‑1 (MQ=255)
                                                                                                              cAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGa                                >  1:202104/1‑150 (MQ=255)
                                                                                                                aGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAAt                              >  2:1009198/1‑150 (MQ=255)
                                                                                                                 gcgcTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATa                             >  2:109102/1‑150 (MQ=255)
                                                                                                                   gcTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCCGCGCCTTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATaaa                           >  1:73013/1‑150 (MQ=255)
                                                                                                                   gcTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATaaa                           <  2:929759/150‑1 (MQ=255)
                                                                                                                         ccTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGc                     <  1:828940/150‑1 (MQ=255)
                                                                                                                                   gATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTg           >  1:404807/1‑150 (MQ=255)
                                                                                                                                   gATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAACAAAGATGGCCCGACTCCTg           >  1:951784/1‑150 (MQ=255)
                                                                                                                                    aTGCGTTCGACTTTGCGGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTg                                  <  1:144618/126‑1 (MQ=255)
                                                                                                                                    aTGCGTTCGACTTTGCGGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTg                                  <  1:145355/126‑1 (MQ=255)
                                                                                                                                    aTGCGTTCGACTTTGCGGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTg                                  >  2:145355/1‑126 (MQ=255)
                                                                                                                                    aTGCGTTCGACTTTGCGGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTg                                  >  2:144618/1‑126 (MQ=255)
                                                                                                                                    aTGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTgg          >  1:778321/1‑150 (MQ=255)
                                                                                                                                         ttCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAg     >  2:195734/1‑150 (MQ=255)
                                                                                                                                            gACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGTATGCGGCGGTGCGCCCCTGAGCCAGGTCGATGGATATCGTGGGCTTGGAATAAAGATGGCCCGATTCCTGGGTCAGTAc  >  1:446450/1‑150 (MQ=255)
                                                                                                                                            gACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTAc  >  1:905740/1‑150 (MQ=255)
                                                                                                                                                    aGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTgg                                 <  2:591397/111‑1 (MQ=255)
                                                                                                                                                    aGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTgg                                 >  1:591397/1‑111 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TTGTTGACACTGAGCAGTGCCTCGGCGCGGCCCTCGATAGGCAACCGCACGGGAGAAGCACCACTGGCTAACACCAGTTGGCCGTAGTCCATTTCGCCGAGGCTGGTGTACAAGCGCTGCGCCTGGGCATCGATGCGTTCGACTTTGCAGCCGGTATAGATACGAATGCCAAGCCTCTGTTCGATCGCCAGCGGCGTTTCGGTGGCCAGGGCTGCGGCGGTGCGCCCCTGAGCCAGGGCGATGGATAGCGCGGGCTTGGAATAAAGATGGCCCGACTCCTGGGTCAGTAC  >  NC_002947/6121595‑6121884

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