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%	161.2 / NA	50	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 (25/25);  total (25/25)

GAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGT  >  NC_000913/4503891‑4504131                                                                                                                    |                                                                                                                              gagTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                        <  1:778326/139‑1 (MQ=255)  agTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGc                                                                                                       <  1:734737/139‑1 (MQ=255)       aaTATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAACAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                                  <  2:820523/139‑1 (MQ=255)              tAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATa                                                                                           >  1:406725/1‑139 (MQ=255)                agaagaAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAac                                                                                         >  2:905332/1‑139 (MQ=255)                       aTGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGaa                                                                                  <  2:215887/139‑1 (MQ=255)                               aaaGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAatat                                                                          >  2:322662/1‑139 (MQ=255)                                 aGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATgg                                                                        <  1:64629/139‑1 (MQ=255)                                        aCTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATGTTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                 <  1:1045246/139‑1 (MQ=255)                                        aCTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                 >  2:171790/1‑139 (MQ=255)                                         cTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                                <  2:489649/139‑1 (MQ=255)                                         cTGATTAAGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                                <  1:401291/139‑1 (MQ=255)                                            aTTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGc                                                                                                       >  1:1026522/1‑97 (MQ=255)                                            aTTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGc                                                                                                       <  2:1026522/97‑1 (MQ=255)                                            aTTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGc                                                                       >  1:528658/1‑129 (MQ=255)                                            aTTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGc                                                                       <  2:528658/129‑1 (MQ=255)                                                   aTTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAAc                                                      >  1:912686/1‑139 (MQ=255)                                                      gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                        <  1:843816/86‑1 (MQ=255)                                                      gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                        >  2:843816/1‑86 (MQ=255)                                                      gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                 <  2:361238/125‑1 (MQ=255)                                                      gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                 >  1:361238/1‑125 (MQ=255)                                                      gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAg                                                   >  2:949079/1‑139 (MQ=255)                                                       aTAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGt                                                  <  2:912686/139‑1 (MQ=255)                                                           caatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTa                                                                                                            >  1:449471/1‑77 (MQ=255)                                                           caatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTa                                                                                                            <  2:449471/77‑1 (MQ=255)                                                            aatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATc                                             <  1:322662/139‑1 (MQ=255)                                                                      aCTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATc                                             <  2:256841/129‑1 (MQ=255)                                                                      aCTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATc                                             >  1:256841/1‑129 (MQ=255)                                                                       cTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACgg                                  >  1:586584/1‑139 (MQ=255)                                                                         ggCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGaa                                >  1:199784/1‑139 (MQ=255)                                                                          gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGaaa                               >  2:762438/1‑139 (MQ=255)                                                                          gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGaaa                               >  2:625861/1‑139 (MQ=255)                                                                               tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGctct                          >  1:992664/1‑139 (MQ=255)                                                                               tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGctct                          >  2:931631/1‑139 (MQ=255)                                                                                     gACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                 <  2:622232/94‑1 (MQ=255)                                                                                     gACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                 >  1:622232/1‑94 (MQ=255)                                                                                         atatAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa                <  2:992664/139‑1 (MQ=255)                                                                                         atatAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa                <  1:931631/139‑1 (MQ=255)                                                                                          tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATcca                                           <  2:377156/111‑1 (MQ=255)                                                                                          tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATcca                                           >  1:377156/1‑111 (MQ=255)                                                                                                   ttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGat      <  2:283847/139‑1 (MQ=255)                                                                                                   ttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGat      <  2:199784/139‑1 (MQ=255)                                                                                                    tCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCg                  >  2:398110/1‑126 (MQ=255)                                                                                                    tCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCg                  <  1:398110/126‑1 (MQ=255)                                                                                                    tCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGata     >  1:61682/1‑139 (MQ=255)                                                                                                       ggATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGt  >  2:936687/1‑139 (MQ=255)                                                                                                         atatTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                 <  2:233334/74‑1 (MQ=255)                                                                                                         atatTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                 >  1:233334/1‑74 (MQ=255)                                                                                                             tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa                >  1:913453/1‑119 (MQ=255)                                                                                                             tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa                <  2:913453/119‑1 (MQ=255)                                                                                                                    |                                                                                                                              GAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGT  >  NC_000913/4503891‑4504131

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