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

ATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGT  >  NC_000913/3456963‑3457221                                                                                                                                    |                                                                                                                                aTCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAg                                                                                                                          >  1:47361/1‑139 (MQ=255)     ggAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGa                                                                                                                      <  2:715100/139‑1 (MQ=255)          atGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGa                                                                                                                 <  2:433950/139‑1 (MQ=255)           tGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGt                                                                                                                         <  1:426463/130‑1 (MQ=255)           tGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGt                                                                                                                         >  2:426463/1‑130 (MQ=255)              gTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGaaaaa                                                                                                             >  1:765730/1‑139 (MQ=255)                cAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGc                                                                                                           >  2:1283711/1‑139 (MQ=255)                 aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                          >  2:800023/1‑139 (MQ=255)                 aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                          >  1:1186487/1‑139 (MQ=255)                   cGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACg                                                                                                        >  1:113865/1‑139 (MQ=255)                         aTCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAgcg                                                                                                  <  2:115545/139‑1 (MQ=255)                            tcGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                               <  2:633370/139‑1 (MQ=255)                             cGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATc                                                                                              >  1:923097/1‑139 (MQ=255)                              gCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCa                                                                                             <  1:946011/139‑1 (MQ=255)                                            tCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGCCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGaaa                                                                               <  1:700885/139‑1 (MQ=255)                                              aaTTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                    >  2:701917/1‑116 (MQ=255)                                              aaTTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                    <  1:701917/116‑1 (MQ=255)                                                ttAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGcc                                                                           >  2:188622/1‑139 (MQ=255)                                                 tAATCAGCGAAAGCCACGGTAAATGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                          <  2:113865/139‑1 (MQ=255)                                                 tAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                          <  2:765730/139‑1 (MQ=255)                                                 tAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                          <  2:1148933/139‑1 (MQ=255)                                                 tAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                          <  1:1283711/139‑1 (MQ=255)                                                 tAATCAACGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                          <  1:449561/139‑1 (MQ=255)                                                  aaTCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCtt                                                                         >  2:351256/1‑139 (MQ=255)                                                   aTCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                        <  2:47361/139‑1 (MQ=255)                                                        cGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                               >  2:836811/1‑111 (MQ=255)                                                        cGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                               <  1:836811/111‑1 (MQ=255)                                                          aaaGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                    >  1:409422/1‑104 (MQ=255)                                                          aaaGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                    <  2:409422/104‑1 (MQ=255)                                                          aaaGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGa                                                                 >  1:1235078/1‑139 (MQ=255)                                                          aaaGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGCGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGa                                                                 >  1:129137/1‑139 (MQ=255)                                                             gCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                        <  2:433099/129‑1 (MQ=255)                                                             gCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                        >  1:433099/1‑129 (MQ=255)                                                                         gCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGa                                                                                                                      <  2:1018566/71‑1 (MQ=255)                                                                         gCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGa                                                                                                                      >  1:1018566/1‑71 (MQ=255)                                                                            aGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGaaa                                               >  1:857964/1‑139 (MQ=255)                                                                             gATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                         <  1:470801/128‑1 (MQ=255)                                                                             gATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                         >  2:470801/1‑128 (MQ=255)                                                                              aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                    >  1:1218929/1‑84 (MQ=255)                                                                              aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                    <  2:1218929/84‑1 (MQ=255)                                                                              aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                               <  2:858911/89‑1 (MQ=255)                                                                              aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                               >  1:858911/1‑89 (MQ=255)                                                                               tGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAc                                            <  1:393069/139‑1 (MQ=255)                                                                                gCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAcc                                           <  2:197091/139‑1 (MQ=255)                                                                                gCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAcc                                           <  1:351256/139‑1 (MQ=255)                                                                                        tcctcTCCGCGAAGATTGTGGCGGATAAGCGGACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTtta                                   <  2:86136/139‑1 (MQ=255)                                                                                        tcctcTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTtta                                   <  2:857964/139‑1 (MQ=255)                                                                                        tcctcTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTtta                                   <  2:923097/139‑1 (MQ=255)                                                                                        tcctcTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTtta                                   <  1:800023/139‑1 (MQ=255)                                                                                          ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCAc                                                                                            <  1:267887/80‑1 (MQ=255)                                                                                          ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCAc                                                                                            >  2:267887/1‑80 (MQ=255)                                                                                          ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                        >  2:237868/1‑100 (MQ=255)                                                                                          ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                        <  1:237868/100‑1 (MQ=255)                                                                                          ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGaa                                                                >  1:135923/1‑108 (MQ=255)                                                                                          ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGaa                                                                <  2:135923/108‑1 (MQ=255)                                                                                          ctctCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTttatt                                 >  2:229813/1‑139 (MQ=255)                                                                                            ctcCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagga                                                     <  2:574677/117‑1 (MQ=255)                                                                                            ctcCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagga                                                     >  1:574677/1‑117 (MQ=255)                                                                                            ctcCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTAc                               <  1:229813/139‑1 (MQ=255)                                                                                               cgcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTg                            >  2:219350/1‑139 (MQ=255)                                                                                                   aaGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAagagag                                                           >  1:1163713/1‑104 (MQ=255)                                                                                                   aaGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAagagag                                                           <  2:1163713/104‑1 (MQ=255)                                                                                                    aGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAccc                                          <  2:27239/120‑1 (MQ=255)                                                                                                    aGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAccc                                          >  1:27239/1‑120 (MQ=255)                                                                                                    aGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCATTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAata                       >  2:561282/1‑139 (MQ=255)                                                                                                                         ccAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCtggt  >  2:951539/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                ATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGT  >  NC_000913/3456963‑3457221

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