Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,145,433 C→T 14.8% intergenic (‑100/+41) 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,4330CT14.8% 129.7 / 12.5 55intergenic (‑100/+41)hemO/PP_1006heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (5/3);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-02

GGTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGC  >  NC_002947/1145286‑1145577
                                                                                                                                                   |                                                                                                                                                 
ggTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTc                                                                                                                                                   <  7:272138/148‑1 (MQ=255)
         gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTg                                                                                                                                         <  8:96832/148‑1 (MQ=255)
              ggCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGt                                                                                                                                          >  4:131043/1‑142 (MQ=255)
              ggCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGt                                                                                                                                          <  3:131043/142‑1 (MQ=255)
                cgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCaaa                                                                                                                                 <  5:77661/149‑1 (MQ=255)
                cgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCaaa                                                                                                                                 <  2:165012/149‑1 (MQ=255)
                     ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAAGAGGc                                                                                                                            <  5:120509/149‑1 (MQ=255)
                     ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAAGAGGc                                                                                                                            >  6:120509/1‑149 (MQ=255)
                     ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGc                                                                                                                            >  5:11557/1‑149 (MQ=255)
                       aGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAAGAGGCCg                                                                                                                          <  1:232476/149‑1 (MQ=255)
                            aGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGtc                                                                                                                       >  7:249351/1‑146 (MQ=255)
                                cTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTCGCGTTGCCAGTACCGGCCTCTTCGCG‑AGCGTGCCCGCAAAGAGGCCGGTACAGACg                                                                                                                 >  6:172030/1‑149 (MQ=255)
                                   ggCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTaa                                                                                                              <  6:11557/149‑1 (MQ=255)
                                   ggCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTaa                                                                                                              >  8:220795/1‑149 (MQ=255)
                                          ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAagag                                                                                                                              <  3:248329/126‑1 (MQ=255)
                                          ggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCG‑GGCACGCCCGCGAagag                                                                                                                              >  4:248329/1‑126 (MQ=255)
                                            tcATGGGGGAGTACTTGGAAAGAGGTCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCACGCAAAGAGGCCTGTACACACGTAACGCCt                                                                                                         >  6:50747/1‑145 (MQ=255)
                                                     aGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                            <  7:220795/149‑1 (MQ=255)
                                                         cTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACca                                                                                         >  3:293660/1‑148 (MQ=255)
                                                            ggAAAGTGGGCTTCTTTAAACAACACTAAACAGCGGGACCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACcaca                                                                                       <  2:173731/147‑1 (MQ=255)
                                                             gAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAtt                                                                                     >  6:4986/1‑148 (MQ=255)
                                                                agagGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                   >  4:110096/1‑147 (MQ=255)
                                                                agagGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGa                                                                                   <  3:110096/147‑1 (MQ=255)
                                                                   gggCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTCGCGTTGCCAGTACCGGCCTCTTCGCG‑AGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCg                                                                              <  5:172030/149‑1 (MQ=255)
                                                                     gCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                            <  3:269796/149‑1 (MQ=255)
                                                                               aCAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGcccc                                                                  >  6:77674/1‑149 (MQ=255)
                                                                               aCAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTAGCGTTCCCAGTACCGCCCTCTACGCG‑GGCGTGCCCGCAAAGAGGCCCGAACAGACGTAACGCCTCAGAAATCACAAACCACAATGAAACCGAAGTAGCGaccc                                                                  >  2:264143/1‑149 (MQ=255)
                                                                                cAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                            <  3:7474/138‑1 (MQ=255)
                                                                                cAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGaa                                                                            >  4:7474/1‑138 (MQ=255)
                                                                                        aaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGt                                                          >  8:223544/1‑148 (MQ=255)
                                                                                               ggAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                            <  7:15569/139‑1 (MQ=255)
                                                                                               ggAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                            >  8:15569/1‑139 (MQ=255)
                                                                                                taGCGAGAAATTAAAAAAAGAAGTTATTTTCGATGCCATTACCGGCCTATAAGCA‑GGCGTGCACGCAAAGAGGCCGGTACAGACGTAACTACTAAGATATCAAAAACAACAATGATACCGAAGTTGCGCACCGGCATGGTCAGGCGGTc                                                 >  6:68286/2‑149 (MQ=255)
                                                                                                gAGCGAGACAGTAAAAAAAGCCGGTATTTTCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGTCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTc                                                 <  5:50747/149‑1 (MQ=255)
                                                                                                  gCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGa                                               >  8:246543/1‑149 (MQ=255)
                                                                                                   cGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt                                              <  1:145985/149‑1 (MQ=255)
                                                                                                           tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                            <  1:152161/127‑1 (MQ=255)
                                                                                                           tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATg                                                            >  2:152161/1‑127 (MQ=255)
                                                                                                                aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                 <  5:77674/149‑1 (MQ=255)
                                                                                                                aaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                 <  5:4986/149‑1 (MQ=255)
                                                                                                                 aaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTggg                                >  1:34574/1‑149 (MQ=255)
                                                                                                                 aaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTggg                                <  4:293660/149‑1 (MQ=255)
                                                                                                                  aGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGt                               <  7:223544/149‑1 (MQ=255)
                                                                                                                    ccGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGTGGGGTCCCCGCAAAGGGGCCGGTACAGCGGTACCGCCGCAGATTTCCCAACCCACTTTGTTACCGAGGTTGCGCCCCGGCATGGTCAGGCGGGCGATTTTGCCTGGCTGGGTc                              >  1:54805/1‑149 (MQ=255)
                                                                                                                      ggTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcaca                           <  2:34574/149‑1 (MQ=255)
                                                                                                                                  gCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCAccc               >  7:62161/1‑149 (MQ=255)
                                                                                                                                            gCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg      >  2:169261/1‑148 (MQ=255)
                                                                                                                                              ctctTTGCG‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg   >  4:222268/1‑149 (MQ=255)
                                                                                                                                              ctctTTGCG‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg   >  8:195131/1‑149 (MQ=255)
                                                                                                                                              ctctTTGCG‑GGCACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg   >  5:159207/1‑149 (MQ=255)
                                                                                                                                              ctctTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAg   >  5:101450/1‑149 (MQ=255)
                                                                                                                                               tctTCGCG‑GGCGTGTCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc  >  6:310689/1‑149 (MQ=255)
                                                                                                                                               tctTCGCG‑GGCGTGCCCGCAAAGAGGCCGCTACAGACGTAACGCCTCAGATATCCCAACCCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc  >  3:299686/1‑149 (MQ=255)
                                                                                                                                               tctTCGCG‑GGCGTGCCAGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGACACAGCCGCCTCACACTCGCCGTCGTAGc  >  7:194241/1‑149 (MQ=255)
                                                                                                                                                 ttCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGtt                                           <  4:252193/106‑1 (MQ=255)
                                                                                                                                                 ttCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGtt                                           >  3:252193/1‑106 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCG‑GGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGC  >  NC_002947/1145286‑1145577

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