Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,873,583 A→G S333S (AGT→AGC recF ← DNA replication/repair protein RecF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,873,5830AG100.0% 63.7 / NA 18S333S (AGT→AGCrecFDNA replication/repair protein RecF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/6);  total (12/6)

ATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGA  >  NZ_CP009273/3873499‑3873666
                                                                                    |                                                                                   
ataCTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATc                                                                               <  2:611558/91‑1 (MQ=255)
     tGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCt                                                                          >  2:668252/1‑91 (MQ=255)
     tGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCt                                                                          <  2:623174/91‑1 (MQ=255)
     tGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCt                                                                          >  2:367603/1‑91 (MQ=255)
     tGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCt                                                                          <  2:49619/91‑1 (MQ=255)
     tGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCt                                                                          >  2:588514/1‑91 (MQ=255)
          tAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACaa                                                                     <  2:222140/91‑1 (MQ=255)
                 ttATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCtg                                                              >  1:771978/1‑91 (MQ=255)
                           cTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTc                                                    <  2:124781/91‑1 (MQ=255)
                             tttCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGc                                                  >  2:740049/1‑91 (MQ=255)
                                          aaTCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTg                                     <  2:467801/90‑1 (MQ=255)
                                                             gACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCt                  >  1:356363/1‑91 (MQ=255)
                                                                           gtgtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCa    >  2:46354/1‑91 (MQ=255)
                                                                           gtgtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCa    >  1:316279/1‑91 (MQ=255)
                                                                            tgtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCag   >  2:574625/1‑91 (MQ=255)
                                                                            tgtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCag   >  2:348245/1‑91 (MQ=255)
                                                                             gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCTAGCTCagc  >  2:741046/1‑90 (MQ=255)
                                                                             gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCaga  >  2:198938/1‑91 (MQ=255)
                                                                                    |                                                                                   
ATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGA  >  NZ_CP009273/3873499‑3873666

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