Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,382,394	T→C	G65G (GGT→GGC)	degS →	serine endoprotease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,382,394	0	T	C	100.0%	194.6 / NA	56	G65G (GGT→GGC)	degS	serine endoprotease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/31);  total (25/31)

CTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGA  >  NC_000913/3382257‑3382528                                                                                                                                          |                                                                                                                                       ctgGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCACGCTGCGACAGCGGTGTTTAATGTTTACAACCGTGGctt                                                                                                                                      >  2:230381/1‑140 (MQ=255)      tGCCATGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAaca                                                                                                                                 >  2:679773/1‑140 (MQ=255)      tGCCATGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAaca                                                                                                                                 >  2:817656/1‑140 (MQ=255)          aTGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCaa                                                                                                                             <  1:817656/140‑1 (MQ=255)          aTGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCaa                                                                                                                             <  1:575952/140‑1 (MQ=255)           tGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCaa                                                                                                                             <  1:777298/139‑1 (MQ=255)           tGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCaa                                                                                                                             >  2:777298/1‑139 (MQ=255)              cTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAActct                                                                                                                         <  1:4380/140‑1 (MQ=255)                   cTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACaa                                                                                                                    <  2:768317/140‑1 (MQ=255)                   cTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACaa                                                                                                                    <  1:55431/140‑1 (MQ=255)                     gcgcAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAAcc                                                                                                                  <  2:633269/140‑1 (MQ=255)                     gcgcAGCCTTAACCCACTTTCCCCTTCGCCATTTGACCGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAAcc                                                                                                                  <  1:230381/140‑1 (MQ=255)                        cAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGc                                                                                                               >  2:43840/1‑140 (MQ=255)                           ccTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTg                                                                                                            >  1:473623/1‑140 (MQ=255)                                aCCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATc                                                                                                       >  2:646222/1‑140 (MQ=255)                                     cTTTCCACTCCGCAGTTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCAc                                                                                                  >  2:543246/1‑140 (MQ=255)                                     cTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACaa                                                                                                                    <  2:474862/122‑1 (MQ=255)                                     cTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACaa                                                                                                                    >  1:474862/1‑122 (MQ=255)                                     cTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCg                                                                                                     <  1:652284/137‑1 (MQ=255)                                     cTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCg                                                                                                     >  2:652284/1‑137 (MQ=255)                                     cTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCCTGAGATCCGCAc                                                                                                  >  2:535100/1‑140 (MQ=255)                                     cTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCACCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACAAACCCTTACAACCTGCTTGAGATCCGCAc                                                                                                  >  1:392691/1‑140 (MQ=255)                                              ccGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATc                                                                                         <  2:728707/140‑1 (MQ=255)                                              ccGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATc                                                                                         <  2:473623/140‑1 (MQ=255)                                              ccGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATc                                                                                         <  1:121154/140‑1 (MQ=255)                                                gCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCg                                                                                       >  2:834846/1‑140 (MQ=255)                                                  aaTTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATccagt                                                                                     <  2:31800/140‑4 (MQ=255)                                                              ccGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGAt                                                                         >  2:759889/1‑140 (MQ=255)                                                                 aTGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCaa                                                                      <  1:679993/140‑1 (MQ=255)                                                                 aTGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCaa                                                                      <  1:725389/140‑1 (MQ=255)                                                                  tGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAAc                                                                     <  1:679773/140‑1 (MQ=255)                                                                   gagaCGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAAcg                                                                    <  1:759889/140‑1 (MQ=255)                                                                   gagaCGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAAcg                                                                    <  1:646222/140‑1 (MQ=255)                                                                   gagaCGCCAGCCAGCTGTAATCCGGTGGTTGCCTTGGCCGCGCCAGCGATGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGTATCAAcg                                                                    <  2:418197/140‑1 (MQ=255)                                                                      aCGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAAcgcg                                                                  >  2:772906/1‑139 (MQ=255)                                                                      aCGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAAcgcg                                                                  <  1:772906/139‑1 (MQ=255)                                                                      aCGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCgg                                                                 >  1:775043/1‑140 (MQ=255)                                                                      aCGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCgg                                                                 >  1:754778/1‑140 (MQ=255)                                                                      aCGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCgg                                                                 >  1:160763/1‑140 (MQ=255)                                                                       cGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGt                                                                <  2:754778/140‑1 (MQ=255)                                                                           tGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACACTGGGATCCGGGGTCATCATGTATCAACGCGGTtata                                                            >  1:421735/1‑140 (MQ=255)                                                                               aGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATatcat                                                        <  1:43840/140‑1 (MQ=255)                                                                               aGCTATAATCTGGCGGTGCGCCACGCCGCGGCACCGGTGCTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATatcat                                                        <  2:421735/140‑1 (MQ=255)                                                                                gCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATatcatc                                                       >  1:136817/1‑140 (MQ=255)                                                                                  tataATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCAc                                                     <  1:706131/140‑1 (MQ=255)                                                                                          ggCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAac                                             <  2:556507/140‑1 (MQ=255)                                                                                           gCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAaca                                            >  2:836319/1‑140 (MQ=255)                                                                                           gCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAaca                                            >  2:571870/1‑140 (MQ=255)                                                                                            cGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAacac                                           >  2:524020/1‑140 (MQ=255)                                                                                              gTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGt                                         >  2:526966/1‑140 (MQ=255)                                                                                                       gccgcgccAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAacga                                <  1:53455/140‑1 (MQ=255)                                                                                                           cgccAGCGGTGGTTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCGTCACCAATAAACAGGTCATCAACGACGcc                            >  1:335879/1‑140 (MQ=255)                                                                                                              cAGCGGTGGTTAACGTTAACAACCGTGGCTTGAACGCCACCTCTCACATCCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCgat                         <  2:335879/140‑1 (MQ=255)                                                                                                                      gtTAACGTTTACAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGatcat                 <  1:836319/140‑1 (MQ=255)                                                                                                                           cGTTTACAACCGTGGCTAGAACCCCAACTCTCAAAACCAACTTGAGATGCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACAACGCCGAGCAGATCAtcgtcg            <  2:392691/140‑1 (MQ=255)                                                                                                                                aCAACCGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGCGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACCCGTCATCCACGTCGACGATCATATCATCGTCGCCTTa       >  1:409284/1‑140 (MQ=255)                                                                                                                                     cGTGGCTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCTGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGa  <  1:535100/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       CTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGTGGTTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGA  >  NC_000913/3382257‑3382528

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