Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,171,725 G→C D185H (GAC→CAC)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,171,7250GC100.0% 73.7 / NA 23D185H (GAC→CAC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (14/9);  total (14/9)

TTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCC  >  NZ_CP009273/4171641‑4171809
                                                                                    |                                                                                    
tttGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTc                                                                                 <  1:605128/90‑1 (MQ=255)
            aaaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCg                                                                     <  2:441917/90‑1 (MQ=255)
                    aaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACaa                                                             <  1:331764/90‑1 (MQ=255)
                        cACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTg                                                         >  1:331360/1‑90 (MQ=255)
                         aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  2:57937/1‑90 (MQ=255)
                         aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  2:441500/1‑90 (MQ=255)
                                      aGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGa                                           >  2:287296/1‑90 (MQ=255)
                                          cTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       <  2:137886/90‑1 (MQ=255)
                                                 aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  1:125207/1‑90 (MQ=255)
                                                 aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  2:230704/1‑90 (MQ=255)
                                                 aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  2:385017/1‑90 (MQ=255)
                                                   gcgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATAGACCGTCGCCGTaaa                              >  2:551817/1‑90 (MQ=255)
                                                      cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAActg                           >  2:248756/1‑90 (MQ=255)
                                                              cccTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGAc                   <  1:230704/90‑1 (MQ=255)
                                                              cccTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGAc                   <  2:440230/90‑1 (MQ=255)
                                                                  tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACcatc               >  1:163093/1‑90 (MQ=255)
                                                                  tACCGTGGTCCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTaaa                              >  2:352939/1‑75 (MQ=255)
                                                                  tACCGTGGTCCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTaaa                              <  1:352939/75‑1 (MQ=255)
                                                                     cGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCAtt            >  2:592096/1‑90 (MQ=255)
                                                                      gTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGCCAACCTGCTCGCACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCCTTc           >  2:281533/1‑90 (MQ=255)
                                                                        ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCg                     <  1:483979/78‑1 (MQ=255)
                                                                        ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCg                     >  2:483979/1‑78 (MQ=255)
                                                                           tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATAGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgc      <  1:551817/90‑1 (MQ=255)
                                                                               ggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGccc  <  2:377782/90‑1 (MQ=255)
                                                                               ggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGccc  <  2:192325/90‑1 (MQ=255)
                                                                               ggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGccc  <  1:385017/90‑1 (MQ=255)
                                                                                    |                                                                                    
TTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCC  >  NZ_CP009273/4171641‑4171809

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