Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,643,336	Δ1 bp	coding (1847/2043 nt)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,334	0	C	.	100.0%	202.1 / NA	47	coding (1849/2043 nt)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (19/28);  total (19/28)

CGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGATGGCACCACGGCAACCAGTTTCTTGATTTCTGCCAGAGTTTCGAGGA  >  NC_000913/3643195‑3643465                                                                                                                                            |                                                                                                                                    cGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGc                                                                                                                                     >  2:897227/1‑140 (MQ=255) cGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGc                                                                                                                                     >  1:404399/1‑140 (MQ=255) cGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGc                                                                                                                                     >  1:949553/1‑140 (MQ=255) cGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGc                                                                                                                                     <  1:76243/140‑1 (MQ=255)  ggAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGcc                                                                                                                                    <  2:840550/140‑1 (MQ=255)     aGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGccaca                                                                                                                                 <  2:404399/140‑4 (MQ=255)     aGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGccaca                                                                                                                                 <  2:561650/140‑4 (MQ=255)        gTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGt                                                                                                                             <  1:160888/140‑1 (MQ=255)        gTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGt                                                                                                                             <  2:519374/140‑1 (MQ=255)         tttGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTa                                                                                                                            >  2:844123/1‑140 (MQ=255)          ttGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAg                                                                                                                           <  1:220874/140‑1 (MQ=255)          ttGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAg                                                                                                                           <  2:92808/140‑1 (MQ=255)          ttGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAg                                                                                                                           <  1:897227/140‑1 (MQ=255)               cAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTgta                                                                                                                      <  1:844123/140‑1 (MQ=255)                    ccATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAAc                                                                                                                 >  1:823525/1‑140 (MQ=255)                      aTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCt                                                                                                               <  2:949553/140‑1 (MQ=255)                      aTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCt                                                                                                               <  2:847866/140‑1 (MQ=255)                       tCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTg                                                                                                              <  2:242956/140‑1 (MQ=255)                            tcttctGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCa                                                                                                         >  2:94089/1‑140 (MQ=255)                              ttctGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATa                                                                                                       >  2:412994/1‑140 (MQ=255)                                  gAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAacca                                                                                                   <  2:601533/140‑1 (MQ=255)                                   aaCCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTg                                                                                                              >  2:791532/1‑128 (MQ=255)                                   aaCCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTg                                                                                                              <  1:791532/128‑1 (MQ=255)                                           cGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGaaa                                                                                          <  2:896819/140‑1 (MQ=255)                                             gCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAat                                                                                        <  2:823525/140‑1 (MQ=255)                                             gCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAat                                                                                        <  1:94089/140‑1 (MQ=255)                                             gCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAat                                                                                        <  1:411769/140‑1 (MQ=255)                                                        tgtCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGaa                                                                               <  1:141169/138‑1 (MQ=255)                                                        tgtCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGaa                                                                               >  2:141169/1‑138 (MQ=255)                                                             aGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGc                                                                        >  1:830267/1‑140 (MQ=255)                                                                 aCGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGa                                                                    >  1:84173/1‑140 (MQ=255)                                                                 aCGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGa                                                                    >  1:696794/1‑140 (MQ=255)                                                                   gACTGCCCGGTTTCACGGTTGAAAATGCCCTCCTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGaaa                                                                  <  2:46516/140‑1 (MQ=255)                                                                    aCTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAAc                                                                 >  1:137272/1‑140 (MQ=255)                                                                    aCTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAAc                                                                 <  2:893169/140‑1 (MQ=255)                                                                        cccGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGaaa                                                                  <  1:630189/135‑1 (MQ=255)                                                                        cccGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGaaa                                                                  >  2:630189/1‑135 (MQ=255)                                                                        cccGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGcc                                                             >  1:671963/1‑140 (MQ=255)                                                                        cccGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGcc                                                             >  2:490136/1‑140 (MQ=255)                                                                         ccGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGccc                                                            <  2:137272/140‑1 (MQ=255)                                                                         ccGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGccc                                                            <  1:227228/140‑1 (MQ=255)                                                                             tttCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAgg                                                        <  2:830267/140‑1 (MQ=255)                                                                                  cGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGga                                                   <  2:671963/140‑1 (MQ=255)                                                                                    gTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGgaga                                                 <  2:49771/140‑1 (MQ=255)                                                                                     ttGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGAt                                                <  2:696794/140‑1 (MQ=255)                                                                                       gAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGcccc                                                           >  2:29874/1‑127 (MQ=255)                                                                                       gAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGcccc                                                           <  1:29874/127‑1 (MQ=255)                                                                                                             cgcgAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGTAGATGGCACCACGGCAACCAGTTTCTTg                        >  1:282730/1‑140 (MQ=255)                                                                                                               cgAGAAAGCATCTGCCGCCAGTGCGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGTGGGAAACGGCCCCAGGACGGAGATGGCCCCACGGCAGCCAGTTTCTTGAt                      <  2:565234/140‑1 (MQ=255)                                                                                                               cgAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGATGGCACCACGGCAACCAGTTTCTTGAt                      <  2:84173/140‑1 (MQ=255)                                                                                                               cgAGAAAGCATCTGCCGCCAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGATGGCACCACGGCAACCAGTTTCTTGAt                      >  1:857692/1‑140 (MQ=255)                                                                                                                             ccgccAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGATGGCACCACGGCAACCAGTTTCTTGATTTCTGCCAGAGttt        >  1:738713/1‑140 (MQ=255)                                                                                                                                ccAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGATGGCACCACGGCAACCAGTTTCTTGATTTCTGCCAGAGTTTCGa     >  2:740278/1‑140 (MQ=255)                                                                                                                                 cAGTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGATGGCACCACGGCAACCAGTTTCTTGATTTCTGCCAGAGTTTCGAg    >  1:88441/1‑140 (MQ=255)                                                                                                                                   gTACGTCAG‑CCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGATGGCACCACGGCAACCAGTTTCTTGATTTCTGCCAGAGTTTCGAGGa  <  2:857692/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                    CGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAAGCGTGCGGGAAACGGCCCCAGGACGGAGATGGCACCACGGCAACCAGTTTCTTGATTTCTGCCAGAGTTTCGAGGA  >  NC_000913/3643195‑3643465

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