Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,276,251 Δ1 bp coding (932/3744 nt) narG → nitrate reductase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,276,2490G.100.0% 59.4 / NA 15coding (930/3744 nt)narGnitrate reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (5/10);  total (5/10)

TCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCC  >  NZ_CP009273/1276161‑1276332
                                                                                        |                                                                                   
tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATgg                                                                                    >  2:465459/1‑90 (MQ=255)
 tacacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  1:510987/89‑2 (MQ=255)
 cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:70314/90‑2 (MQ=255)
 cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:503159/90‑2 (MQ=255)
 cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:159467/90‑2 (MQ=255)
 cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                                   <  2:869242/90‑2 (MQ=255)
             cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAAtgctg                                                                      >  2:827275/1‑90 (MQ=255)
                 gAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTg                                                                  >  1:22634/1‑90 (MQ=255)
                      cGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggcca                                                                                 >  1:671525/1‑68 (MQ=255)
                      cGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggcca                                                                                 <  2:671525/71‑4 (MQ=255)
                      cGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggcca                                                                                 >  1:905025/1‑68 (MQ=255)
                      cGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggcca                                                                                 <  2:905025/71‑4 (MQ=255)
                       gCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTcc                                                            <  1:856006/90‑1 (MQ=255)
                               gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                                    <  1:122101/90‑1 (MQ=255)
                                         tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                                          <  2:146041/90‑1 (MQ=255)
                                         tggatggcACCTAAACGGGGCCCCGATGCGGCGATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACACCCGAAGcc                                          <  1:667494/90‑1 (MQ=255)
                                              ggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt                                     >  1:955180/1‑90 (MQ=255)
                                                   cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTacac                                <  2:955180/90‑1 (MQ=255)
                                                              aCCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc                     <  2:782574/90‑1 (MQ=255)
                                                              aCCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc                     <  1:465459/90‑1 (MQ=255)
                                                              aCCGATGCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc                     <  1:304205/90‑1 (MQ=255)
                                                              aCCGATGCGGCAATGGCGCTGGCAATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc                     <  1:215587/90‑1 (MQ=255)
                                                                    gCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCt               >  1:742345/1‑90 (MQ=255)
                                                                    gCGGCAATGGCGCTGGCGATGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCt               >  1:403279/1‑90 (MQ=255)
                                                                                 tGGCGATGGCCACGTAATGCTGCGTGAATTCCCCCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGcc  <  1:804349/90‑1 (MQ=255)
                                                                                        |                                                                                   
TCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCC  >  NZ_CP009273/1276161‑1276332

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