Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,145,446 G→A 15.1% intergenic (‑107/+34) hemO ← / ← PP_1006 heme oxygenase/heme receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,145,4460GA15.1% 71.6 / 6.6 33intergenic (‑107/+34)hemO/PP_1006heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base G (12/16);  new base A (3/2);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGC  >  NC_002947/1145307‑1145586
                                                                                                                                           |                                                                                                                                             
gcgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCaaa                                                                                                                                     <  1:304462/150‑1 (MQ=255)
  gcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAaga                                                                                                                                   <  2:643803/150‑1 (MQ=255)
    ggggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGCCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAgg                                                                                                                                 <  2:613256/149‑1 (MQ=255)
             aGGGCGGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAg                                                                                                                        >  2:673922/1‑150 (MQ=255)
                                 gggggAGTCCTTGGAAAGAGGGGGTCTGGATACAAGACGAAACAGCGGGAGGGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCTTCTTCGCGGGCGTGCCCTCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATcc                                                                                                    <  2:149581/150‑1 (MQ=255)
                                 gggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCACCGATTAGGCCGGTACAGACGTAACGCCTCAGATATcc                                                                                                    <  1:839752/150‑1 (MQ=255)
                                     gAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAgg                                                                                                                                 <  2:720717/117‑1 (MQ=38)
                                     gAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAgg                                                                                                                                 >  1:720717/1‑117 (MQ=38)
                                           ttGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                         >  1:585150/1‑150 (MQ=255)
                                            tGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAtt                                                                                        <  2:851829/150‑1 (MQ=255)
                                              gagagggggCGTGTAAAAACGAAACGAAACAGCGGTAACGAGACAGTAAAAAAAGCCGGTATGTACGGTGCCAGTACCGGCCTCTTCGCGTGCGTGCCTGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                      <  2:799264/144‑1 (MQ=255)
                                                            aGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgcgc                                                                        <  1:673922/150‑1 (MQ=255)
                                                                    gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCAt                                                                >  2:721122/1‑150 (MQ=255)
                                                                            aGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGc                                                        <  1:721122/150‑1 (MQ=255)
                                                                             gCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCg                                                       >  2:144966/1‑150 (MQ=255)
                                                                                   gCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt                                                 <  1:144966/150‑1 (MQ=255)
                                                                                     gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGAc                                                                                                                      <  1:74621/80‑1 (MQ=255)
                                                                                     gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGAc                                                                                                                      >  2:74621/1‑80 (MQ=255)
                                                                                     gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTATTCGAGGGAGTGCCCGCAAAAAGGCCGTTACAGACGTAAATGCCCAAATATCCCAAACCAGATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGt                                               >  2:736539/1‑150 (MQ=255)
                                                                                       gaaaGTCAAAAGAGCCGGTATTTGCTTTTCCGGTACCGGCCTCGTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGttt                                             >  2:692735/4‑149 (MQ=255)
                                                                                       gaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                             <  2:577545/150‑1 (MQ=255)
                                                                                                           tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcac                                                                                            <  2:52459/84‑1 (MQ=255)
                                                                                                           tttGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcac                                                                                            >  1:52459/1‑84 (MQ=255)
                                                                                                                    tctgTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGATAACAGTCGCCTCAccata                >  2:819560/4‑147 (MQ=37)
                                                                                                                     cAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTgg               >  1:572314/1‑150 (MQ=255)
                                                                                                                      aGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGCAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGc              <  2:572314/150‑1 (MQ=255)
                                                                                                                             gCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTa       >  1:551277/1‑150 (MQ=255)
                                                                                                                             gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTc                                                    <  1:46624/105‑1 (MQ=255)
                                                                                                                             gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTc                                                    >  2:46624/1‑105 (MQ=255)
                                                                                                                               ctctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCACACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc     >  2:797911/1‑150 (MQ=255)
                                                                                                                                tctTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCt    >  2:353389/1‑150 (MQ=255)
                                                                                                                                  ttCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAATCGACGACAATTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgc  <  1:692735/150‑1 (MQ=255)
                                                                                                                                  tgcgcgGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgc  <  2:378370/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGC  >  NC_002947/1145307‑1145586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: