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_000913	4,504,006	0	G	A	100.0%	134.7 / NA	43	intergenic (‑575/‑52)	insG/yjhB	IS4 transposase/putative MFS transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/21);  total (22/21)

TATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGAT  >  NC_000913/4503871‑4504127                                                                                                                                        |                                                                                                                          tATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTc                                                                                                                        >  1:431217/1‑139 (MQ=255)     aGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTa                                                                                                                    <  2:643926/139‑1 (MQ=255)                          aaaTATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCtata                                                                                               >  2:1388565/1‑139 (MQ=255)                                    agaagaAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTcc                                                                                                                       >  2:819124/1‑105 (MQ=255)                                    agaagaAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTcc                                                                                                                       <  1:819124/105‑1 (MQ=255)                                       agaaATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAAcaca                                                                                  <  2:1320462/139‑1 (MQ=255)                                                            aCTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                             <  1:1388565/139‑1 (MQ=255)                                                                     gTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                    >  2:319725/1‑91 (MQ=255)                                                                     gTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                    <  1:319725/91‑1 (MQ=255)                                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAg                                               >  2:729648/1‑139 (MQ=255)                                                                               caatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAAt                                          >  1:720982/1‑139 (MQ=255)                                                                                  tcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATcca                                       >  1:186600/1‑139 (MQ=255)                                                                                     atAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                              <  2:361673/81‑1 (MQ=255)                                                                                     atAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                              >  1:361673/1‑81 (MQ=255)                                                                                           cTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGt                                              <  2:371322/123‑1 (MQ=255)                                                                                           cTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGt                                              >  1:371322/1‑123 (MQ=255)                                                                                            tGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACaa                                                <  1:733318/120‑1 (MQ=255)                                                                                            tGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACaa                                                >  2:733318/1‑120 (MQ=255)                                                                                              gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                              >  2:1136075/1‑72 (MQ=255)                                                                                              gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                              <  1:1136075/72‑1 (MQ=255)                                                                                              gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGt                                                        <  1:1255870/110‑1 (MQ=255)                                                                                              gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGt                                                        >  2:1255870/1‑110 (MQ=255)                                                                                                gTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                    >  1:1301365/1‑64 (MQ=255)                                                                                                gTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                    <  2:1301365/64‑1 (MQ=255)                                                                                                   tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                            >  1:488165/1‑101 (MQ=255)                                                                                                   tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                            <  2:488165/101‑1 (MQ=255)                                                                                                   tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATccacca                                    <  2:1329846/125‑1 (MQ=255)                                                                                                   tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATccacca                                    >  1:1329846/1‑125 (MQ=255)                                                                                                        aGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                            <  1:1464129/96‑1 (MQ=255)                                                                                                        aGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                            >  2:1464129/1‑96 (MQ=255)                                                                                                        aGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCtttttt                 >  1:318986/1‑139 (MQ=255)                                                                                                           cAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCg              <  2:186600/139‑1 (MQ=255)                                                                                                            aaTATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGc             <  2:318986/139‑1 (MQ=255)                                                                                                            aaTATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGc             >  1:327436/1‑139 (MQ=255)                                                                                                              tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa            >  2:922151/1‑138 (MQ=255)                                                                                                              tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa            <  1:922151/138‑1 (MQ=255)                                                                                                              tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGATAATCCACCACAACGGAAAGCTCTTTTTTCCGCAt           >  1:910292/1‑139 (MQ=255)                                                                                                                  attattTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCt       <  1:729648/139‑1 (MQ=255)                                                                                                                   ttattTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCtt      >  2:747966/1‑139 (MQ=255)                                                                                                                      tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGc                                                                   >  2:262271/1‑75 (MQ=255)                                                                                                                      tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGc                                                                   <  1:262271/75‑1 (MQ=255)                                                                                                                      tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGa   <  2:431217/139‑1 (MQ=255)                                                                                                                       ttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGat  <  1:747966/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                          TATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGAT  >  NC_000913/4503871‑4504127

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