Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,021,663 (A)5→4 coding (1737/2055 nt) helD → DNA helicase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,021,6590A.100.0% 62.1 / NA 16coding (1733/2055 nt)helDDNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/6);  total (10/6)

TTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACG  >  NZ_CP009273/1021571‑1021739
                                                                                        |                                                                                
ttGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaa                                                                                 >  1:266046/1‑90 (MQ=255)
  gCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaa                                                                               <  2:566024/90‑1 (MQ=255)
          aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAac                                                                      >  1:643164/1‑90 (MQ=255)
               tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCt                                                                 >  1:800142/1‑90 (MQ=255)
                     tATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGa                                                           >  1:425015/1‑90 (MQ=255)
                                 gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                               <  2:861220/90‑1 (MQ=255)
                                 gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                               <  2:266046/90‑1 (MQ=255)
                                                aTCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTc                                <  2:36048/90‑1 (MQ=255)
                                                     ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                           >  2:220858/1‑90 (MQ=255)
                                                          gTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa                      >  2:449803/1‑90 (MQ=255)
                                                          gTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa                      >  1:862053/1‑90 (MQ=255)
                                                           ttACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaaa                     >  2:373145/1‑90 (MQ=255)
                                                                  cacaTGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAAc              <  2:643164/90‑1 (MQ=255)
                                                                      tGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGc          <  2:792629/90‑1 (MQ=255)
                                                                      tGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGc          <  1:628825/90‑1 (MQ=255)
                                                                       gAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCg         >  2:115885/1‑90 (MQ=255)
                                                                            cTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa    >  2:133570/1‑90 (MQ=255)
                                                                              gCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACg  >  2:441461/1‑90 (MQ=255)
                                                                                        |                                                                                
TTGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACG  >  NZ_CP009273/1021571‑1021739

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