Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,485,462 A→G intergenic (‑33/+21) oxc ← / ← frc oxalyl‑CoA decarboxylase/formyl‑CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,485,4620AG100.0% 78.1 / NA 22intergenic (‑33/+21)oxc/frcoxalyl‑CoA decarboxylase/formyl‑CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/9);  total (13/9)

AATGCTTCAACGATGATATGCATACCATCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGAGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCTATAACCCAGCT  >  NZ_CP009273/2485377‑2485538
                                                                                     |                                                                            
aaTGCTTCAACGATGATATGCATACCATCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCa                                                                          >  2:159325/1‑90 (MQ=255)
       cAACGATGATATGCATACCATCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGccc                                                                   <  1:298871/90‑1 (MQ=255)
        aaCGATGATATGCATACCATCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGcccc                                                                  <  1:18962/90‑1 (MQ=255)
                      tACCATCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGAt                                                    >  2:389580/1‑90 (MQ=255)
                      tACCATCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGAt                                                    <  1:159325/90‑1 (MQ=255)
                          aTCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCg                                                >  2:251889/1‑90 (MQ=255)
                          aTCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCg                                                >  1:269078/1‑90 (MQ=255)
                           tCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGt                                               >  1:108622/1‑90 (MQ=255)
                           tCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGt                                               >  2:22313/1‑90 (MQ=255)
                           tCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGt                                               >  1:315326/1‑90 (MQ=255)
                           tCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGt                                               >  1:108804/1‑90 (MQ=255)
                            cTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTg                                              <  1:251889/90‑1 (MQ=255)
                                      aaGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTc                                    >  1:503180/1‑90 (MQ=255)
                                          tGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTg                                >  2:119955/1‑90 (MQ=255)
                                               tGACATGATTAATAACCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTgcagca                           <  1:119955/90‑1 (MQ=255)
                                                          ataaCCTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCAtcgtc                <  1:389580/90‑1 (MQ=255)
                                                              ccTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCta            <  2:108804/90‑1 (MQ=255)
                                                               cTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCtat           <  1:268743/90‑1 (MQ=255)
                                                               cTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCtat           >  2:171045/1‑90 (MQ=255)
                                                               cTCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCtat           >  2:241674/1‑90 (MQ=255)
                                                                tCGAAATTTTTATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCtata          >  1:105451/1‑90 (MQ=255)
                                                                        tttATAAATTGGGGGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCTATAACCCAGCt  <  1:434962/90‑1 (MQ=255)
                                                                                  ggggGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCtat           <  2:270004/71‑1 (MQ=255)
                                                                                     |                                                                            
AATGCTTCAACGATGATATGCATACCATCTGTCATTTGAAGTTGATCTGACATGATTAATAACCTCGAAATTTTTATAAATTGGGAGCCAGAAGCCCCCGTTAATATCAGATGGCGTGGTTTTGCTTCATTGCAGCAATTTCATCGTCGCTATAACCCAGCT  >  NZ_CP009273/2485377‑2485538

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