Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 493,309 (C)5→4 intergenic (+34/‑202) adk → / → hemH adenylate kinase/ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273493,3050C.97.7% 187.4 / ‑1.9 44intergenic (+30/‑206)adk/hemHadenylate kinase/ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (14/29);  total (15/29)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACC  >  NZ_CP009273/493168‑493436
                                                                                                                                         |                                                                                                                                   
ggCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGGAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGc                                                                                                                                     <  1:543016/138‑1 (MQ=255)
ggCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGcc                                                                                                                                    <  2:569458/139‑1 (MQ=255)
ggCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGc                                                                                                                                     <  2:295752/138‑1 (MQ=255)
ggCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGc                                                                                                                                     <  1:999818/138‑1 (MQ=255)
  ctactaCTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGcccc                                                                                                                                  <  2:1277583/139‑1 (MQ=255)
  ctactaCTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGcccc                                                                                                                                  >  1:980681/1‑139 (MQ=255)
           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  2:693667/137‑1 (MQ=255)
           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  2:336276/137‑1 (MQ=255)
           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  2:633241/137‑1 (MQ=255)
           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  1:1277058/137‑1 (MQ=255)
           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  1:974637/137‑1 (MQ=255)
           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          >  1:56229/1‑137 (MQ=255)
           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  1:755667/137‑1 (MQ=255)
           cAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  1:627833/137‑1 (MQ=255)
                  gcagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          >  1:1119705/1‑130 (MQ=255)
                  gcagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  2:1119705/130‑1 (MQ=255)
                   cagaagcGGGTAATACCGAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTa                                                                                                                  >  1:468676/1‑137 (MQ=255)
                   cagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTa                                                                                                                  <  2:468676/137‑1 (MQ=255)
                   cagaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAgg                                                                                                                <  1:1264202/139‑1 (MQ=255)
                     gaagcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGt                                                                                                              <  1:4021/139‑1 (MQ=255)
                        gcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAg                                                                                                           <  2:828198/139‑1 (MQ=255)
                        gcGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAg                                                                                                           <  2:18669/139‑1 (MQ=255)
                             tAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                      <  1:699624/139‑1 (MQ=255)
                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          <  2:509006/111‑1 (MQ=255)
                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGga                                                                                                                          >  1:509006/1‑111 (MQ=255)
                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAgg                                                                                                                <  1:597028/121‑1 (MQ=255)
                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAgg                                                                                                                >  2:597028/1‑121 (MQ=255)
                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCg                                                                                              <  2:980681/139‑1 (MQ=255)
                                     aaTACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAAGCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGAAGAGGGTTAGGGTGAGGGGAACAGGCCCg                                                                                              <  2:326475/139‑1 (MQ=255)
                                                   gACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACtt                                                                                 <  2:1205364/138‑1 (MQ=255)
                                                       gCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATGCGAAAGCGCGCACGGACAGCCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAg                                                                            <  1:1272814/139‑1 (MQ=255)
                                                         aCCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGc                                                                           >  1:51066/1‑138 (MQ=255)
                                                          ccAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCaa                                                                         >  2:921043/1‑139 (MQ=255)
                                                                 ggTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAgg                                                                  <  2:237673/139‑1 (MQ=255)
                                                                            ttCGCTCTGATCTTGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGAGAGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATattca                                                       >  2:1287555/1‑139 (MQ=38)
                                                                              cgcgCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCATCGCCCCTCGGGGAGAGGGCTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATattcatt                                                     >  2:961555/1‑139 (MQ=255)
                                                                                 gcTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAg                                                                                                           >  2:1189450/1‑82 (MQ=255)
                                                                                 gcTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAg                                                                                                           <  1:1189450/82‑1 (MQ=255)
                                                                                               aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGgaga                                                                                                                        <  2:1242467/55‑1 (MQ=255)
                                                                                               aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGgaga                                                                                                                        >  1:1242467/1‑55 (MQ=255)
                                                                                               aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCa                                                                                      >  2:71312/1‑89 (MQ=255)
                                                                                               aaTCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCa                                                                                      <  1:71312/89‑1 (MQ=255)
                                                                                                        cTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTg                           <  2:490532/139‑1 (MQ=255)
                                                                                                                 aaGCGCGCACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTcc                  <  2:794740/139‑1 (MQ=255)
                                                                                                                      cgcACGGGCAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttttt             <  1:1097042/139‑1 (MQ=255)
                                                                                                                      cgcACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCtttt              >  2:361156/1‑139 (MQ=255)
                                                                                                                      cgcACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttttt             >  1:426776/1‑139 (MQ=255)
                                                                                                                      cgcACGGACAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCttttt             <  1:563110/139‑1 (MQ=255)
                                                                                                                              cAGTCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTg     <  2:56229/139‑1 (MQ=255)
                                                                                                                                 tCCCCTCGCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGAcc  >  2:1218947/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                   
GGCTACTACTCCAAAGAAGCAGAAGCGGGTAATACCAAATACGCGAAAGTTGACGGCACCAAGCCGGTTGCTGAAGTTCGCGCTGATCTGGAAAAAATCCTCGGCTAATTCGAAAGCGCGCACGGACAGTCCCCTCGCCCCCTCGGGGAGAGGGTTAGGGTGAGGGGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACC  >  NZ_CP009273/493168‑493436

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