Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 499,014 C→T R320R (CGC→CGT
*320* (TAG→TAA
hemH →
aes ←
ferrochelatase
acetyl esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913499,0140CT97.8% 137.7 / ‑4.1 46R320R (CGC→CGT
*320* (TAG→TAA
hemH
aes
ferrochelatase
acetyl esterase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/19);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

ATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCC  >  NC_000913/498892‑499146
                                                                                                                          |                                                                                                                                    
aTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTaa                                                                                                                      >  1:80085/1‑139 (MQ=255)
                   aGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTcgcg                                                                                                   <  2:80085/139‑1 (MQ=255)
                        ttCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                            <  2:1312021/109‑1 (MQ=255)
                        ttCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                            >  1:1312021/1‑109 (MQ=255)
                                                        tatTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaataa                                                              >  2:1422961/1‑139 (MQ=255)
                                                               gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAg                                                                                                                                 <  1:361923/65‑1 (MQ=255)
                                                               gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAg                                                                                                                                 >  2:361923/1‑65 (MQ=255)
                                                               gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                            <  1:758867/70‑1 (MQ=255)
                                                               gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                            >  2:758867/1‑70 (MQ=255)
                                                                             gCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTTAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTg                                         <  1:1099910/139‑1 (MQ=255)
                                                                             gCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTa                                                                                                                       >  2:419108/1‑61 (MQ=255)
                                                                             gCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTa                                                                                                                       <  1:419108/61‑1 (MQ=255)
                                                                             gCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTg                                         <  1:844425/139‑1 (MQ=255)
                                                                                           aaaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                           <  2:496208/91‑1 (MQ=255)
                                                                                           aaaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                           >  1:496208/1‑91 (MQ=255)
                                                                                           aaaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTagag                              >  1:1205854/1‑136 (MQ=255)
                                                                                           aaaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTagag                              <  2:1205854/136‑1 (MQ=255)
                                                                                            aaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                           >  1:1337121/1‑90 (MQ=255)
                                                                                            aaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                           <  2:1337121/90‑1 (MQ=255)
                                                                                            aaTGATGGCTAATCTGGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCCGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTg                          >  2:84053/1‑139 (MQ=255)
                                                                                              tgatgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                            <  2:857905/39‑1 (MQ=255)
                                                                                              tgatgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                            >  1:857905/1‑39 (MQ=255)
                                                                                              tgatgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGaa                        >  1:837663/1‑139 (MQ=255)
                                                                                                atgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACt                      >  1:1141211/1‑139 (MQ=255)
                                                                                                atgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACt                      >  1:341326/1‑139 (MQ=255)
                                                                                                 tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTATAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATTATGCaaa                                                        <  2:1006271/104‑1 (MQ=255)
                                                                                                 tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGaa                                                                                                                  >  1:881222/1‑46 (MQ=255)
                                                                                                 tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGaa                                                                                                                  <  2:881222/46‑1 (MQ=255)
                                                                                                 tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCaaa                                                        >  1:1006271/1‑104 (MQ=255)
                                                                                                  ggCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTca                    >  1:1536567/1‑139 (MQ=255)
                                                                                                   gCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAg                                                                                                                    >  2:148118/1‑42 (MQ=255)
                                                                                                   gCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAg                                                                                                                    <  1:148118/42‑1 (MQ=255)
                                                                                                       aTCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGc                                        <  1:746477/114‑1 (MQ=255)
                                                                                                       aTCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGc                                        >  2:746477/1‑114 (MQ=255)
                                                                                                        tCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGgc              >  2:1393655/1‑139 (MQ=255)
                                                                                                            gttgCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGgctgct          >  1:1161479/1‑139 (MQ=255)
                                                                                                               gCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                                                                                        >  2:1475277/1‑58 (MQ=255)
                                                                                                               gCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                                                                                        <  1:1475277/58‑1 (MQ=255)
                                                                                                               gCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGAt       >  1:544572/1‑139 (MQ=255)
                                                                                                                  gcgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTca                    >  1:1145712/1‑123 (MQ=255)
                                                                                                                  gcgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTca                    <  2:1145712/123‑1 (MQ=255)
                                                                                                                    gTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGcc                                       <  1:968492/102‑1 (MQ=255)
                                                                                                                    gTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGcc                                       >  2:968492/1‑102 (MQ=255)
                                                                                                                    gTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGcgcc  >  1:327890/1‑139 (MQ=255)
                                                                                                                     tATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAGAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgcg    >  2:41996/1‑136 (MQ=255)
                                                                                                                     tATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgcg    <  1:41996/136‑1 (MQ=255)
                                                                                                                          |                                                                                                                                    
ATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCC  >  NC_000913/498892‑499146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: