Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,457,094	A→G	I240V (ATC→GTC)	gspD →	general secretory pathway component, cryptic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,457,094	0	A	G	100.0%	158.8 / NA	42	I240V (ATC→GTC)	gspD	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/23);  total (19/23)

TTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTG  >  NC_000913/3456960‑3457228                                                                                                                                       |                                                                                                                                       ttCAGCTGGAATATGCGTCAGCGGAAGATCTCGCCGCGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAGAGCCAGATGCCATCCCTCCTCTCCGCGAAGTTTGTGGCGGATAAGCGAACCAACTCTCTTgtcat                                                                                                                                    >  1:218440/1‑139 (MQ=255)      cTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTg                                                                                                                               >  1:857096/1‑139 (MQ=255)             atGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGa                                                                                                                        <  1:891218/139‑1 (MQ=255)              tGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGaa                                                                                                                       >  2:871852/1‑139 (MQ=255)                tgTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGaaaa                                                                                                                     >  1:978485/2‑139 (MQ=255)                cGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGaaaa                                                                                                                     >  2:301460/1‑139 (MQ=255)                   cAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGc                                                                                                                  >  2:1047297/1‑139 (MQ=255)                    aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                                 >  1:290615/1‑139 (MQ=255)                          aGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                           <  1:301460/139‑1 (MQ=255)                             tctcGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAgcgc                                                                                                        <  1:559241/139‑1 (MQ=255)                             tctcGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAgcgc                                                                                                        <  2:350776/139‑1 (MQ=255)                                 gCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCa                                                                                                    <  2:471150/139‑1 (MQ=255)                                   cGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACt                                                                                                  >  2:105152/1‑139 (MQ=255)                                      gaTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCa                                                                                               <  1:267643/139‑1 (MQ=255)                                      gaTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCa                                                                                               <  2:916008/139‑1 (MQ=255)                                              aTCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGaa                                                                                       <  2:897455/139‑1 (MQ=255)                                                   ttatccAGCGAAGCCCACGTTAAAAGCCAAATGCCAG‑CCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                                 <  2:587613/134‑1 (MQ=255)                                                     aaTCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACg                                                                                                               <  2:310251/108‑1 (MQ=255)                                                     aaTCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACg                                                                                                               >  1:310251/1‑108 (MQ=255)                                                      aTCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                               <  1:726284/139‑1 (MQ=255)                                                       tCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGa                                                                              >  2:666911/1‑139 (MQ=255)                                                             aaaGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGa                                                                        >  2:47591/1‑139 (MQ=255)                                                                     ggTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                                >  2:533562/1‑139 (MQ=255)                                                                     ggTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                                <  2:768755/139‑1 (MQ=255)                                                                        aaaaGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagg                                                             <  2:790808/139‑1 (MQ=255)                                                                        aaaaGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagg                                                             <  1:137685/139‑1 (MQ=255)                                                                                 aTGCCAGCCCTCCTCTCCGCGGAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATa                                                    >  1:369531/1‑139 (MQ=255)                                                                                 aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                      >  2:100159/1‑89 (MQ=255)                                                                                 aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                      <  1:100159/89‑1 (MQ=255)                                                                                  tGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAc                                                   <  1:317207/139‑1 (MQ=255)                                                                                          ctcctcTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGttt                                           <  1:47591/139‑1 (MQ=255)                                                                                               ctcCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTAc                                      >  2:156346/1‑139 (MQ=255)                                                                                               ctcCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTAc                                      <  1:533562/139‑1 (MQ=255)                                                                                                       aGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                                 >  1:443763/1‑88 (MQ=255)                                                                                                       aGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                                 <  2:443763/88‑1 (MQ=255)                                                                                                           tgtgGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagga                                                            <  2:201939/105‑1 (MQ=255)                                                                                                           tgtgGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagga                                                            >  1:201939/1‑105 (MQ=255)                                                                                                                   aTAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAcc                                                  >  1:1076496/1‑107 (MQ=255)                                                                                                                   aTAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAcc                                                  <  2:1076496/107‑1 (MQ=255)                                                                                                                   aTAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCAc                  >  2:295754/1‑139 (MQ=255)                                                                                                                    tAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACg                 <  2:972656/139‑1 (MQ=255)                                                                                                                       gCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATATCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAAt              >  1:709111/1‑139 (MQ=255)                                                                                                                                   tctTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTg  <  1:259901/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       TTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTG  >  NC_000913/3456960‑3457228

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