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%	123.9 / NA	39	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 (19/20);  total (19/20)

CATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATAT  >  NC_000913/4503887‑4504129                                                                                                                        |                                                                                                                            cATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTa                                                                                                          >  2:11748/1‑139 (MQ=255)      agTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTGCCTAATATTTACTTTAAGGGc                                                                                                     <  2:789879/139‑1 (MQ=255)      agTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGc                                                                                                     <  1:11748/139‑1 (MQ=255)          aaaTATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCtata                                                                                                 >  1:720389/1‑139 (MQ=255)                   aagaagAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTa                                                                                                          >  1:477503/1‑121 (MQ=255)                   aagaagAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTa                                                                                                          <  2:477503/121‑1 (MQ=255)                       agaaATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAAcaca                                                                                    <  1:484850/139‑1 (MQ=255)                           aTGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGaa                                                                                           <  1:924570/128‑1 (MQ=255)                           aTGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGaa                                                                                           >  2:924570/1‑128 (MQ=255)                                      gATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGc                                                                     <  2:84211/139‑1 (MQ=255)                                            aCTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                               <  2:666394/139‑1 (MQ=255)                                             cTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                              <  2:720389/139‑1 (MQ=255)                                                    ggTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGtata                                                       <  2:168323/139‑1 (MQ=255)                                                     gTATTGATAACAATCAATAGTACTGGCTTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATaa                                                      >  2:828468/1‑139 (MQ=255)                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                      <  2:820359/86‑1 (MQ=255)                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                      >  2:141718/1‑86 (MQ=255)                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                      <  1:141718/86‑1 (MQ=255)                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                      >  1:820359/1‑86 (MQ=255)                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                              >  1:323588/1‑126 (MQ=255)                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                              <  2:323588/126‑1 (MQ=255)                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAg                                                 >  1:238430/1‑139 (MQ=255)                                                          gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAg                                                 >  1:238906/1‑139 (MQ=255)                                                                           cTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCAAAACgg                                >  1:359593/1‑139 (MQ=255)                                                                                gTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                      >  2:402221/1‑64 (MQ=255)                                                                                gTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                      <  1:402221/64‑1 (MQ=255)                                                                                   tttGAGGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGctct                        >  1:453581/1‑139 (MQ=255)                                                                                   tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGctct                        >  1:983528/1‑139 (MQ=255)                                                                                   tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGctct                        >  1:225592/1‑139 (MQ=255)                                                                                              tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa              <  1:856396/138‑1 (MQ=255)                                                                                              tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa              >  2:856396/1‑138 (MQ=255)                                                                                                 aaTTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGc          <  2:983528/139‑1 (MQ=255)                                                                                                  attattTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCt         <  2:453581/139‑1 (MQ=255)                                                                                                  attattTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCt         <  2:238906/139‑1 (MQ=255)                                                                                                  attattTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCt         <  1:828468/139‑1 (MQ=255)                                                                                                      tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGa     >  2:49556/1‑139 (MQ=255)                                                                                                       ttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGat    >  2:301431/1‑139 (MQ=255)                                                                                                         cTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGatat  <  2:298469/139‑1 (MQ=255)                                                                                                                 tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa              <  2:310963/119‑1 (MQ=255)                                                                                                                 tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa              >  1:310963/1‑119 (MQ=255)                                                                                                                        |                                                                                                                            CATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATAT  >  NC_000913/4503887‑4504129

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