Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,145,448 A→G 17.4% intergenic (‑115/+26) 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,4480AG17.4% 95.4 / 13.8 46intergenic (‑115/+26)hemO/PP_1006heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base A (18/20);  new base G (4/4);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCAC  >  NC_002947/1145301‑1145594
                                                                                                                                                   |                                                                                                                                                   
gcgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCa                                                                                                                                                     >  6:344714/1‑148 (MQ=255)
gcgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCACTACCGGCCTCTTCGCGGGCGTGCCCGCaa                                                                                                                                                    <  3:62507/149‑1 (MQ=255)
   cTGGGAGCGCAGGGCTGGGCGGTCGGTCAGGGGGGAGCGCTTGGAAAGAGCGCGTCTAGATACAACACCAAACAGCGGGAGCGAGAGAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAaga                                                                                                                                                 <  7:146089/149‑1 (MQ=255)
      ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAagag                                                                                                                                                <  4:266598/147‑1 (MQ=255)
      ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGc                                                                                                                                              <  5:121169/149‑1 (MQ=255)
      ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGc                                                                                                                                              >  4:140625/1‑149 (MQ=255)
        aGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCg                                                                                                                                            <  5:269130/149‑1 (MQ=255)
         gcgcAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCgg                                                                                                                                           >  5:231315/1‑149 (MQ=255)
         gcgcAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCgg                                                                                                                                           >  3:35728/1‑149 (MQ=255)
          cgcAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCgg                                                                                                                                           >  6:349792/1‑148 (MQ=255)
           gcAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTa                                                                                                                                         <  7:285255/149‑1 (MQ=255)
                           ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCg                                                                                                                                            <  4:285439/130‑1 (MQ=255)
                           ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCa                                                                                                                        <  5:52946/149‑1 (MQ=255)
                           ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGACCGCAAAGAGGCCg                                                                                                                                            >  3:285439/1‑130 (MQ=255)
                                  ggggAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATcc                                                                                                                 >  2:125843/1‑149 (MQ=255)
                                             ggAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCTGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACAt                                                                                                       <  1:274018/148‑1 (MQ=255)
                                             ggAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACAt                                                                                                       <  4:35728/148‑1 (MQ=255)
                                             ggAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACAt                                                                                                       <  5:344714/148‑1 (MQ=255)
                                             ggAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACAt                                                                                                       <  6:231315/148‑1 (MQ=255)
                                                 agagGGCGTCTATATACAAGACGAAACAGCCGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATa                                                                                                  <  3:320673/149‑1 (MQ=255)
                                                       cGTCTAGATACAAGACGAAACAGCGGCAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAg                                                                                            <  8:174138/149‑1 (MQ=255)
                                                             gATACAAGACGAAACAGCGCGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTgcgc                                                                                      >  7:242185/1‑149 (MQ=255)
                                                                    gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCa                                                                               >  3:318271/1‑149 (MQ=255)
                                                                    caCGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCCCGGGCGTGCCCGCAAACAGGCCCGTACAGGCGAAACG‑CATCAGATATCACCAACCACAGTGATACCGAAGTTGCGCCCCGGCa                                                                               >  4:153487/2‑149 (MQ=255)
                                                                      cGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                                             <  3:140625/149‑1 (MQ=255)
                                                                        aaaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGt                                                                           >  6:364772/1‑149 (MQ=255)
                                                                              cGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCg                                                                     >  6:17857/1‑149 (MQ=255)
                                                                                     gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATg                                                              >  3:92112/1‑149 (MQ=255)
                                                                                      agaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGt                                                             <  5:17857/149‑1 (MQ=255)
                                                                                         cAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtgg                                                          <  1:125843/149‑1 (MQ=255)
                                                                                            tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGa                                                                <  4:114134/140‑1 (MQ=255)
                                                                                            tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGa                                                                >  3:114134/1‑140 (MQ=255)
                                                                                                 aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                                  <  4:92112/149‑1 (MQ=255)
                                                                                                       ggTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGGGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcaca                                            <  5:349792/149‑1 (MQ=255)
                                                                                                                gTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCa                                   >  2:324106/1‑149 (MQ=255)
                                                                                                                    tctgTACCGGCCTCTTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGt                                                                   <  4:115363/110‑1 (MQ=37)
                                                                                                                    ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTc                                                                          >  1:322233/1‑106 (MQ=255)
                                                                                                                    ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTc                                                                          <  2:322233/106‑1 (MQ=255)
                                                                                                                           cGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtc                        <  1:324106/149‑1 (MQ=255)
                                                                                                                            ggCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACAACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGAGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtct                       >  6:359913/1‑148 (MQ=255)
                                                                                                                             gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                                  >  5:339024/1‑121 (MQ=255)
                                                                                                                             gCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                                  <  6:339024/121‑1 (MQ=255)
                                                                                                                               ctctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg                    >  1:119109/1‑149 (MQ=255)
                                                                                                                                tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg                    >  1:68573/1‑149 (MQ=255)
                                                                                                                                tctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc                   >  7:47799/1‑149 (MQ=255)
                                                                                                                                 ctTTGCGGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCt                  >  4:329220/1‑149 (MQ=255)
                                                                                                                                      cgGGCACGCCCGCGAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTgcgc              <  3:329220/148‑1 (MQ=255)
                                                                                                                                                 gCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCAc  >  4:66680/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCCCAC  >  NC_002947/1145301‑1145594

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