Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,155 C→T 15.5% intergenic (+109/‑204) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1550CT15.5% 108.4 / 12.2 52intergenic (+109/‑204)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base T (3/5);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

CGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955013‑955294
                                                                                                                                               |                                                                                                                                           
cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACACCAACTACAACGTTccaggag                                                                                                                                        >  4:33797/1‑143 (MQ=255)
     cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                    <  4:262873/147‑1 (MQ=38)
     cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                  <  5:251289/149‑1 (MQ=255)
      ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                  <  8:86430/149‑1 (MQ=37)
             cTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAccc                                                                                                                          <  8:9436/149‑1 (MQ=255)
              tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCAACTAGAAGGTTCCGGTAGGCGCGGGTTTAGCcg                                                                                                                         >  7:168521/1‑149 (MQ=255)
                    ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACACCGCCTCCAAGGGTGCGGGAGGAGCGGGGTTAACCGGGaa                                                                                                                     >  6:224986/1‑147 (MQ=255)
                         gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                         >  2:111348/1‑122 (MQ=255)
                         gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                      >  4:260971/1‑141 (MQ=255)
                         gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                      <  3:260971/141‑1 (MQ=255)
                              cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAGTTTTCCTGTGCCGGACTCTTCGCGGGTAAACCCGCACCAAACAAGTTTCTGGATGACCCGGTTTACCCCCCGATAGGCCAGcac                                                                                                          >  7:203660/1‑148 (MQ=255)
                               tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          <  1:135871/131‑1 (MQ=25)
                               tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          >  2:135871/1‑131 (MQ=25)
                                           cTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTACTACAAGGTTGCTGTAGGAGCGGGTTTTCCCGCGGAGAGGGCCGCACAGaaaaaaaaaaa                                                                                             >  8:244119/1‑148 (MQ=37)
                                                tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          <  5:158257/114‑1 (MQ=18)
                                                tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          >  6:158257/1‑114 (MQ=18)
                                                tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                            >  7:203768/1‑96 (MQ=255)
                                                tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                            <  8:203768/96‑1 (MQ=255)
                                                acACGAAACCCCGGCCCCCGGGGCCTGGGTTTTTTTTTGGGGGTTTTTTACCTGTGCCGGCCTCTTGGCGGGAAAACCCGCTCCTACACGGTTCCTGTAGTTGGGGGTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCt                                                                                        <  8:168521/147‑1 (MQ=25)
                                                    gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACACGCAACTACAAGGTTCCTGTAGGGACGGGTTTAACCGCGAAGAGGCCAGCCCACAAAAAAAAAAATAACTGAatt                                                                                   >  2:217971/1‑149 (MQ=37)
                                                    gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCGACAAGGTTGCTGTTGGAGCGGGGTTTGCCGCGAAGCGGCCCGCACAGAAAAATAAAAATAACTGAttt                                                                                   >  1:17962/1‑149 (MQ=38)
                                                            cGGCCACTGGGCCGGGGTTTTTTTTTTGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCATACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                           <  3:33797/149‑1 (MQ=255)
                                                                 aCTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCGCAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt                                                                      <  5:295750/149‑1 (MQ=255)
                                                                   tgtgCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTAGCGGGTAAAACCGCTACTACAAGGTTTCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAAAGAAAAATAAAAATACATGATTTTTAAAAAATAGAtat                                                                    >  5:7526/1‑147 (MQ=37)
                                                                         cggggTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCAACAAGGTTACTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCAGAATTTAAAAAAATAGAAATTAAAAg                                                              >  6:147668/2‑149 (MQ=25)
                                                                            gTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                           <  7:244119/149‑1 (MQ=255)
                                                                            gTTTTTTTTTGGGGCTATTTTTCCTGTCCCGGCGTTTTCGCGGGTAACCCCGCTCGTACAAGGTGCCTGTTGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                           <  2:214514/149‑1 (MQ=37)
                                                                              tttttttttGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                 >  8:52506/1‑78 (MQ=255)
                                                                              tttttttttGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                 <  7:52506/78‑1 (MQ=255)
                                                                                 ttttttGGGGGTATTTTTCCTGGCCCGCCCTCTTCGCGGGTAACCCCGCTCCAACACGGTTCCTTTCGTCGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                        <  4:137332/148‑1 (MQ=255)
                                                                                 ttttttGGGGGCATTTTTCCTGTGCCGGCCTCTTCGCGGGAAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                        <  8:6879/148‑1 (MQ=255)
                                                                                   ttttGGTGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCATACAATTTTCCTGTAGTAGCGTGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaata                                                      <  5:147668/148‑1 (MQ=37)
                                                                                      tGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACAAGGTTCCTTTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaaa                                                   <  2:17962/148‑1 (MQ=255)
                                                                                          ggTATTTTTCCTGTGCCGGCCTCTTAGCGGGTTAAACCCCACCAAAAAAGTTGCAGTAGGAGCGGGGTTAACCCCGAAGAGGCCAGCACAGAAAAATAAAAAAAACAGATTTTAAAAAAAAATATCTGAAAAATTATAATAAAACAgg                                               >  4:60584/1‑148 (MQ=255)
                                                                                                  tcatgtcCCGGCCTCTTCTCGTGTAAACCCGCTCCTACAAGGTTACTTTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCAAAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                       <  5:224986/145‑1 (MQ=37)
                                                                                                          cGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTTTAGTAGCGTGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAgg                              <  3:60584/149‑1 (MQ=255)
                                                                                                                 ttCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTg                       >  7:273539/1‑149 (MQ=25)
                                                                                                                    gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                    >  4:92942/1‑149 (MQ=255)
                                                                                                                    gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                    >  4:78032/1‑149 (MQ=255)
                                                                                                                       ggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc                                 >  7:33450/1‑133 (MQ=255)
                                                                                                                       ggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc                                 <  8:33450/133‑1 (MQ=255)
                                                                                                                          aaaCCCGCTACTACAAGGTTCATGATGGATCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGTCa               <  1:217971/148‑1 (MQ=255)
                                                                                                                            aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa            >  3:126001/1‑149 (MQ=255)
                                                                                                                              ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc          >  4:62251/1‑149 (MQ=255)
                                                                                                                                gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt        <  8:273539/149‑1 (MQ=38)
                                                                                                                                gCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt        >  8:255255/1‑149 (MQ=255)
                                                                                                                                    cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAg                               <  6:278964/122‑1 (MQ=255)
                                                                                                                                    cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAg                               >  5:278964/1‑122 (MQ=255)
                                                                                                                                     tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  3:15425/1‑149 (MQ=255)
                                                                                                                                     tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  1:317371/1‑149 (MQ=255)
                                                                                                                                       cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt  <  3:78032/148‑1 (MQ=255)
                                                                                                                                             tCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCAt                           <  6:180797/117‑1 (MQ=255)
                                                                                                                                             tCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCAt                           >  5:180797/1‑117 (MQ=255)
                                                                                                                                               |                                                                                                                                           
CGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955013‑955294

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