Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,293,106 (G)7→8 coding (1595/1659 nt) nrfE → putative cytochrome c‑type biogenesis protein NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,293,0991.G96.8% 102.7 / ‑3.2 31coding (1588/1659 nt)nrfEputative cytochrome c‑type biogenesis protein NrfE
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (16/14);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACC  >  NC_000913/4292965‑4293231
                                                                                                                                        |                                                                                                                                       
cACGCCGTCAACAAATGATGGATCCGTCAATACGCTGGAACGGCATCCCTGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTggggg                                                                                                                                   <  2:317334/140‑1 (MQ=255)
  cGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTCTGTATGTACAACGCGGTGTGCGCTGGATCTggggggg                                                                                                                                 >  1:321662/1‑140 (MQ=255)
     cGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTggggggggag                                                                                                                              <  1:775452/140‑4 (MQ=255)
                 atgGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTtgatga                                                                                                                  >  1:663103/1‑140 (MQ=255)
                     aaCCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGc                                                                                                              >  1:18400/1‑140 (MQ=255)
                     aaCCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGc                                                                                                              >  1:691664/1‑140 (MQ=255)
                      aCCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGc                                                                                                              >  1:135178/1‑139 (MQ=255)
                      aCCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGc                                                                                                              <  2:135178/139‑1 (MQ=255)
                                          gCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATg                                                                                        <  2:45472/140‑1 (MQ=255)
                                             tCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGAGGGGGGAGAATTGTTGCTGATTGCGGGCGCATTGTTAAGCGGATGGCg                                                                                     >  2:210302/1‑140 (MQ=255)
                                                   aCTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGgaa                                                                               >  1:492047/1‑140 (MQ=255)
                                                     tGGTATGCGGTCATGGGGGAGAAACTGGGCCGGGTCCGTAACGTTTTTCGTTGGTTGTACAAAGCCGGGGGGGGTGGGTCTTGGGGGGGGGAGTTGTTTAAGGTTTCGAGGCGCTTTCTTAAACGGATGCCggggaagga                                                                             >  2:388520/1‑135 (MQ=255)
                                                             ggTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGAt                                                                     <  1:597477/140‑1 (MQ=255)
                                                              gTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGAt                                                                     >  2:479814/1‑139 (MQ=255)
                                                              gTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGAt                                                                     <  1:479814/139‑1 (MQ=255)
                                                               tCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGa                                                                   >  1:793546/1‑140 (MQ=255)
                                                                  tGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGACGGGGGGGGATGAAGCGCGATGAATa                                                                >  1:393117/1‑140 (MQ=255)
                                                                  tGGGGGAGAAAACTGGGCCGGATCGTTACGCATTTTGTTCGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATa                                                                <  2:275575/140‑1 (MQ=255)
                                                                  tGGGGGAGAAAACTGGGCCGGATCGTTACACTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATa                                                                <  1:501445/140‑1 (MQ=255)
                                                                      ggAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAgg                                                            <  2:663103/140‑1 (MQ=255)
                                                                       gagaAAACTGGGCCGGAACGTTACGCTTTTCTTTTGTATCTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAggg                                                           <  2:321662/140‑1 (MQ=255)
                                                                               tGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCAc                                                   >  1:13348/1‑140 (MQ=255)
                                                                                 ggCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGc                                                 >  1:434485/1‑140 (MQ=255)
                                                                                   ccGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTa                                               <  2:793546/140‑1 (MQ=255)
                                                                                      gATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCt                                                <  2:555397/136‑1 (MQ=255)
                                                                                      gATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCt                                                >  1:555397/1‑136 (MQ=255)
                                                                                                    ttCGTTTGTATGTACTAAGCGGAGTACGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCt                               <  1:388520/140‑1 (MQ=255)
                                                                                                              tgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCa                     >  2:507841/1‑140 (MQ=255)
                                                                                                                 aCAAAGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGc                  >  1:195827/1‑140 (MQ=255)
                                                                                                                   aaaGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTc                <  2:18400/140‑1 (MQ=255)
                                                                                                                   aaaGCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTc                <  2:691664/140‑1 (MQ=255)
                                                                                                                      gCGGTGTGCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAgg             <  2:195827/140‑1 (MQ=255)
                                                                                                                          tgtgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGt         <  2:594602/140‑1 (MQ=255)
                                                                                                                          tgtgCGCTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGt         >  1:637228/1‑140 (MQ=255)
                                                                                                                               gcTGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCCCGCTACTGCTGTAATTAACCTGTTTTGCGCACGCTCAGGTCGTAGaca    >  2:281246/1‑140 (MQ=255)
                                                                                                                                 tGGATCTGGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACAcc  >  2:451123/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACC  >  NC_000913/4292965‑4293231

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