Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,042,407 T→G K131N (AAA→AAC h16_A0947 ← Coproporphyrinogen III oxidase, Fe‑S oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,042,4070TG100.0% 51.0 / NA 14K131N (AAA→AACh16_A0947Coproporphyrinogen III oxidase, Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (7/7);  total (7/7)

CTTCGGCACCGCCGTGGATGCGCCCGAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGTTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTGCCGCCGCCGATAAAGACCGTATGCACC  >  AM260479/1042277‑1042536
                                                                                                                                  |                                                                                                                                 
cTTCGGCACCGCCGTGGATGCGCCCGAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTg                                                                                                                                <  1:625461/134‑1 (MQ=255)
     gCACCGCCGTGGATGCGCCCGAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTc                                                                                                                           >  1:434116/1‑134 (MQ=255)
               ggATGCGCCCGAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGAtc                                                                                                                 >  2:995554/1‑134 (MQ=255)
                          aGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCgg                                                                                                       >  1:57557/1‑133 (MQ=255)
                                  gAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGc                                                                                               <  2:800679/133‑1 (MQ=255)
                                                           aTGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGAcag                                                                     <  2:434116/134‑1 (MQ=255)
                                                                      cAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCa                                                          >  1:595320/1‑134 (MQ=255)
                                                                                aTGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGAcag                                                                     >  1:606715/1‑113 (MQ=255)
                                                                                aTGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGAcag                                                                     <  2:606715/113‑1 (MQ=255)
                                                                                                  cTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGccgccg                                              >  1:364980/1‑118 (MQ=255)
                                                                                                  cTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGccgccg                                              <  2:364980/118‑1 (MQ=255)
                                                                                                                         ggTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCagcagc                                                                                      <  1:420502/55‑1 (MQ=255)
                                                                                                                         ggTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCagcagc                                                                                      >  2:420502/1‑55 (MQ=255)
                                                                                                                              cGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTGCCGCCGCCGATAAAGACCGTATGCAcc  <  1:398823/134‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                 
CTTCGGCACCGCCGTGGATGCGCCCGAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGTTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTGCCGCCGCCGATAAAGACCGTATGCACC  >  AM260479/1042277‑1042536

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