Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	524,888	T→C	L543S (TTG→TCG)	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	524,888	0	T	C	100.0%	170.6 / NA	47	L543S (TTG→TCG)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/18);  total (29/18)

GCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACC  >  NC_000913/524754‑525018                                                                                                                                       |                                                                                                                                   gCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTgg                                                                                                                               <  1:819804/140‑1 (MQ=255)           gagaCATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCGCCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGa                                                                                                                     <  1:752020/140‑1 (MQ=255)           gagaCATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGa                                                                                                                     <  2:738361/140‑1 (MQ=255)           gagaCATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGa                                                                                                                     >  1:475126/1‑140 (MQ=255)               cATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGc                                                                                                                 >  2:599674/1‑140 (MQ=255)                aTCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGc                                                                                                                                  >  2:832120/1‑122 (MQ=255)                aTCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCt                                                                                                                >  1:284306/1‑140 (MQ=255)                aTCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCt                                                                                                                <  1:157256/140‑1 (MQ=255)                  cgcgCAGCGGGGAGGCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCt                                                                                                                <  2:769806/138‑1 (MQ=255)                  cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCt                                                                                                                <  2:53325/138‑1 (MQ=255)                  cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCt                                                                                                                >  1:53325/1‑138 (MQ=255)                  cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCt                                                                                                                >  1:769806/1‑138 (MQ=255)                  cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCg                                                                                                              <  1:171888/140‑1 (MQ=255)                   gcgcAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCgg                                                                                                             >  2:664725/1‑140 (MQ=255)                      cAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCgg                                                                                                             >  1:878876/1‑137 (MQ=255)                      cAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCgg                                                                                                             <  2:878876/137‑1 (MQ=255)                        gCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTAc                                                                                                        <  1:599674/140‑1 (MQ=255)                            ggAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGa                                                                                                    >  1:1101736/1‑140 (MQ=255)                            ggAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGa                                                                                                    <  2:290952/140‑1 (MQ=255)                             gagaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGAc                                                                                                   <  1:664725/140‑1 (MQ=255)                             gagaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGAc                                                                                                   <  1:719789/140‑1 (MQ=255)                                  aGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGtt                                                                                              <  1:784901/140‑1 (MQ=255)                                   gTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTa                                                                                             >  1:1146165/1‑140 (MQ=255)                                      cgctaccgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGa                                                                                          >  1:82311/1‑140 (MQ=255)                                          accgctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATAc                                                                                      >  1:438971/1‑140 (MQ=255)                                              ctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGAcc                                                                                  >  1:853898/1‑140 (MQ=255)                                              ctacGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGAcc                                                                                  >  1:959000/1‑140 (MQ=255)                                                   aTGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTc                                                                             >  1:1061363/1‑140 (MQ=255)                                                   aTGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTc                                                                             >  2:38685/1‑140 (MQ=255)                                                     gACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCgg                                                                           >  1:556330/1‑140 (MQ=255)                                                      aCGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGc                                                                          >  2:985952/1‑140 (MQ=255)                                                        gcgcACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGTGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCa                                                                        >  1:691269/1‑140 (MQ=255)                                                           cacaGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGAt                                                                     <  1:165981/140‑1 (MQ=255)                                                           cacaGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGAt                                                                     >  2:463775/1‑140 (MQ=255)                                                             caGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGa                                                                   >  1:662609/1‑140 (MQ=255)                                                                    ttACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGt                                                            <  1:463775/140‑1 (MQ=255)                                                                    ttACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGt                                                            <  2:959000/140‑1 (MQ=255)                                                                      aCCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTcc                                                          >  1:132466/1‑140 (MQ=255)                                                                            gacgacAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACcgcg                                                    >  2:1163303/1‑140 (MQ=255)                                                                                    cGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGc                                            >  2:891071/1‑140 (MQ=255)                                                                                     ggATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCa                                           <  1:905171/140‑1 (MQ=255)                                                                                           gACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGcc                                     >  1:700991/1‑140 (MQ=255)                                                                                                 cagcaCCCGGCAGATGACCTGGAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTAcc                               <  1:38685/140‑1 (MQ=255)                                                                                                    cccccGGCAGATGACCGGGAGCCGCTACGGGCAGTCGTTGGCTTTCACCGACTGCTCGGGCCACCAGCCCCGTTATGAAGACGGCCGCTTCGGCCAGATGTGGGCGGTCCACCGCGAGGAAGGCATCCGCCTTTTccgcc                            <  2:1061363/138‑1 (MQ=255)                                                                                                                   cgggAGCCGCTGCGGGCAGTCGCTGGCGTCCACCGACTGCTGGGGCTACCAGACCCTTTATGAATACGCCCGCTTCGGCCAGATGCCGGGGGTCCCCCGCGAGGAAGGCATCAGCCTTTCCCGCCGCTATGACAACCGTg             <  2:552264/138‑1 (MQ=255)                                                                                                                    gggAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGTCTTTACCGCCGCTATGACAACCGTgg            <  1:985952/139‑1 (MQ=255)                                                                                                                     ggAGCCGCTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGc           >  1:1167827/1‑140 (MQ=255)                                                                                                                            cTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTaa    >  1:552264/1‑140 (MQ=255)                                                                                                                            cTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTaa    >  1:651384/1‑140 (MQ=255)                                                                                                                            cTACGGGCAGTCGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTaa    >  1:2007/1‑140 (MQ=255)                                                                                                                              aCGGGCAGCCGCTGGCTTTCACCGACTGCTCGGGCTCCAAGACCCGTTATGAATACGCCCGCTTCGGCAAGATGACGGCGGTCCCCCGCGGGGAAGGCATCAGCCTTTCCCGCCGCTATGACAACCGTGGCCGGTTAAcc  <  2:284306/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   GCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACC  >  NC_000913/524754‑525018

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: ‑