Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,145,433 C→T 14.2% 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.2% 116.7 / 9.1 49intergenic (‑100/+41)hemO/PP_1006heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base T (4/3);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GGTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCT  >  NC_002947/1145286‑1145578
                                                                                                                                                   |                                                                                                                                                  
ggTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTc                                                                                                                                                    <  3:38515/148‑1 (MQ=255)
ggTGATCTGGATCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGCGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTc                                                                                                                                                    <  7:179289/148‑1 (MQ=255)
         gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGc                                                                                                                                          >  6:121169/1‑149 (MQ=255)
         gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGc                                                                                                                                          >  6:52946/1‑149 (MQ=255)
         gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGGCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGc                                                                                                                                          <  8:247867/149‑1 (MQ=255)
               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                                                                                       <  5:344714/148‑1 (MQ=255)
                                                            ggAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACAt                                                                                       <  4:35728/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                                                          <  2:322233/106‑1 (MQ=255)
                                                                                                                                   ccAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTc                                                          >  1:322233/1‑106 (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)
                                                                                                                                                   |                                                                                                                                                  
GGTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGCGTGCCCGCAAAGAGGCCGGTACAGACGTAACG‑CCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCT  >  NC_002947/1145286‑1145578

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: ‑