Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,152 T→A 14.5% intergenic (+106/‑207) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1520TA14.5% 107.7 / 7.6 52intergenic (+106/‑207)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base A (3/5);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955006‑955294
                                                                                                                                                   |                                                                                                                                              
gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                              <  3:57649/149‑1 (MQ=255)
       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                                                                                                                          >  6:158257/1‑114 (MQ=18)
                                                       tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          <  5:158257/114‑1 (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                                                                                                                                 <  7:52506/78‑1 (MQ=255)
                                                                                     tttttttttGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                 >  8:52506/1‑78 (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:78032/1‑149 (MQ=255)
                                                                                                                           gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                    >  4:92942/1‑149 (MQ=255)
                                                                                                                              ggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc                                 <  8:33450/133‑1 (MQ=255)
                                                                                                                              ggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc                                 >  7:33450/1‑133 (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)
                                                                                                                                                   |                                                                                                                                              
GCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955006‑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: