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% 121.8 / NA 39*234* (TAG→TAA) yhaKredox‑sensitive bicupin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/19);  total (20/19)

ATTTCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAGTTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAGC  >  NC_000913/3254885‑3255155
                                                                                                                                       |                                                                                                                                       
aTTTCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGATCTGCGGTGATGGGGCGTTTATTCGTGACGAGACTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAAATTTGCCTGTCTAattg                                                                                                                                      >  1:659659/1‑139 (MQ=255)
   tCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGttt                                                                                                                                   <  2:1050935/139‑1 (MQ=255)
      tGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTa                                                                                                                                <  2:359607/139‑1 (MQ=255)
        cgcTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                              <  2:659659/139‑1 (MQ=255)
         gcTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                             <  2:828939/139‑1 (MQ=255)
                     tGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                             >  1:329322/1‑127 (MQ=255)
                     tGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAg                                                                                                                             <  2:329322/127‑1 (MQ=255)
                             aaGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTaaaa                                                                                                         >  2:161807/1‑139 (MQ=255)
                              aGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTaaaaa                                                                                                        >  2:209995/1‑139 (MQ=255)
                                   cgcTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTgg                                                                                                   >  1:1075/1‑139 (MQ=255)
                                                  aTGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGc                                                                                    >  2:531796/1‑139 (MQ=255)
                                                      ggCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGt                                                                                <  1:531796/139‑1 (MQ=255)
                                                          tttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAattgt                                                                                                                                     <  2:483579/82‑1 (MQ=255)
                                                          tttATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAattgt                                                                                                                                     >  1:483579/1‑82 (MQ=255)
                                                                cGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCgtg                                                                      <  1:209995/139‑1 (MQ=255)
                                                                        ggCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                              >  2:470210/1‑75 (MQ=255)
                                                                        ggCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                              <  1:470210/75‑1 (MQ=255)
                                                                        ggCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGaa                                                              <  2:836787/139‑1 (MQ=255)
                                                                            aaCATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCgtg                                                                      >  1:401875/1‑127 (MQ=255)
                                                                            aaCATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCgtg                                                                      <  2:401875/127‑1 (MQ=255)
                                                                                   cGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTc                                                   >  1:422363/1‑139 (MQ=255)
                                                                                         ttGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                              <  2:775392/58‑1 (MQ=255)
                                                                                         ttGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACa                                                                                                                              >  1:775392/1‑58 (MQ=255)
                                                                                                  ccccACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGt                                                                                <  1:485885/95‑1 (MQ=255)
                                                                                                  ccccACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGt                                                                                >  2:485885/1‑95 (MQ=255)
                                                                                                          cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAgg                                                                                                                            <  1:654679/43‑1 (MQ=255)
                                                                                                          cgcgcTTTGCTGATAGATTTGCCTGTCTAATTGTTTTTACAgg                                                                                                                            >  2:654679/1‑43 (MQ=255)
                                                                                                                        aGATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTgag              >  2:966426/1‑139 (MQ=255)
                                                                                                                         gATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAg                                                                                     <  2:326411/67‑1 (MQ=255)
                                                                                                                         gATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAg                                                                                     >  1:326411/1‑67 (MQ=255)
                                                                                                                         gATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTgag              >  2:296069/1‑138 (MQ=255)
                                                                                                                         gATTTGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTgag              <  1:296069/138‑1 (MQ=255)
                                                                                                                             tGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                                 <  1:746941/35‑1 (MQ=255)
                                                                                                                             tGCCTGTCTAATTGTTTTTACAGGGAgatgatgat                                                                                                                 >  2:746941/1‑35 (MQ=255)
                                                                                                                             tGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAg                                                                                     <  2:1071689/63‑1 (MQ=255)
                                                                                                                             tGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAg                                                                                     >  1:1071689/1‑63 (MQ=255)
                                                                                                                             tGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAAttt                                   >  2:781817/1‑113 (MQ=255)
                                                                                                                             tGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAAttt                                   <  1:781817/113‑1 (MQ=255)
                                                                                                                             tGCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTgagcgagc         >  1:900150/1‑139 (MQ=255)
                                                                                                                              gCCTGTCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCt        <  1:966426/139‑1 (MQ=255)
                                                                                                                                   tCTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAg   >  2:862621/1‑139 (MQ=255)
                                                                                                                                    cTAATTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAGc  >  1:1024100/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                       
ATTTCATGCGCTTACGCATCATGAAGAGAAAGCGGCGCTGACCTGCGGTGATGGGGCGTTTATTCGTGACGAGGCTAACATTACGCTGGTTGCCGATTCCCCACTGCGCGCTTTGCTGATAGATTTGCCTGTCTAGTTGTTTTTACAGGGAGATGATGATGAGTAAAAAATTGGCCAAAAAGCGCCAGCCGGTGAAGCCCGTGGTGGCGAAAGAACCTGCTCGCACCGCCAAAAATTTTGGCTATGAAGAGATGTTGAGCGAGCTGGAAGC  >  NC_000913/3254885‑3255155

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