Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,650,173	T→C	K705K (AAA→AAG)	yfhM ←	alpha2‑macroglobulin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,650,173	0	T	C	100.0%	167.0 / NA	47	K705K (AAA→AAG)	yfhM	alpha2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/23);  total (24/23)

GCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCTTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCT  >  NC_000913/2650038‑2650308                                                                                                                                        |                                                                                                                                        gcTACCGTCGCTGTTGTCCTGCCAGCTATAGCCATCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCa                                                                                                                                     <  2:458561/140‑1 (MQ=255)  cTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAg                                                                                                                                    >  2:264978/1‑140 (MQ=255)      cGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATcc                                                                                                                                <  1:264978/140‑1 (MQ=255)             gttgtCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGttt                                                                                                                         >  2:577239/1‑140 (MQ=255)                gtCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGtt                                                                                                                      <  1:341804/140‑1 (MQ=255)                           tataGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATc                                                                                                           >  2:616733/1‑140 (MQ=255)                                  agccCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCtt                                                                                                    <  2:707384/140‑1 (MQ=255)                                    cccAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCtttt                                                                                                  >  2:735636/1‑140 (MQ=255)                                    cccAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCtttt                                                                                                  >  2:403976/1‑140 (MQ=255)                                              aaCACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACt                                                                                        >  2:602754/1‑140 (MQ=255)                                              aaCACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACt                                                                                        >  2:719370/1‑140 (MQ=255)                                              aaCACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACt                                                                                        >  1:501592/1‑140 (MQ=255)                                              aaCACTACTGACCGCTTCATTCGGCGCTTTGACTACCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACt                                                                                        >  2:50868/1‑140 (MQ=255)                                                   tactGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGAc                                                                                   >  2:911471/1‑140 (MQ=255)                                                      tGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGc                                                                                <  1:719370/140‑1 (MQ=255)                                                       gACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGcc                                                                               <  2:1013562/140‑1 (MQ=255)                                                        aCCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCa                                                                              >  1:533384/1‑140 (MQ=255)                                                        aCCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCa                                                                              <  1:911471/140‑1 (MQ=255)                                                         ccGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAg                                                                             <  1:616733/140‑1 (MQ=255)                                                          cGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGAt                                                                                              <  2:613982/122‑1 (MQ=255)                                                          cGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGAt                                                                                              >  1:613982/1‑122 (MQ=255)                                                             ttCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACtgt                                                                                      <  2:986105/127‑1 (MQ=255)                                                             ttCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACtgt                                                                                      >  1:986105/1‑127 (MQ=255)                                                             ttCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGc                                                                                <  1:601122/133‑1 (MQ=255)                                                             ttCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGc                                                                                >  2:601122/1‑133 (MQ=255)                                                             ttCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCtt                                                                         >  2:203435/1‑140 (MQ=255)                                                             ttCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCtt                                                                         >  1:497886/1‑140 (MQ=255)                                                             ttCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTATGAGTGCCAGCCtt                                                                         >  2:626868/1‑140 (MQ=255)                                                              tCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTc                                                                        <  1:841619/140‑1 (MQ=255)                                                              tCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTc                                                                        <  1:403976/140‑1 (MQ=255)                                                                    ggCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCAtcttc                                                                  >  2:914914/1‑140 (MQ=255)                                                                       gcTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGa                                                               <  1:203435/140‑1 (MQ=255)                                                                        cTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGAc                                                              <  1:602754/140‑1 (MQ=255)                                                                         tttGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGAcc                                                             <  1:735636/140‑1 (MQ=255)                                                                         tttGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGAcc                                                             <  2:533384/140‑1 (MQ=255)                                                                         tttGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGAcc                                                             <  1:914914/140‑1 (MQ=255)                                                                         tttGACTACCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGAcc                                                             <  1:50868/140‑1 (MQ=255)                                                                            gACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGt                                                          <  2:501592/140‑1 (MQ=255)                                                                                                 ccACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGtt                                     >  2:541359/1‑140 (MQ=255)                                                                                                 ccACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGtt                                     >  2:611667/1‑140 (MQ=255)                                                                                                         ccGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATc                             <  1:577239/140‑1 (MQ=255)                                                                                                          cGGAAAACTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCa                            >  2:745513/1‑140 (MQ=255)                                                                                                               aaCTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGgcgc                       >  1:542916/1‑140 (MQ=255)                                                                                                                 cTGACCTTGCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCAc                     >  2:203266/1‑140 (MQ=255)                                                                                                                         gCCGGTTTCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGc             >  1:788657/1‑140 (MQ=255)                                                                                                                               ttCGTCCGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGaca       <  1:745513/140‑1 (MQ=255)                                                                                                                                    ccGCCTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCt  <  2:497886/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        GCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCTTTCAGATCCAGAGTTTGTTCATTTTCGATCAGATCTTTTTGATCAAACTGTGACTGCCAGCCTTCATCTTCTGACCAGTTCCAGTAGTAATCGCGGCGTTCGCGAATCAGGCGCACCTGCAAGCCCGACACGGCT  >  NC_000913/2650038‑2650308

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