Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,043,500	A→G	Y28C (TAC→TGC)	dsbA →	thiol:disulfide oxidoreductase DsbA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,043,500	0	A	G	100.0%	164.4 / NA	48	Y28C (TAC→TGC)	dsbA	thiol:disulfide oxidoreductase DsbA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/21);  total (27/21)

TACAAGAACCCCCTTTGCAATTAACACCTATGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAA  >  NC_000913/4043364‑4043630                                                                                                                                         |                                                                                                                                   tACAAGAACCCCCTTTGCAATTAACACCTATGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGATTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTgcac                                                                                                                                 >  2:205728/1‑140 (MQ=255)        aCCCCCTTTGCAATTAACACCTATGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTg                                                                                                                          >  1:112933/1‑140 (MQ=255)         cccccTTTGCAATTAACACCTATGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTgg                                                                                                                         >  2:77674/1‑140 (MQ=255)                tGCAATTAACACCTATGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAAcc                                                                                                                  >  2:407168/1‑140 (MQ=255)                             tATGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGgcgcg                                                                                                     >  1:764273/1‑140 (MQ=255)                               tGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGcgcc                                                                                                   >  1:746123/1‑140 (MQ=255)                               tGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGcgcc                                                                                                   >  1:518764/1‑140 (MQ=255)                                     aaTCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGt                                                                                             >  1:27039/1‑140 (MQ=255)                                      aTCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCaa                                                                                               >  1:342033/1‑137 (MQ=255)                                      aTCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCaa                                                                                               <  2:342033/137‑1 (MQ=255)                                      aTCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTg                                                                                            >  2:776108/1‑140 (MQ=255)                                      aTCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTg                                                                                            >  1:661905/1‑140 (MQ=255)                                         ggAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTg                                                                                         <  2:746123/140‑1 (MQ=255)                                         ggAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTg                                                                                         <  2:63543/140‑1 (MQ=255)                                           agagagTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTggg                                                                                       >  2:404404/1‑139 (MQ=255)                                           agagagTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTg                                                                                         <  2:37423/138‑1 (MQ=255)                                           agagagTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTg                                                                                         >  1:37423/1‑138 (MQ=255)                                           agagagTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGa                                                                                       >  2:461707/1‑140 (MQ=255)                                           agagagTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGa                                                                                       >  2:298093/1‑140 (MQ=255)                                             agagTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGt                                                                                     >  1:203340/1‑140 (MQ=255)                                                gTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGtttt                                                                                  >  1:710834/1‑140 (MQ=255)                                                 tAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGttttt                                                                                 >  2:142767/1‑140 (MQ=255)                                                    aTCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTtctc                                                                              >  1:391016/1‑140 (MQ=255)                                                         gAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCtctttct                                                                         <  2:391016/140‑1 (MQ=255)                                                         gAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCtctttct                                                                         >  2:597690/1‑140 (MQ=255)                                                         gAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCtctttct                                                                         >  2:531622/1‑140 (MQ=255)                                                         gAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCtctttct                                                                         <  1:531622/140‑1 (MQ=255)                                                         gAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCtctttct                                                                         >  1:503729/1‑140 (MQ=255)                                                              aGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGc                                                                    >  1:94291/1‑140 (MQ=255)                                                                aTTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGccc                                                                  <  1:77674/140‑1 (MQ=255)                                                                 tttGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCg                                                                 <  2:764273/140‑1 (MQ=255)                                                                  ttGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGc                                                                <  2:148620/140‑1 (MQ=255)                                                                  ttGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGc                                                                <  1:776108/140‑1 (MQ=255)                                                                  ttGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGc                                                                <  2:518764/140‑1 (MQ=255)                                                                   tggctggcGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGcg                                                               <  2:203340/140‑2 (MQ=255)                                                                   tggctggcGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCa                                                               <  1:461707/140‑1 (MQ=255)                                                                    ggctggcGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCAc                                                              >  1:789616/1‑140 (MQ=255)                                                                          cgctggcTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTAt                                                        <  1:407168/140‑1 (MQ=255)                                                                          cgctgccTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTAt                                                        <  2:710834/140‑1 (MQ=255)                                                                               gctggTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGtt                                                   <  1:340973/140‑1 (MQ=255)                                                                                    tttAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTgaag                                              <  1:142767/140‑1 (MQ=255)                                                                                    tttAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTgaag                                              <  2:112933/140‑1 (MQ=255)                                                                                              gcgtttagcgCATCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCa                                    <  1:298093/140‑1 (MQ=255)                                                                                                          tCGGCGGCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATaa                        >  1:529956/1‑140 (MQ=255)                                                                                                           cggcggCGCAGTATGAAGATGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAAt                       <  1:597690/140‑1 (MQ=255)                                                                                                                              tGGTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCgg    <  2:94291/140‑1 (MQ=255)                                                                                                                                gTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGaa  >  2:561909/1‑140 (MQ=255)                                                                                                                                gTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGaa  >  1:643587/1‑140 (MQ=255)                                                                                                                                gTAAACAGTGCACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGaa  >  1:301701/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                   TACAAGAACCCCCTTTGCAATTAACACCTATGTATTAATCGGAGAGAGTAGATCATGAAAAAGATTTGGCTGGCGCTGGCTGGTTTAGTTTTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAA  >  NC_000913/4043364‑4043630

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