Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,834,119 G→T intergenic (+68/‑46) osmC → / → USA300HOU_RS09140 Peroxiredoxin OsmC/Soluble hydrogenase 42 kDa subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,834,1190GT100.0% 67.6 / NA 20intergenic (+68/‑46)osmC/USA300HOU_RS09140Peroxiredoxin OsmC/Soluble hydrogenase 42 kDa subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (9/11);  total (9/11)

TAATAATTGCATGATTTCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCAGAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCATATAAATGTATTATCATCAACCGTTGTTATTAACACCTGGCCCAACCCCTGTACCTGATGCCATTATGAGAGAA  >  USA300TCH1516_ALE/1833976‑1834233
                                                                                                                                               |                                                                                                                  
taataaTTGCATGATTTCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATaaaa                                                                                                                <  2:6246/148‑1 (MQ=255)
taataaTTGCATGATTTCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATaaaa                                                                                                                <  2:56439/148‑1 (MQ=255)
taataaTTGCATGATTTCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATaaaa                                                                                                                <  2:44419/148‑1 (MQ=255)
taataaTTGCATGATTTCAAATGCTGTAAGAAATAATGTGAACATAATAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCAACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATaaaa                                                                                                                <  2:41650/148‑1 (MQ=255)
  ataatTGCATGATTTCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAAtt                                                                                                              <  1:90178/148‑1 (MQ=255)
          aTGATTTCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTaaca                                                                                                      >  2:206608/1‑146 (MQ=255)
          aTGATTTCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTaaaa                                                                                                      <  1:68695/148‑1 (MQ=255)
              tttCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGt                                                                                                  <  2:205250/148‑1 (MQ=255)
                       cTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCtta                                                                                         <  2:157967/148‑1 (MQ=255)
                       cTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCAATCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCtta                                                                                         <  1:149882/148‑1 (MQ=255)
                        tGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCttat                                                                                        <  2:231715/148‑1 (MQ=255)
                        tGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCttat                                                                                        <  1:140795/148‑1 (MQ=255)
                                      tgGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCata                                                                          >  2:134135/1‑148 (MQ=255)
                                                                     cAAATAAAATCGTTAAGTGATAACGTTGGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCATATAAATGTATTATCATCAACCGTTGTTATTaa                                           >  2:222332/1‑148 (MQ=255)
                                                                                                  ttCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTTAAATCTTGTTTAGAGGAGTCATATAAATTTATGATCATCAACCGGTGTTATTAACACCTGGCCCAACCCCTGTACCTGATGcc              >  1:324410/1‑148 (MQ=255)
                                                                                                       cAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCATATAAATGTATTATCATCAACCGTTGTTATTATCACCTGGCCCAACCCCTGTACCTGATGCCATTAt         >  1:354710/1‑148 (MQ=255)
                                                                                                       cAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCATATAAATGTATTATCATCAACCGTTGTTATTAACACCTGGCCCAACCCCTGTACCTGATGCCATTAt         >  1:101671/1‑148 (MQ=255)
                                                                                                       cAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCATATAAATGTATTATCATCAACCGTTGTTATTAACACCTGGCCCAACCCCTGTACCTGATGCCATTAt         >  1:126236/1‑148 (MQ=255)
                                                                                                           ataGCGTTATCACTTTTTTAGGCTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCATATAAATGTATTATCATCAACCGTTGTTATTAACACCTGGCCCAACCCCTGTACCTGATGCCATTATgaga     >  2:44780/1‑148 (MQ=255)
                                                                                                              gCGTTATCACTTTTTTTTTTTTTTAAATTTTCATAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCATATAAATGTATTATCATCAACCGTTGTTATTAACACCTGGCCCAACCCCTGTACCTGATGCCATTATGAGAGaa  >  2:150767/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                  
TAATAATTGCATGATTTCAAATGCTGTAAGAAATAATGTGGACATAAAAATCTATCCCATCATTCAGGCCAAATAAAATCGTTAAGTGATAACGTTTGTTCGACAATATAGCGTTATCACTTTTTTAGGCTTTTAAATTTTCAGAAAATTTGCTAAAATAGTGAAATCTTATTTATAGGAGGCATATAAATGTATTATCATCAACCGTTGTTATTAACACCTGGCCCAACCCCTGTACCTGATGCCATTATGAGAGAA  >  USA300TCH1516_ALE/1833976‑1834233

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