Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 800,807 G→T 100% A17S (GCG→TCG)  ybhH → putative isomerase YbhH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913800,8070GT100.0% 74.8 / NA 27A17S (GCG→TCG) ybhHputative isomerase YbhH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (10/17);  total (10/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTAATCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGC  >  NC_000913/800670‑800940
                                                                                                                                          |                                                                                                                                     
ttaATCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAAAGATGCATGATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTc                                                                                                                                  <  1:612351/143‑1 (MQ=255)
 taATCGTGAAGGTAAAGGCATGAGAGATTCGTTGGCGTTAATGTATACATATCGCAATAAGTGATGTGTTTAAAGTGAGGGCATTATGATGAAAAAAATACCCTGCGTGAGGATGCGAGGTGGAACCTCGAGGGGCTCGTTcc                                                                                                                                 <  2:468104/143‑1 (MQ=255)
       ttGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAAGCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTAGGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAAGCTCGAGGGGCTCGTCCCTGTTAg                                                                                                                           >  2:316519/1‑143 (MQ=255)
        tGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTCTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGc                                                                                                                          <  1:325949/143‑1 (MQ=255)
           gCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGCGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGa                                                                                                                       >  1:373394/1‑143 (MQ=255)
                  tCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTCGCGGACCAtttc                                                                                                                >  2:75669/1‑142 (MQ=255)
                     ttAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTAccc                                                                                                             >  2:604625/1‑143 (MQ=255)
                        gagaTTCTTTGACCTTAATGTAACCATAGCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGaa                                                                                                          <  2:394860/143‑1 (MQ=255)
                                         aTGTAACCGTATCGCAATAAGCTATGTTTTTAAATTGAGGGCAATAATAAAAAATAAAGACCCGGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGa                                                                                         <  2:559624/143‑1 (MQ=255)
                                           gTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAGGATTAAACGCCGTGCGATa                                                                                       >  1:292873/1‑143 (MQ=255)
                                                      cGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGAt                                                                             <  1:113079/143‑1 (MQ=255)
                                                        cAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGACAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAGGATCAGACGCAGCGCGATAAAGTATTGATgg                                                                           >  1:222245/1‑143 (MQ=255)
                                                                          aaTTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAAAGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTa                                                         >  1:436778/1‑143 (MQ=255)
                                                                            ttGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAAc                                                       <  2:557862/143‑1 (MQ=255)
                                                                                gggCATGATTATGGAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATc                                                   <  2:23488/143‑1 (MQ=255)
                                                                                     ttattaTGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGaa                                              <  1:322076/143‑1 (MQ=255)
                                                                                          atGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGTGCTCGTTCCTGTTAACGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGa                                         <  1:2369/143‑1 (MQ=255)
                                                                                               aaaaaTACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTa                                    >  1:566466/1‑143 (MQ=255)
                                                                                                        cTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATAGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGgcggcg                           <  2:628107/143‑1 (MQ=255)
                                                                                                        cTGCGAAATGATGCGAGGTGGAAACCCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGgcggcg                           <  1:396554/143‑1 (MQ=255)
                                                                                                                tgaggCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCa                   <  2:212408/143‑1 (MQ=255)
                                                                                                                 gatgCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCAc                  >  1:319377/1‑143 (MQ=255)
                                                                                                                  gtggggggTGGAACCGCGGGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACt                 <  2:310983/137‑1 (MQ=255)
                                                                                                                   tgCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTg                >  1:129239/1‑143 (MQ=255)
                                                                                                                     cGATGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGAc              <  1:237900/143‑1 (MQ=255)
                                                                                                                            gAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTaaa       <  1:23391/143‑1 (MQ=255)
                                                                                                                                 gcGAGGCGCTCGTGCCTGTTAGCGGAACATTTACACGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGc  <  2:236980/142‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                     
TTAATCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGC  >  NC_000913/800670‑800940

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