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%	219.8 / NA	68	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 (33/35);  total (33/35)

TCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGA  >  NC_000913/4503873‑4504137                                                                                                                                      |                                                                                                                                    tCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTccc                                                                                                                                >  1:617828/1‑139 (MQ=255)             ttCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACtt                                                                                                                    >  2:780249/1‑139 (MQ=255)              tCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACttt                                                                                                                   <  1:885913/139‑1 (MQ=255)                   gagTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                              <  1:780249/139‑1 (MQ=255)                    agTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGc                                                                                                             <  1:1022878/139‑1 (MQ=255)                    agTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGc                                                                                                             <  1:1030711/139‑1 (MQ=255)                      tGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAatat                                                                                                                          >  2:1351880/1‑124 (MQ=255)                      tGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAatat                                                                                                                          <  1:1351880/124‑1 (MQ=255)                         aaTATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAAGAATCAATAGGACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                                        <  1:203383/139‑1 (MQ=255)                         aaTATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATGTACTTTAAGGGCTatat                                                                                                        <  1:756201/139‑1 (MQ=255)                                tAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAatat                                                                                                                          >  2:1356534/1‑114 (MQ=255)                                tAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAatat                                                                                                                          <  1:1356534/114‑1 (MQ=255)                                tAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACt                                                                                                                     <  1:52157/119‑1 (MQ=255)                                tAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACt                                                                                                                     >  2:52157/1‑119 (MQ=255)                                  agaagaAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAAt                                                                                                  >  1:189202/1‑136 (MQ=255)                                  agaagaAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAAt                                                                                                  <  2:189202/136‑1 (MQ=255)                                    aagaaATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAacac                                                                                             <  2:79092/139‑1 (MQ=255)                                       aaaTGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAgg                                                                                          >  2:636513/1‑139 (MQ=255)                                                    gATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGc                                                                             >  1:383959/1‑139 (MQ=255)                                                       tGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAatat                                                                                                                          <  2:457299/91‑1 (MQ=255)                                                       tGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAatat                                                                                                                          >  1:457299/1‑91 (MQ=255)                                                       tGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAAc                                                                          >  2:1375227/1‑139 (MQ=255)                                                          aCTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                       <  1:1375227/139‑1 (MQ=255)                                                           cTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                                      <  1:636513/139‑1 (MQ=255)                                                           cTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                                      <  2:5393/139‑1 (MQ=255)                                                           cTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCa                                                                      <  1:230926/139‑1 (MQ=255)                                                              aTTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAg                                                                                                                <  2:804901/94‑1 (MQ=255)                                                              aTTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAg                                                                                                                >  1:804901/1‑94 (MQ=255)                                                               ttAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGt                                                                  >  1:306512/1‑139 (MQ=255)                                                                        gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                              <  2:866016/86‑1 (MQ=255)                                                                        gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                       >  1:504886/1‑125 (MQ=255)                                                                        gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                       <  2:504886/125‑1 (MQ=255)                                                                        gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAg                                                         >  2:74303/1‑139 (MQ=255)                                                                        gATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAg                                                         >  2:1248881/1‑139 (MQ=255)                                                                        gATAACAATCAATAGTACTGACGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAggg                                                                                                              >  1:866016/1‑86 (MQ=255)                                                                             caatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAg                                                                                                                <  2:1418528/79‑1 (MQ=255)                                                                             caatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAg                                                                                                                >  1:1418528/1‑79 (MQ=255)                                                                             caatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAAt                                                    >  1:1205326/1‑139 (MQ=255)                                                                              aatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                       >  1:1128476/1‑119 (MQ=255)                                                                              aatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                       <  2:1128476/119‑1 (MQ=255)                                                                              aatcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATc                                                   <  1:1248881/139‑1 (MQ=255)                                                                                tcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                                        >  2:1330550/1‑84 (MQ=255)                                                                                tcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                                        <  1:1330550/84‑1 (MQ=255)                                                                                tcaatAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATcca                                                 >  1:672771/1‑139 (MQ=255)                                                                                        aCTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                                        >  2:1384103/1‑76 (MQ=255)                                                                                        aCTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTatat                                                                                                        <  1:1384103/76‑1 (MQ=255)                                                                                            gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAAcaca                                                                                            >  2:688652/1‑84 (MQ=255)                                                                                            gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAAcaca                                                                                            >  2:269366/1‑84 (MQ=255)                                                                                            gCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAAcaca                                                                                            <  1:269366/84‑1 (MQ=255)                                                                                            gCGTATTTGAAGAAAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAAcaca                                                                                            <  1:688652/84‑1 (MQ=255)                                                                                              gTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGAGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACATCGGAAAGc                                   >  2:871981/1‑139 (MQ=255)                                                                                                 tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCaa                                                                           <  1:1350188/96‑1 (MQ=255)                                                                                                 tttGAAGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCaa                                                                           >  2:1350188/1‑96 (MQ=255)                                                                                                      aGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                       <  1:328059/95‑1 (MQ=255)                                                                                                      aGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGc                                                                       >  2:328059/1‑95 (MQ=255)                                                                                                      aGACAATATAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCtttttt                           >  1:100002/1‑139 (MQ=255)                                                                                                           atatAATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCa                      <  2:672771/139‑1 (MQ=255)                                                                                                            tataATTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCAt                     >  2:1421195/1‑139 (MQ=255)                                                                                                               aaTTATTTCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGc                  <  2:100002/139‑1 (MQ=255)                                                                                                                    tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGc                                   >  1:655117/1‑117 (MQ=255)                                                                                                                    tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGc                                   <  2:655117/117‑1 (MQ=255)                                                                                                                    tttCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGa             >  2:766381/1‑139 (MQ=255)                                                                                                                      tCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACg                                         <  1:1299419/109‑1 (MQ=255)                                                                                                                      tCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACg                                         >  2:1299419/1‑109 (MQ=255)                                                                                                                      tCTGGATATTGTTGAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGata           <  1:74303/139‑1 (MQ=255)                                                                                                                               tgttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGa  <  2:306512/139‑1 (MQ=255)                                                                                                                                 ttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCTACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGat            >  2:933323/1‑127 (MQ=255)                                                                                                                                 ttgAGACTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGat            <  1:933323/127‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    TCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTACTGATTAGGTATTGATAACAATCAATAGTACTGGCGTATTTGAAGACAATATAATTATTTCTGGATATTGTTGAGGCTCCCTAATATTTACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGA  >  NC_000913/4503873‑4504137

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