Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,646,214 (G)6→7 coding (362/846 nt) sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,646,2081.G100.0% 51.6 / NA 15L119R (CTG→CGG) sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (7/8);  total (7/8)

ACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTA  >  NZ_CP009273/2646072‑2646335
                                                                                                                                         |                                                                                                                               
aCGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgg                                                                                                                               <  2:27295/139‑1 (MQ=255)
                      aaTTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCCCCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTgg                                                                                                           <  2:35920/138‑1 (MQ=255)
                      aaTTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTgg                                                                                                           <  2:266569/138‑1 (MQ=255)
                                                                           ttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                                                        <  2:448030/136‑1 (MQ=255)
                                                                                 gCCCCACGAGCATTGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                                                              >  1:549082/1‑92 (MQ=255)
                                                                                 gCCCCACGAGCATTGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                                                              <  2:549082/92‑1 (MQ=255)
                                                                                                tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCt                                >  2:577130/1‑139 (MQ=255)
                                                                                                         cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCgtg                       >  2:566233/1‑139 (MQ=255)
                                                                                                                  ggTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGt               <  1:67432/138‑1 (MQ=255)
                                                                                                                    tgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAAc            <  1:577130/139‑1 (MQ=255)
                                                                                                                        gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGAt        >  1:187722/1‑139 (MQ=255)
                                                                                                                        gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGAt        <  1:566233/139‑1 (MQ=255)
                                                                                                                        gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGAt        >  2:146639/1‑139 (MQ=255)
                                                                                                                             aGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGa         <  1:416964/133‑1 (MQ=255)
                                                                                                                             aGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGa         >  2:416964/1‑133 (MQ=255)
                                                                                                                              gtgtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGtatta  >  2:580103/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                               
ACGTTCGCCGTGGCGATGCGTGAATTAGGCGTTAACCAGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTA  >  NZ_CP009273/2646072‑2646335

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