Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,255,020 G→A *234* (TAG→TAA yhaK → redox‑sensitive bicupin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,255,0200GA100.0% 86.8 / NA 28*234* (TAG→TAAyhaKredox‑sensitive bicupin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/16);  total (12/16)

GCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAGTTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATG  >  NC_000913/3254894‑3255138
                                                                                                                              |                                                                                                                      
gcTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                            <  2:867159/139‑1 (MQ=255)
  ttACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAggg                                                                                                          <  2:991931/139‑1 (MQ=255)
          catGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatg                                                                                                  <  1:204787/139‑1 (MQ=255)
           atGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatga                                                                                                 <  1:576128/139‑1 (MQ=255)
            tGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                <  2:76887/139‑1 (MQ=255)
                        ggcgCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAattgt                                                                                                                    >  2:345992/1‑107 (MQ=255)
                        ggcgCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAattgt                                                                                                                    <  1:345992/107‑1 (MQ=255)
                           gcTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                            >  1:567305/1‑112 (MQ=255)
                           gcTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                            <  2:672768/112‑1 (MQ=255)
                           gcTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                            <  2:567305/112‑1 (MQ=255)
                           gcTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                            >  1:672768/1‑112 (MQ=255)
                                         aTGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                             <  2:950959/97‑1 (MQ=255)
                                         aTGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                             >  1:950959/1‑97 (MQ=255)
                                                   tATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatg                                                                                                  >  1:1053106/1‑98 (MQ=255)
                                                   tATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatg                                                                                                  <  2:1053106/98‑1 (MQ=255)
                                                   tATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGcc                                                         >  1:799044/1‑139 (MQ=255)
                                                          gACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTAc                                                                                                              >  2:989339/1‑79 (MQ=255)
                                                          gACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTAc                                                                                                              <  1:989339/79‑1 (MQ=255)
                                                               ggCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                <  2:299318/88‑1 (MQ=255)
                                                               ggCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                >  1:299318/1‑88 (MQ=255)
                                                                                         ccccACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                             <  2:1085699/49‑1 (MQ=255)
                                                                                         ccccACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                             >  1:1085699/1‑49 (MQ=255)
                                                                                         ccccACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCgtggt                                                   >  2:249647/1‑107 (MQ=255)
                                                                                         ccccACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCgtggt                                                   <  1:249647/107‑1 (MQ=255)
                                                                                                   cgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGa                                              <  2:240262/102‑1 (MQ=255)
                                                                                                   cgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGa                                              >  1:240262/1‑102 (MQ=255)
                                                                                                    gcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGaa        >  1:202575/1‑139 (MQ=255)
                                                                                                          cTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATg  <  1:573028/139‑1 (MQ=255)
                                                                                                                    tGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGcc                                                                      >  1:611197/1‑61 (MQ=255)
                                                                                                                    tGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGcc                                                                      <  2:611197/61‑1 (MQ=255)
                                                                                                                              |                                                                                                                      
GCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAGTTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATG  >  NC_000913/3254894‑3255138

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