Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 932,104 A→C intergenic (‑54/‑491) trxB ← / → lrp thioredoxin reductase/DNA‑binding transcriptional dual regulator Lrp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913932,1040AC100.0% 69.7 / NA 20intergenic (‑54/‑491)trxB/lrpthioredoxin reductase/DNA‑binding transcriptional dual regulator Lrp
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (8/12);  total (8/12)

GTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGATTTTGTTAACAATTTGTGCAATCGGCAGCATCGATAAGCAGGTCAAA  >  NC_000913/932026‑932172
                                                                              |                                                                    
gTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCg                                                                    <  1:192454/81‑1 (MQ=255)
            tGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAagag                                                         <  1:192455/80‑1 (MQ=255)
             ggTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGGAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTa                                                       >  2:192486/1‑81 (MQ=255)
               tCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATg                                                     <  2:192456/81‑1 (MQ=255)
                 gTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGa                                                    <  2:192457/80‑1 (MQ=255)
                  tGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACg                                                  <  1:192458/81‑1 (MQ=255)
                       aTGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGAtttt                                             >  2:192495/1‑81 (MQ=255)
                       aTGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGAtttt                                             <  1:192459/81‑1 (MQ=255)
                          agaTCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGtt                                          <  2:192461/81‑1 (MQ=255)
                          agaTCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGt                                           >  2:192483/1‑80 (MQ=255)
                          agaTCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGt                                           <  1:192460/80‑1 (MQ=255)
                               cccATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACa                                      <  2:192462/80‑1 (MQ=255)
                                  aTAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACAAttt                                  <  1:192463/81‑1 (MQ=255)
                                   tAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACAATTtg                                 >  2:192491/1‑81 (MQ=255)
                                     gttgttGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACAATTtgtg                               <  1:192464/81‑1 (MQ=255)
                                      ttgttgGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACAATTtgtg                               >  2:192493/1‑80 (MQ=255)
                                       tgttgGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACAATTTGTGc                              >  1:192492/1‑80 (MQ=255)
                                          tgGCAGACAATGGGCAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACAATTTGTGCAATc                          >  1:192494/1‑81 (MQ=255)
                                                        cAGGATTGTAGGGAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACAATTTGTGCAATCGGCAGCATCGATaa            <  1:192465/81‑1 (MQ=255)
                                                                  gggAATTTACAGCCGTAAAAAAAGAGTATGACGATTTTGTTAACAATTTGTGCAATCGGCAGCATCGATAAGCAGGTCaaa  >  1:192564/1‑81 (MQ=255)
                                                                              |                                                                    
GTTTACTGTGTTTGGTCGTGCCCATGAGATCCCCATAGTTGTTGGCAGACAATGGGCAGGATTGTAGGGAATTTACAGACGTAAAAAAAGAGTATGACGATTTTGTTAACAATTTGTGCAATCGGCAGCATCGATAAGCAGGTCAAA  >  NC_000913/932026‑932172

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