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%	190.4 / NA	50	I240V (ATC→GTC)	gspD	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/27);  total (23/27)

ATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGA  >  NC_000913/3456972‑3457224                                                                                                                           |                                                                                                                                   atGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGa                                                                                                                    <  1:782344/139‑1 (MQ=255)        aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                             >  2:449098/1‑139 (MQ=255)              aGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                       <  1:196468/139‑1 (MQ=255)                aTCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAgcg                                                                                                     <  2:300166/139‑1 (MQ=255)                   tcGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                  <  1:662552/139‑1 (MQ=255)                     gCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCa                                                                                                <  2:370422/139‑1 (MQ=255)                       cGAGATTCTTAATCAATCAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACt                                                                                              >  1:835718/1‑139 (MQ=255)                          gaTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAACCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCa                                                                                           >  2:659848/1‑139 (MQ=255)                           aTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                  <  1:127316/131‑1 (MQ=255)                           aTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                  >  2:127316/1‑131 (MQ=255)                            ttCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCAtt                                                                                         >  1:537674/1‑139 (MQ=255)                                   tCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGaaa                                                                                  <  2:953655/139‑1 (MQ=255)                                       ttAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGcc                                                                              >  2:955235/1‑139 (MQ=255)                                               cGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                  >  1:373076/1‑111 (MQ=255)                                               cGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                  <  2:373076/111‑1 (MQ=255)                                               cGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTc                                                                      >  2:334565/1‑139 (MQ=255)                                                 aaaGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGa                                                                    <  1:955235/139‑1 (MQ=255)                                                      cACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAagaga                                                               <  1:104917/139‑1 (MQ=255)                                                         ggTAAAAGCCAGATGCCAGCCGTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                            <  1:659848/139‑1 (MQ=255)                                                         ggTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                            >  1:940549/1‑139 (MQ=255)                                                         ggTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                            <  1:850259/139‑1 (MQ=255)                                                         ggTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                            <  1:353212/139‑1 (MQ=255)                                                                   aGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGaaa                                                  >  2:579457/1‑139 (MQ=255)                                                                       gCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAcc                                              <  1:579457/139‑1 (MQ=255)                                                                          aGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCgg                                           >  2:183712/1‑139 (MQ=255)                                                                              ctcctcTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGttt                                       <  1:319297/139‑1 (MQ=255)                                                                                cctcTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTttat                                     <  1:334565/139‑1 (MQ=255)                                                                                 ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGaa                                                                                   <  1:1065564/92‑1 (MQ=255)                                                                                 ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGaa                                                                                   >  2:1065564/1‑92 (MQ=255)                                                                                    tcCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTAcc                                 >  2:788558/1‑139 (MQ=255)                                                                                      cgcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                             <  2:430937/93‑1 (MQ=255)                                                                                      cgcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                             >  1:430937/1‑93 (MQ=255)                                                                                      cgcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTg                               >  2:501710/1‑139 (MQ=255)                                                                                        cgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGaa                             <  1:183712/139‑1 (MQ=255)                                                                                                gtgGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                            >  2:205830/1‑100 (MQ=255)                                                                                                gtgGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                            <  1:205830/100‑1 (MQ=255)                                                                                                      gATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCa               <  1:979146/139‑1 (MQ=255)                                                                                                        tAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACg             <  1:788558/139‑1 (MQ=255)                                                                                                        tAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGAACGCCAGCGCATCACTTCATTGCTAAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTGATTACCTGAAATATGCTAAAGCCACg             >  2:89016/1‑139 (MQ=255)                                                                                                           gCGAACCCACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATAACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTAAAATATGCTAAAGCCACGAAt          <  2:695557/139‑1 (MQ=255)                                                                                                           gCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                             <  1:15920/72‑1 (MQ=255)                                                                                                           gCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                             >  2:15920/1‑72 (MQ=255)                                                                                                           gCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAAt          >  1:16454/1‑139 (MQ=255)                                                                                                           gCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAAt          >  1:373944/1‑139 (MQ=255)                                                                                                                ccAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAAt          <  1:524265/134‑1 (MQ=255)                                                                                                                ccAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAAt          >  2:524265/1‑134 (MQ=255)                                                                                                                ccAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCtggt     >  2:620018/1‑139 (MQ=255)                                                                                                                 cAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCtggtg    <  1:620018/139‑1 (MQ=255)                                                                                                                  aaCTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCtggtgg   <  1:501710/139‑1 (MQ=255)                                                                                                                   aCTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGa  >  1:382553/1‑139 (MQ=255)                                                                                                                           |                                                                                                                                   ATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGA  >  NC_000913/3456972‑3457224

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