Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,031,241	A→G	Y60H (TAT→CAT)	yccA ←	Modulator of FtsH protease, inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,031,241	0	A	G	100.0%	180.1 / NA	51	Y60H (TAT→CAT)	yccA	Modulator of FtsH protease, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/25);  total (26/25)

ACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAG  >  NC_000913/1031124‑1031376                                                                                                                      |                                                                                                                                        acagaTAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAAt                                                                                                                    <  2:228363/139‑1 (MQ=255) acagaTAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAAt                                                                                                                    >  1:228363/1‑139 (MQ=255) acagaTAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAAt                                                                                                                    <  1:19223/139‑1 (MQ=255)     aTAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGa                                                                                                                >  2:268671/1‑139 (MQ=255)         gtgtTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCg                                                                                                            <  1:464165/139‑1 (MQ=255)            ttCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGgag                                                                                                         <  1:1446136/139‑1 (MQ=255)             tCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGgaga                                                                                                        >  2:154741/1‑139 (MQ=255)                 aaTAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGgc                                                                                                    >  1:1229622/1‑139 (MQ=255)                      gTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCAt                                                                                               <  2:904946/139‑1 (MQ=255)                      gTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCAt                                                                                               >  1:904946/1‑139 (MQ=255)                      gTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCAt                                                                                               <  1:1167595/139‑1 (MQ=255)                           aGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                          <  1:700591/139‑1 (MQ=255)                           aGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                          >  2:597884/1‑139 (MQ=255)                               tataACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGt                                                                                      >  2:390436/1‑139 (MQ=255)                               tataACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGt                                                                                      >  2:269994/1‑139 (MQ=255)                                    cccAGAAAACCGGTAAAGGCGAATGCAGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTgg                                                                                 <  2:1002285/139‑1 (MQ=255)                                       aGAAAACCGGTAAAGGCGAATGCGGAGATCATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTc                                                                                                                         <  1:541010/96‑1 (MQ=255)                                       aGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTc                                                                                                                         >  2:541010/1‑96 (MQ=255)                                        gAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGt                                                                             >  2:79974/1‑139 (MQ=255)                                         aaaaCCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGtt                                                                            >  1:313015/1‑139 (MQ=255)                                         aaaaCCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGtt                                                                            >  1:1243392/1‑139 (MQ=255)                                               ggTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTa                                                                      >  1:1405390/1‑139 (MQ=255)                                                     ggCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGcc                                                                >  2:1132013/1‑139 (MQ=255)                                                        gAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGt                                                                                      >  2:201553/1‑114 (MQ=255)                                                        gAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGt                                                                                      <  1:201553/114‑1 (MQ=255)                                                        gAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGaa                                                             >  1:1241098/1‑139 (MQ=255)                                                             cGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGc                                                        <  1:1132013/139‑1 (MQ=255)                                                                      tCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAg                                               >  1:590956/1‑139 (MQ=255)                                                                        ccGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGc                                             <  1:269994/139‑1 (MQ=255)                                                                         cGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCt                                            <  1:1086931/139‑1 (MQ=255)                                                                         cGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCt                                            <  1:191915/139‑1 (MQ=255)                                                                            tCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTcag                                         >  1:376137/1‑139 (MQ=255)                                                                               gCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTcagcag                                      <  2:1243392/139‑1 (MQ=255)                                                                                 ttATTCGCCGTTTTATAGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGaa                                    <  1:268671/139‑1 (MQ=255)                                                                                               ataGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACgc                      <  1:209707/139‑1 (MQ=255)                                                                                                 aGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACgcag                    <  2:8578/139‑1 (MQ=255)                                                                                                 aGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACgcag                    <  2:376137/139‑1 (MQ=255)                                                                                                 aGGTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACgcag                    <  2:313015/139‑1 (MQ=255)                                                                                                   gTCAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACgcagca                  >  1:1250018/1‑139 (MQ=255)                                                                                                     cAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                          <  1:860095/65‑1 (MQ=255)                                                                                                     cAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCa                                                                                          >  2:860095/1‑65 (MQ=255)                                                                                                     cAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTcagca                                       <  2:314013/116‑1 (MQ=255)                                                                                                     cAGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTcagca                                       >  1:314013/1‑116 (MQ=255)                                                                                                      aGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGc                                                        >  1:443638/1‑98 (MQ=255)                                                                                                      aGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGc                                                        <  2:443638/98‑1 (MQ=255)                                                                                                      aGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGc                                                        <  1:1249979/98‑1 (MQ=255)                                                                                                      aGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGc                                                        >  2:1249979/1‑98 (MQ=255)                                                                                                      aGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCt               >  2:1089877/1‑139 (MQ=255)                                                                                                      aGGAACATCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCt               >  2:892046/1‑139 (MQ=255)                                                                                                             tCAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGt        >  2:1257967/1‑139 (MQ=255)                                                                                                              cAAACCATGCATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTg       <  1:40838/139‑1 (MQ=255)                                                                                                                   catgcataCCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTTTGGGTGCTAAg  <  2:1428047/139‑1 (MQ=255)                                                                                                                      |                                                                                                                                        ACAGATAGGTGTTCAGAATAGGTCCGAGGATATAACCCAGAAAACCGGTAAAGGCGAATGCGGAGATAATCCCGGTCGGCTTATTCGCCGTTTTATAGGTCAGGAACATCAAACCATACATACCCACCAGCGTCAGAATCAGACCCGGAGATGGCAGCATCAGCACCGTACTGGCAGTTGCGGTAATCGCCGAAAAGGCCAGCGTCAGGCTCAGCAGAAAATAGGTATTACGCAGCACCTTATGGGTGCTAAG  >  NC_000913/1031124‑1031376

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