Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,504,006 G→A intergenic (‑575/‑52) insG ← / → yjhB IS4 transposase/putative MFS transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,504,0060GA100.0% 94.3 / NA 30intergenic (‑575/‑52)insG/yjhBIS4 transposase/putative MFS transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/15);  total (15/15)

GAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGA  >  NC_000913/4503891‑4504137
                                                                                                                   |                                                                                                                                   
gagTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                              <  1:355321/139‑1 (MQ=255)
             tAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACGATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                                        <  1:1131928/132‑1 (MQ=255)
             tAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACGATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                                        >  2:1131928/1‑132 (MQ=255)
                                           aTTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAatat                                                                                                                          <  1:1718245/84‑1 (MQ=255)
                                           aTTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAatat                                                                                                                          >  2:1718245/1‑84 (MQ=255)
                                                     gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAg                                                                                                                >  1:1775582/1‑84 (MQ=255)
                                                     gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAg                                                                                                                <  2:1775582/84‑1 (MQ=255)
                                                     gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAg                                                         >  1:430197/1‑139 (MQ=255)
                                                          caatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTcc                                                                                                                                 >  1:1702129/1‑62 (MQ=255)
                                                          caatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTcc                                                                                                                                 <  2:1702129/62‑1 (MQ=255)
                                                                       tGGAGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTGAATCCACCACAACGGa                                       >  2:1395821/1‑139 (MQ=255)
                                                                              tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                                      <  1:1820304/101‑1 (MQ=255)
                                                                              tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                                      >  2:1820304/1‑101 (MQ=255)
                                                                                   aGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCtttttt                           >  1:722769/1‑139 (MQ=255)
                                                                                        atatAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa                      <  1:1395821/139‑1 (MQ=255)
                                                                                         tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCAt                     >  1:277732/1‑139 (MQ=255)
                                                                                             attattTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCt                 <  2:722769/139‑1 (MQ=255)
                                                                                                attTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTgg              <  1:970308/139‑1 (MQ=255)
                                                                                                 tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGc                  >  1:255241/1‑134 (MQ=255)
                                                                                                 tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGc                  <  2:255241/134‑1 (MQ=255)
                                                                                                 tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGa             >  1:208043/1‑139 (MQ=255)
                                                                                                 tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGa             >  2:1927040/1‑139 (MQ=255)
                                                                                                  ttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGat            <  1:1927040/139‑1 (MQ=255)
                                                                                                  ttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGat            <  2:208043/139‑1 (MQ=255)
                                                                                                            tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGTCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGc                  >  2:1223691/1‑123 (MQ=255)
                                                                                                            tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCGACCACAACGGAAAGCTCTTTTTTCCGCATGGc                  <  1:1223691/123‑1 (MQ=255)
                                                                                                            tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCttttt                            <  1:2004789/113‑1 (MQ=255)
                                                                                                            tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCttttt                            >  2:2004789/1‑113 (MQ=255)
                                                                                                            tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGa  >  1:224735/1‑139 (MQ=255)
                                                                                                            tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGa  <  2:224735/139‑1 (MQ=255)
                                                                                                                   |                                                                                                                                   
GAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGA  >  NC_000913/4503891‑4504137

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