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% 77.2 / NA 21K131N (AAA→AACh16_A0947Coproporphyrinogen III oxidase, Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (11/10);  total (11/10)

CACCGCCGTGGATGCGCCCGAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGTTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTGCCGCCGCCGATAAAGACCGTATGCAC  >  AM260479/1042283‑1042535
                                                                                                                            |                                                                                                                                
cACCGCCGTGGATGCGCCCGAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTc                                                                                                                          >  1:519676/1‑133 (MQ=255)
                   gAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCgg                                                                                                      >  1:373009/1‑134 (MQ=255)
                    aGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCgg                                                                                                      >  1:828862/1‑133 (MQ=255)
                                  ccGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGgcagc                                                                                        <  1:366355/133‑1 (MQ=255)
                                      tCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTg                                                                                                                               >  2:920323/1‑90 (MQ=255)
                                      tCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTg                                                                                                                               <  1:920323/90‑1 (MQ=255)
                                      tCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAgcgc                                                                                   >  2:853021/1‑134 (MQ=255)
                                        gTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAgcgcg                                                                                  <  1:853021/133‑1 (MQ=255)
                                             aGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGata                                                                            <  2:519676/134‑1 (MQ=255)
                                                        ccGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGAcagcag                                                                 >  2:255594/1‑134 (MQ=255)
                                                                  gCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCAtg                                                       <  1:255594/134‑1 (MQ=255)
                                                                      gTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCAtgcctg                                                   <  2:1024554/134‑1 (MQ=255)
                                                                       ttGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCAtgcctgc                                                  >  2:125164/1‑134 (MQ=255)
                                                                                    cgcgGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAgg                                     >  1:60045/1‑134 (MQ=255)
                                                                                                       gTCGGCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTGccgccgccg                   <  2:528760/132‑1 (MQ=255)
                                                                                                           gCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCagcagc                                                                                     >  1:69263/1‑63 (MQ=255)
                                                                                                           gCCTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCagcagc                                                                                     <  2:69263/63‑1 (MQ=255)
                                                                                                            ccTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTgcc                          <  1:756757/121‑1 (MQ=255)
                                                                                                            ccTCGAAGGTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTgcc                          >  2:756757/1‑121 (MQ=255)
                                                                                                                   ggTGCCGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTGCCGCCGCCGATAAAGACCGTAt      >  2:469394/1‑134 (MQ=255)
                                                                                                                       ccGGGGTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTGCCGCCGCCGATAAAGACCGTATGCAc  <  2:60045/134‑1 (MQ=255)
                                                                                                                            |                                                                                                                                
CACCGCCGTGGATGCGCCCGAGCGCCTGGAGGTGCCGGTCGTTGAAGCTCTGGATGCCGATCGACAGCCGGTTGATGCCGCTGGCGCGGTAGCTGGCGAACTTGTCGGCCTCGAAGGTGCCGGGTTTGGCCTCCATCGTGATCTCGGCGTCGGCATCCAGCGGCAGCAGCGCGCGGATATCCGACAGCAGCCGGTCCATGCCTGCCGCCGACAGCAGGCTGGGCGTGCCGCCGCCGATAAAGACCGTATGCAC  >  AM260479/1042283‑1042535

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