Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,882,413 A→C V349G (GTT→GGT)  cysI ← assimilatory sulfite reductase (NADPH) hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,882,4130AC100.0% 22.1 / NA 8V349G (GTT→GGT) cysIassimilatory sulfite reductase (NADPH) hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (2/6);  total (2/6)

CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGC  >  NZ_CP009273/2882326‑2882482
                                                                                       |                                                                     
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                     >  1:109006/1‑90 (MQ=255)
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                     >  2:197866/1‑90 (MQ=255)
            ggACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGa                                                         >  2:6797/1‑90 (MQ=255)
                       aTAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTcc                                              <  2:117544/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                       <  2:109006/90‑1 (MQ=255)
                                gATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                       <  1:335368/88‑1 (MQ=255)
                                gATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                       >  2:335368/1‑88 (MQ=255)
                                   gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                  <  2:6617/90‑1 (MQ=255)
                                             tCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                        <  2:8893/90‑1 (MQ=255)
                                                                   aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc  <  1:148595/90‑1 (MQ=255)
                                                                                       |                                                                     
CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGC  >  NZ_CP009273/2882326‑2882482

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