Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,151 G→A 15.8% intergenic (+105/‑208) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1510GA15.8% 107.0 / 10.9 52intergenic (+105/‑208)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (23/19);  new base A (3/5);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

ACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955003‑955294
                                                                                                                                                     |                                                                                                                                               
aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAgg                                                                                                                                                 >  7:9436/1‑149 (MQ=255)
aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGCATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAgg                                                                                                                                                 >  7:6879/1‑149 (MQ=255)
   gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                              <  3:57649/149‑1 (MQ=255)
          cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACACCAACTACAACGTTccaggag                                                                                                                                        >  4:33797/1‑143 (MQ=255)
               cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                    <  4:262873/147‑1 (MQ=38)
               cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                  <  5:251289/149‑1 (MQ=255)
                ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                  <  8:86430/149‑1 (MQ=37)
                       cTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAccc                                                                                                                          <  8:9436/149‑1 (MQ=255)
                        tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCAACTAGAAGGTTCCGGTAGGCGCGGGTTTAGCcg                                                                                                                         >  7:168521/1‑149 (MQ=255)
                              ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACACCGCCTCCAAGGGTGCGGGAGGAGCGGGGTTAACCGGGaa                                                                                                                     >  6:224986/1‑147 (MQ=255)
                                   gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                         >  2:111348/1‑122 (MQ=255)
                                   gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                      <  3:260971/141‑1 (MQ=255)
                                   gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                      >  4:260971/1‑141 (MQ=255)
                                        cTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTAGTTTTCCTGTGCCGGACTCTTCGCGGGTAAACCCGCACCAAACAAGTTTCTGGATGACCCGGTTTACCCCCCGATAGGCCAGcac                                                                                                          >  7:203660/1‑148 (MQ=255)
                                         tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          >  2:135871/1‑131 (MQ=25)
                                         tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          <  1:135871/131‑1 (MQ=25)
                                                     cTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTACTACAAGGTTGCTGTAGGAGCGGGTTTTCCCGCGGAGAGGGCCGCACAGaaaaaaaaaaa                                                                                             >  8:244119/1‑148 (MQ=37)
                                                          tCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          >  6:158257/1‑114 (MQ=18)
                                                          tCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                          <  5:158257/114‑1 (MQ=18)
                                                          tCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                            <  8:203768/96‑1 (MQ=255)
                                                          tCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                            >  7:203768/1‑96 (MQ=255)
                                                          acACGAAACCCCGGCCCCCGGGGCCTGGG‑TTTTTTTTTGGGGGTTTTTTACCTGTGCCGGCCTCTTGGCGGGAAAACCCGCTCCTACACGGTTCCTGTAGTTGGGGGTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCt                                                                                        <  8:168521/147‑1 (MQ=25)
                                                              gAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACACGCAACTACAAGGTTCCTGTAGGGACGGGTTTAACCGCGAAGAGGCCAGCCCACAAAAAAAAAAATAACTGAatt                                                                                   >  2:217971/1‑149 (MQ=37)
                                                              gAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCGACAAGGTTGCTGTTGGAGCGGGGTTTGCCGCGAAGCGGCCCGCACAGAAAAATAAAAATAACTGAttt                                                                                   >  1:17962/1‑149 (MQ=38)
                                                                      cGGCCACTGGGCCGGGG‑TTTTTTTTTTGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCATACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                           <  3:33797/149‑1 (MQ=255)
                                                                           aCTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCGCAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt                                                                      <  5:295750/149‑1 (MQ=255)
                                                                             tgtgCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTAGCGGGTAAAACCGCTACTACAAGGTTTCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAAAGAAAAATAAAAATACATGATTTTTAAAAAATAGAtat                                                                    >  5:7526/1‑147 (MQ=37)
                                                                                   cggg‑gTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCAACAAGGTTACTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCAGAATTTAAAAAAATAGAAATTAAAAg                                                              >  6:147668/2‑149 (MQ=25)
                                                                                      g‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                           <  7:244119/149‑1 (MQ=255)
                                                                                      g‑TTTTTTTTTGGGGCTATTTTTCCTGTCCCGGCGTTTTCGCGGGTAACCCCGCTCGTACAAGGTGCCTGTTGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                           <  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:78032/1‑149 (MQ=255)
                                                                                                                              gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                    >  4:92942/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   >  1:317371/1‑149 (MQ=255)
                                                                                                                                               tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  3:15425/1‑149 (MQ=255)
                                                                                                                                                 cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt  <  3:78032/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955003‑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: ‑