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

CGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTG  >  NC_000913/3456975‑3457228                                                                                                                          |                                                                                                                                       cGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGaaaa                                                                                                                     >  1:354449/1‑139 (MQ=255) cGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGaaa                                                                                                                      <  1:826467/138‑1 (MQ=255) cGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGaaa                                                                                                                      >  2:826467/1‑138 (MQ=255)     aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGa                                                                                                                             <  2:611576/127‑1 (MQ=255)     aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGa                                                                                                                             >  1:611576/1‑127 (MQ=255)     aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                                 >  1:1118873/1‑139 (MQ=255)     aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                                 >  1:229499/1‑139 (MQ=255)     aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                                 >  1:912113/1‑139 (MQ=255)     aGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCa                                                                                                                 <  2:753317/139‑1 (MQ=255)           aGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                           <  1:880343/139‑1 (MQ=255)                     gagaTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACtt                                                                                                 <  2:1097840/139‑1 (MQ=255)                       gaTTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCa                                                                                               <  2:642189/139‑1 (MQ=255)                          tCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTa                                                                                            <  2:229499/139‑1 (MQ=255)                                  aaTTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGaaa                                                                                      <  1:70308/139‑1 (MQ=255)                                   aTTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                      <  2:568583/120‑1 (MQ=255)                                   aTTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                      >  1:568583/1‑120 (MQ=255)                                     tAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCt                                                                                 <  1:677748/139‑1 (MQ=255)                                       aTCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                               <  2:1114496/139‑1 (MQ=255)                                          aGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATg                                                                            >  2:1112533/1‑139 (MQ=255)                                              aaaGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGa                                                                        >  1:1103616/1‑139 (MQ=255)                                              aaaGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGa                                                                        >  2:1024559/1‑139 (MQ=255)                                               aaGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGaa                                                                       <  2:800388/139‑1 (MQ=255)                                                   cACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAagaga                                                                   <  2:912113/139‑1 (MQ=255)                                                     cGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                           <  2:984750/97‑1 (MQ=255)                                                     cGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                           >  1:984750/1‑97 (MQ=255)                                                      ggTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                                >  1:559580/1‑139 (MQ=255)                                                      ggTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                                <  1:208012/139‑1 (MQ=255)                                                      ggTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                                <  2:559580/139‑1 (MQ=255)                                                         aaaaGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagg                                                             >  2:384429/1‑139 (MQ=255)                                                         aaaaGCCAGATGCCAGCCCTC‑‑CACTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagg                                                             <  2:608555/139‑1 (MQ=255)                                                          aaaGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGagga                                                            <  1:1112533/139‑1 (MQ=255)                                                                  aTGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATa                                                    >  1:974156/1‑139 (MQ=255)                                                                  aTGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATa                                                    >  1:55604/1‑139 (MQ=255)                                                                  aTGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACGCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTCATGTCGAACCGTGCGAgtaaggaacaa                                                    >  1:1079230/1‑136 (MQ=255)                                                                    gCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAcc                                                  <  2:354449/139‑1 (MQ=255)                                                                     ccAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATAccc                                                 <  1:1024559/139‑1 (MQ=255)                                                                      cAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCg                                                <  2:848886/139‑1 (MQ=255)                                                                      cAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCg                                                <  1:257476/139‑1 (MQ=255)                                                                           ctc‑‑ctcTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGttt                                           <  2:874801/139‑1 (MQ=255)                                                                            tc‑‑ctcTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTtta                                          <  2:1103616/139‑1 (MQ=255)                                                                             c‑‑ctcTCCGCGAAGATTGTGGCGGATAAGTGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTttat                                         >  1:12316/1‑139 (MQ=255)                                                                                  ctcCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTAc                                      <  2:1118873/139‑1 (MQ=255)                                                                                            aTTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGttt                                           <  1:953174/124‑1 (MQ=255)                                                                                            aTTGTGGCGGATAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGttt                                           >  2:953174/1‑124 (MQ=255)                                                                                                      aTAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGaggaagga                                                        <  1:521410/101‑1 (MQ=255)                                                                                                      aTAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGaggaagga                                                        >  2:521410/1‑101 (MQ=255)                                                                                                      aTAAGCGAACCAACTCTCTTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCAc                  >  2:633949/1‑139 (MQ=255)                                                                                                                   ctctctTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGAt                                                                             >  1:52526/1‑67 (MQ=255)                                                                                                                   ctctctTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGAt                                                                             <  2:52526/67‑1 (MQ=255)                                                                                                                      tctTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTg  <  2:1079230/139‑1 (MQ=255)                                                                                                                      tctTGTCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTg  <  2:12316/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                       CGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTC‑‑CTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTG  >  NC_000913/3456975‑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: ‑