Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,176,858 C→T *141* (TAG→TAA)  yqiB ← DUF1249 protein YqiB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,8580CT100.0% 188.6 / NA 59*141* (TAG→TAA) yqiBDUF1249 protein YqiB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/27);  total (32/27)

GTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAG  >  NC_000913/3176723‑3176990
                                                                                                                                       |                                                                                                                                    
gtTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaa                                                                                                                                   >  1:414991/1‑139 (MQ=255)
   tACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaataa                                                                                                                                <  2:372258/139‑1 (MQ=255)
              gggCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCATATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGc                                                                                                                     <  2:27378/139‑1 (MQ=255)
                              tgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaa                                                                                                                                   <  2:1170139/109‑1 (MQ=255)
                              tgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaa                                                                                                                                   >  1:1170139/1‑109 (MQ=255)
                              tgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaa                                                                                                                                   >  1:54283/1‑109 (MQ=255)
                              tgCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTtaa                                                                                                                                   <  2:54283/109‑1 (MQ=255)
                                     aGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGt                                                                                              >  1:344703/1‑139 (MQ=255)
                                               cGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgtg                                                                                                       <  2:405514/120‑1 (MQ=255)
                                               cGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAtgtg                                                                                                       >  1:405514/1‑120 (MQ=255)
                                                 gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATaaa                                                                                                                               <  2:820006/94‑1 (MQ=255)
                                                 gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATaaa                                                                                                                               >  1:820006/1‑94 (MQ=255)
                                                 gTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGt                                                                                  >  1:519114/1‑139 (MQ=255)
                                                                 tGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              <  1:1082783/95‑1 (MQ=255)
                                                                 tGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              >  2:1082783/1‑95 (MQ=255)
                                                                       tGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  2:1314380/1‑91 (MQ=255)
                                                                       tGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  1:1314380/91‑1 (MQ=255)
                                                                        ggCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              <  1:296237/88‑1 (MQ=255)
                                                                        ggCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCt                                                                                                              >  2:296237/1‑88 (MQ=255)
                                                                              aCCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  2:583383/84‑1 (MQ=255)
                                                                              aCCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  1:583383/1‑84 (MQ=255)
                                                                                   ccagAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAAt                                                                                                                        >  2:256557/1‑67 (MQ=255)
                                                                                   ccagAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTtcgtc                                                >  2:805101/1‑139 (MQ=255)
                                                                                   ccagAGGAAGGGTTAACACGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAAt                                                                                                                        <  1:256557/67‑1 (MQ=255)
                                                                                           aGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGAt                                        >  2:1059352/1‑139 (MQ=255)
                                                                                            gggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATg                                       >  2:745624/1‑139 (MQ=255)
                                                                                             ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAt                                                                                                          <  1:1129980/71‑1 (MQ=255)
                                                                                             ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCAt                                                                                                          >  2:1129980/1‑71 (MQ=255)
                                                                                             ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGc                                      >  1:116502/1‑139 (MQ=255)
                                                                                             ggTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGc                                      <  1:522467/139‑1 (MQ=255)
                                                                                                 aaCAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACt                                  >  1:348509/1‑139 (MQ=255)
                                                                                                   cAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  2:259669/63‑1 (MQ=255)
                                                                                                   cAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  1:259669/1‑63 (MQ=255)
                                                                                                      gCTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTtta                             <  2:116502/139‑1 (MQ=255)
                                                                                                       cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGa                                                                                                                          >  1:887981/1‑45 (MQ=255)
                                                                                                       cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGa                                                                                                                          <  2:887981/45‑1 (MQ=255)
                                                                                                       cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  1:1126389/59‑1 (MQ=255)
                                                                                                       cTTTCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  2:1126389/1‑59 (MQ=255)
                                                                                                        tttCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                           <  1:1059352/139‑1 (MQ=255)
                                                                                                          tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  2:348509/139‑1 (MQ=255)
                                                                                                          tCCAAATGGTGTCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  1:745624/139‑1 (MQ=255)
                                                                                                                 ggTGTCCTTAGGTTTCACGACGTTACTAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATc                  >  1:587772/1‑139 (MQ=255)
                                                                                                                    gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                           >  2:1122234/1‑127 (MQ=255)
                                                                                                                    gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttatt                           <  1:1122234/127‑1 (MQ=255)
                                                                                                                    gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttat                            <  2:411491/126‑1 (MQ=255)
                                                                                                                    gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTttat                            >  1:411491/1‑126 (MQ=255)
                                                                                                                    gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAg                         <  1:1045868/129‑1 (MQ=255)
                                                                                                                    gtCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTCTTCGTCGCGTTGATGCAACTTTttatttg                         >  2:1045868/1‑127 (MQ=255)
                                                                                                                     tCCTTAGGTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATAc              >  2:220048/1‑139 (MQ=255)
                                                                                                                            gTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                <  2:277345/50‑1 (MQ=255)
                                                                                                                            gTTTCACGACGTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACa                                                                                                >  1:277345/1‑50 (MQ=255)
                                                                                                                             tttCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  2:512540/37‑1 (MQ=255)
                                                                                                                             tttCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  1:512540/1‑37 (MQ=255)
                                                                                                                               tCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            <  1:672298/35‑1 (MQ=255)
                                                                                                                               tCACGACGTTAATAAACCGGAATCGCCATCGCTcc                                                                                                            >  2:672298/1‑35 (MQ=255)
                                                                                                                                 acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTaaa                                                                                                  <  2:1098903/43‑1 (MQ=255)
                                                                                                                                 acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTaaa                                                                                                  >  1:1098903/1‑43 (MQ=255)
                                                                                                                                 acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAg  >  2:291286/1‑139 (MQ=255)
                                                                                                                                 acgacgTTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAg  >  2:692568/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                    
GTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGAATCGCCATCGCTCCATGTGCTAAACAGTATCGCAACCAGTCCGCTAAAAACTGATTAATTTGATGCTTTTCGTCGCGTTGATGCAACTTTTTATTAGGATAATCATACCGCGCTTTGAAG  >  NC_000913/3176723‑3176990

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