Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,151 G→A 16.6% 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,1510GA16.6% 90.9 / 10.8 42intergenic (+105/‑208)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (21/14);  new base A (3/4);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

CAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955008‑955294
                                                                                                                                                |                                                                                                                                               
cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                            >  1:269519/1‑149 (MQ=255)
cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                            <  2:269519/149‑1 (MQ=255)
         gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                  <  3:65512/149‑1 (MQ=25)
          cgtaTGTGCTGATGACCAACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                  <  8:177015/146‑1 (MQ=255)
          cGGATGTGCTGATGACCTACGCAAGCTGAAACTGCTTCCTGTTTAACGAAAACCACGGCCACTGTGCCGGGGTTTTTTTTTGGTGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                  <  6:193968/149‑1 (MQ=25)
            gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt                                                                                                                                <  1:355330/149‑1 (MQ=255)
                                acgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                         <  1:29133/120‑1 (MQ=255)
                                acgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                         >  2:29133/1‑120 (MQ=255)
                                                   tttCACGAAAACCCCGGCCACTGTGCCGGGGTTTTATTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTGCTACAAGGTTGCGGTAGGAGCGGGGTTACCCGCGGAGTGGCCAGCACAGAAAAATAAAAATTcc                                                                                         >  1:59425/1‑149 (MQ=37)
                                                         gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                            <  1:246120/92‑1 (MQ=255)
                                                         gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                            >  2:246120/1‑92 (MQ=255)
                                                            aaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCACTTAGCGGGTAAACCCGCTACTACAAGGTTGCGGTAGGAGCGGGTTCACCCGGGAAGAGGGCACCACAGAAAAATAAAAATAACTGAATTTc                                                                                 >  8:152452/1‑148 (MQ=25)
                                                                   gCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGGTACTACAAGGTTCCTGTAGGAGCCGGTTTAACAGCGAAAAGGCCAGCACAGAAAAATAAGAATACATGATTTTCAAAAAATa                                                                         >  5:371418/1‑149 (MQ=37)
                                                                                    ttttttttGGGGGTATTTTCCCTGTGCCGCCCTCTTCGCGGGAAAACCCGCTCCAAAAAGGTTCCTGTAGGAGGGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                                                          <  3:192179/148‑1 (MQ=38)
                                                                                               ggTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAAACCCCACCGAAAAGGGTGGGGGTTGGGCGGGGTTGACCGCGGCGAGGGCCGCCCCAAAAAAAAAAAAAACCGGAGTTTGAAAAAAAAGATATGGAAAATTGTAATAAAAAAgc                                               >  4:224610/1‑147 (MQ=255)
                                                                                                            tgCCGGCCTCTTCGCGGTTAAACCCGCTCCTACAAGGTTCCTGTAGCAGCGGGTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTa                                  <  2:59425/148‑1 (MQ=255)
                                                                                                             gCCGGCCTCTTCGCGGGTAAACCCGCTACTACAAGGTTGCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCAGATTTTCAAAAAATAAATATGAAAAGTTATAATAAAACAGGAAATTATCCTTTa                                  >  8:272486/1‑147 (MQ=37)
                                                                                                                          cgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                   >  4:330129/1‑149 (MQ=255)
                                                                                                                           gggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGCTATACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATAATCCATTACAGGCATCCTGTGTCg                  >  8:115478/1‑149 (MQ=39)
                                                                                                                             gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                 >  7:364096/1‑148 (MQ=255)
                                                                                                                             gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                 <  8:364096/148‑1 (MQ=255)
                                                                                                                               aaaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc              >  1:66024/1‑149 (MQ=255)
                                                                                                                                aaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt             <  3:330129/149‑1 (MQ=255)
                                                                                                                                aaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt             >  1:145802/1‑149 (MQ=255)
                                                                                                                                 aCCCTCTCCTACAAGGTTCCTGTAGGAGCGGTTTTACCCGCGAATAGTCCATCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa            <  7:1354/149‑1 (MQ=255)
                                                                                                                                 aCCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa                                            >  5:172596/1‑117 (MQ=25)
                                                                                                                                 aCCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa                                            <  6:172596/117‑1 (MQ=25)
                                                                                                                                  cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                   >  8:264136/1‑141 (MQ=37)
                                                                                                                                  cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                   <  7:264136/141‑1 (MQ=37)
                                                                                                                                   ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc          >  1:11859/1‑149 (MQ=38)
                                                                                                                                   ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc          >  3:75950/1‑149 (MQ=255)
                                                                                                                                      cTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGATGCCATCACATAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc       <  7:272486/149‑1 (MQ=255)
                                                                                                                                         cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttct      <  2:66024/147‑1 (MQ=255)
                                                                                                                                          tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaaa                                                   >  7:80354/1‑101 (MQ=255)
                                                                                                                                          tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaaa                                                   <  8:80354/101‑1 (MQ=255)
                                                                                                                                          tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  6:133041/1‑149 (MQ=255)
                                                                                                                                          tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  1:12392/1‑149 (MQ=255)
                                                                                                                                           tcAACGTTCCTTTATGATGGGGTTTAACCGCGAACTCGCCATCACATAAAAATAAAAATACCTTATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt  <  2:21127/148‑1 (MQ=255)
                                                                                                                                           aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt  >  6:210944/1‑149 (MQ=255)
                                                                                                                                           aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  5:136370/1‑148 (MQ=255)
                                                                                                                                             aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt  >  2:46769/1‑147 (MQ=255)
                                                                                                                                             aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt  >  2:307652/1‑147 (MQ=255)
                                                                                                                                               ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAg                                                                        >  6:47067/1‑75 (MQ=255)
                                                                                                                                               ggTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAg                                                                        <  5:47067/75‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955008‑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: