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.94.7% 67.7 / ‑1.7 19coding (1733/2055 nt)helDDNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (14/4);  total (15/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATC  >  NZ_CP009273/1021576‑1021744
                                                                                   |                                                                                     
ggATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaagcg                                                                                 >  1:941973/1‑87 (MQ=255)
     aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaa                                                                                    >  2:904987/1‑82 (MQ=255)
     aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaa                                                                                    <  1:904987/82‑1 (MQ=255)
          tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCt                                                                      >  1:505011/1‑90 (MQ=255)
                       aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCa                                                         >  2:883423/1‑90 (MQ=255)
                       aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCa                                                         >  2:302116/1‑90 (MQ=255)
                            gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    <  2:783797/90‑1 (MQ=255)
                                 gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                               >  1:213716/1‑90 (MQ=255)
                                      ttCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGa                                           >  2:366308/1‑90 (MQ=255)
                                              cTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATg                                  >  1:302517/1‑90 (MQ=255)
                                              cTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATg                                  >  1:456072/1‑90 (MQ=255)
                                               tGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGc                                 <  2:802090/90‑1 (MQ=255)
                                                ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATg                                            >  2:347801/1‑78 (MQ=255)
                                                ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATg                                            <  1:347801/78‑1 (MQ=255)
                                                     gTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa                           >  2:633798/1‑90 (MQ=255)
                                                                    ggCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGa            >  2:982299/1‑90 (MQ=255)
                                                                       cTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa         >  1:975323/1‑90 (MQ=255)
                                                                       cTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa         >  2:621202/1‑90 (MQ=255)
                                                                       cTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa         <  1:621202/90‑1 (MQ=255)
                                                                         gCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACg       >  1:939881/1‑90 (MQ=255)
                                                                             gCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGtcat   >  1:77483/1‑90 (MQ=255)
                                                                              ccTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGtcatc  >  2:698394/1‑90 (MQ=255)
                                                                                   |                                                                                     
GGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATC  >  NZ_CP009273/1021576‑1021744

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