Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	536,986	(G)5→6	intergenic (+38/‑149)	clpX → / → lon	ATP‑dependent protease ATP‑binding subunit ClpX/endopeptidase La

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	536,981	1	.	G	100.0%	132.7 / NA	35	intergenic (+33/‑154)	clpX/lon	ATP‑dependent protease ATP‑binding subunit ClpX/endopeptidase La
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (29/6);  total (29/6)

GGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑GGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTT  >  NC_017635/536933‑537026                                                  |                                              ggCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCaaaaat                                               >  2:1778853/1‑49 (MQ=255) ggCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑g                                              <  2:923408/50‑1 (MQ=255)  gCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑gg                                             >  1:1969682/1‑50 (MQ=255)   cATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑ggg                                            <  1:943127/50‑1 (MQ=255)   cATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑ggg                                            <  2:1194390/50‑1 (MQ=255)       tGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑gggggga                                        >  1:1685589/1‑48 (MQ=255)       tGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑gggggga                                        >  1:1874425/1‑48 (MQ=255)       tGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑gggggga                                        >  2:1130458/1‑48 (MQ=255)         gTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑gggggga                                        <  2:194015/48‑3 (MQ=255)         gTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑gggggga                                        <  2:698450/48‑3 (MQ=255)            aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                    >  1:93945/1‑50 (MQ=255)            aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                    <  2:929409/50‑1 (MQ=255)            aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                    >  1:711515/1‑50 (MQ=255)            aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                    >  1:641886/1‑50 (MQ=255)            aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTa                                    >  2:1990119/1‑50 (MQ=255)             ataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAt                                   >  2:1757088/1‑50 (MQ=255)             ataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAt                                   >  1:1309241/1‑50 (MQ=255)              taatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtc                                  >  1:1034904/1‑50 (MQ=255)              taatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtc                                  >  1:1134929/1‑50 (MQ=255)               aatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtct                                 >  1:399282/1‑50 (MQ=255)               aatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtct                                 >  1:952546/1‑50 (MQ=255)                atTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTAtctc                                <  2:474255/50‑1 (MQ=255)                 ttAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTcc                               >  2:284370/1‑50 (MQ=255)                  tAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTccc                              >  1:1801353/1‑50 (MQ=255)                    aCCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCt                            >  2:1597191/1‑50 (MQ=255)                      cATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCtt                           >  1:1632853/1‑49 (MQ=39)                       aTTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCttt                          >  1:61551/1‑49 (MQ=39)                       aTTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCttt                          >  2:1644580/1‑49 (MQ=39)                       aTTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTa                         <  1:1080090/50‑1 (MQ=255)                          cccATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAtt                      >  2:761416/1‑50 (MQ=255)                             aTACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAttttt                   <  1:1539275/50‑1 (MQ=255)                                cAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCct                <  2:270575/50‑1 (MQ=255)                                cAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCct                <  1:925619/50‑1 (MQ=255)                                 aaTTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCctc               >  1:1623638/1‑50 (MQ=255)                                 aaTTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCctc               >  1:510758/1‑50 (MQ=255)                                   ttAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTa             >  2:1185606/1‑50 (MQ=255)                                   ttAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTa             >  1:445350/1‑50 (MQ=255)                                   ttAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTa             >  1:528246/1‑50 (MQ=255)                                    tAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTAt            >  2:1242096/1‑50 (MQ=255)                                    tAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTAt            >  2:1993779/1‑50 (MQ=255)                                      gTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtc          >  2:119227/1‑50 (MQ=255)                                      gTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtc          >  1:1867988/1‑50 (MQ=255)                                       ttAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtct         >  1:1035091/1‑50 (MQ=255)                                       ttAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtct         >  1:584755/1‑50 (MQ=255)                                       ttAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATtct         >  1:119336/1‑50 (MQ=255)                                            aaaaag‑GGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGt   >  1:531919/2‑50 (MQ=255)                                             aaaag‑GGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGtt  >  2:1971595/1‑50 (MQ=255)                                             aaaag‑GGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGtt  >  1:918997/1‑50 (MQ=255)                                              aaag‑GGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGtt  >  1:875140/1‑49 (MQ=255)                                                  |                                              GGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAA‑GGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTT  >  NC_017635/536933‑537026

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