Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,229,967	G→A	P59L (CCG→CTG)	pepE ←	peptidase E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,229,967	0	G	A	100.0%	215.3 / NA	60	P59L (CCG→CTG)	pepE	peptidase E
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/38);  total (22/38)

GTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGG  >  NC_000913/4229830‑4230099                                                                                                                                          |                                                                                                                                     gTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGa                                                                                                                                    <  1:353324/140‑1 (MQ=255)   gCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGc                                                                                                                                  >  2:518852/1‑140 (MQ=255)    ccAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCg                                                                                                                                 >  1:508868/1‑140 (MQ=255)    ccAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCg                                                                                                                                 <  1:1125752/140‑1 (MQ=255)       gcagcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCgaga                                                                                                                              >  1:366508/1‑140 (MQ=255)           cccGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCg                                                                                                                          >  2:1043647/1‑140 (MQ=255)            ccGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGc                                                                                                                         >  2:645539/1‑140 (MQ=255)             cgcgCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCa                                                                                                                        <  1:1143730/140‑1 (MQ=255)             cgcgCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCa                                                                                                                        <  2:667055/140‑1 (MQ=255)                   cgcgGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGtt                                                                                                                  >  1:862210/1‑140 (MQ=255)                    gcgGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGttt                                                                                                                 <  1:230861/140‑1 (MQ=255)                       gCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCg                                                                                                              <  1:131681/140‑1 (MQ=255)                       gCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCg                                                                                                              <  1:1043647/140‑1 (MQ=255)                       gCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCg                                                                                                              <  2:366508/140‑1 (MQ=255)                       gCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCg                                                                                                              <  2:547776/140‑1 (MQ=255)                       gCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCg                                                                                                              <  1:357321/140‑1 (MQ=255)                          cTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCtg                                                                                                           <  2:561946/140‑1 (MQ=255)                                cAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAat                                                                                                     <  1:518852/140‑1 (MQ=255)                                          aaaGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGt                                                                                           <  1:26163/140‑1 (MQ=255)                                           aaGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTc                                                                                          >  1:499856/1‑140 (MQ=255)                                                     gccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGc                                                                                >  2:878172/1‑140 (MQ=255)                                                           gACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGa                                                                          <  1:710612/140‑1 (MQ=255)                                                                tCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAggg                                                                     <  1:645539/140‑1 (MQ=255)                                                                    gATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATaa                                                                 >  1:1164929/1‑140 (MQ=255)                                                                       ctcAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAaca                                                              <  1:158977/140‑1 (MQ=255)                                                                        tcAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAacac                                                             >  2:94951/1‑140 (MQ=255)                                                                                ttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgc                                                     >  1:803433/1‑140 (MQ=255)                                                                                ttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgc                                                     >  2:1102400/1‑140 (MQ=255)                                                                                 tCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcg                                                    >  1:369063/1‑140 (MQ=255)                                                                                 tCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcg                                                    >  1:289650/1‑140 (MQ=255)                                                                                   aaTCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGATAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcggc                                                  <  2:1107929/140‑1 (MQ=255)                                                                                   aaTCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcggc                                                  <  1:1102400/140‑1 (MQ=255)                                                                                   aaTCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcggc                                                  <  2:803433/140‑1 (MQ=255)                                                                                    aTCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAgcggcg                                                 >  2:540009/1‑140 (MQ=255)                                                                                           cAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGc                                          <  1:682305/140‑1 (MQ=255)                                                                                             aCGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCaa                                        <  1:620227/140‑1 (MQ=255)                                                                                             aCGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCaa                                        <  1:205604/140‑1 (MQ=255)                                                                                             aCGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCaa                                        <  1:456499/140‑1 (MQ=255)                                                                                             aCGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCaa                                        <  1:466170/140‑1 (MQ=255)                                                                                                          cGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAtt                           >  1:386109/1‑140 (MQ=255)                                                                                                          cGCTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAtt                           >  2:444405/1‑140 (MQ=255)                                                                                                            cTATGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAg                         <  2:508868/140‑1 (MQ=255)                                                                                                               tGAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCgg                      >  1:330835/1‑140 (MQ=255)                                                                                                                gAATACCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGc                     <  2:1164929/140‑1 (MQ=255)                                                                                                                    aCCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTg                 <  1:954998/140‑1 (MQ=255)                                                                                                                    aCCGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTg                 >  1:1083144/1‑140 (MQ=255)                                                                                                                     ccGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGc                <  1:597886/140‑1 (MQ=255)                                                                                                                     ccGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGc                <  2:369063/140‑1 (MQ=255)                                                                                                                     ccGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGc                <  2:492479/140‑1 (MQ=255)                                                                                                                      cGGTGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCa               >  1:146447/1‑140 (MQ=255)                                                                                                                         tGACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGt            >  1:1144211/1‑140 (MQ=255)                                                                                                                          gACAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGtt           <  1:444405/140‑1 (MQ=255)                                                                                                                           aCAGAAACACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTc          >  1:1106672/1‑140 (MQ=255)                                                                                                                               aaaCACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTccagc      <  1:94951/140‑1 (MQ=255)                                                                                                                               aaaCACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTccagc      <  1:540009/140‑1 (MQ=255)                                                                                                                               aaaCACCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTccagc      <  1:878172/140‑1 (MQ=255)                                                                                                                                  cacCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTccagc      <  1:39380/137‑1 (MQ=255)                                                                                                                                  cacCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTccagc      >  2:39380/1‑137 (MQ=255)                                                                                                                                   acCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAgg  <  2:146447/140‑1 (MQ=255)                                                                                                                                   acCCAGCAGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAgg  <  2:1083144/140‑1 (MQ=255)                                                                                                                                   acCCAGCAGACCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCAGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAgg  <  2:424212/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     GTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAGCGGCAGTGCATGTTCCAGCCAGG  >  NC_000913/4229830‑4230099

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