Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,927,735	Δ1 bp	coding (581/993 nt)	asnA →	asparagine synthetase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,927,734	0	G	.	100.0%	363.6 / NA	77	coding (580/993 nt)	asnA	asparagine synthetase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (31/46);  total (31/46)

CTGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTGGCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCACCCCGT  >  NC_000913/3927635‑3927827                                                                                                    |                                                                                              cTGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTg                                                                                               >  2:410279/1‑100 (MQ=255)   gCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTgcg                                                                                             >  2:870921/1‑98 (MQ=255)      ggATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCgg                                                                                         >  2:365684/1‑100 (MQ=255)       gATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGt                                                                                        <  1:329615/100‑1 (MQ=255)         tCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTAt                                                                                      <  2:1343430/100‑1 (MQ=255)          cAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTAtt                                                                                     >  2:1033329/1‑100 (MQ=255)          cAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGACGTGACCGGGCGATGGCGAAAGATCTT‑GCGCGGTAtt                                                                                     >  2:1066100/1‑100 (MQ=255)            gATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTcc                                                                                   >  1:670786/1‑100 (MQ=255)             aTCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCt                                                                                  <  1:1344154/100‑1 (MQ=255)               ccACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTg                                                                                >  2:732089/1‑100 (MQ=255)               ccACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTg                                                                                <  2:532002/100‑1 (MQ=255)               ccACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTg                                                                                >  2:249451/1‑100 (MQ=255)               ccACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTg                                                                                <  2:1341868/100‑1 (MQ=255)                  cTTCGTACACAGCCAGGATTTACTGCCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTc                                                                              <  1:1425865/99‑1 (MQ=255)                   ttCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCgg                                                                            <  2:143136/100‑1 (MQ=255)                    tCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGAGCTT‑GCGCGGTATTCCTTGTCggg                                                                           <  2:1436524/100‑1 (MQ=255)                     cGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGa                                                                          <  1:139451/100‑1 (MQ=255)                     cGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGAGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGa                                                                          >  1:500698/1‑100 (MQ=255)                       tACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGAtt                                                                        >  2:873437/1‑100 (MQ=255)                        acacAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTg                                                                       >  2:550511/1‑100 (MQ=255)                        acacAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTg                                                                       <  1:1269754/100‑1 (MQ=255)                         cacaGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGCAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTg                                                                       >  2:749146/1‑99 (MQ=255)                         cacaGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTgg                                                                      <  1:167709/100‑1 (MQ=255)                         cacaGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTgg                                                                      >  1:1571698/1‑100 (MQ=255)                           caGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGTGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTggcg                                                                    <  1:1106568/100‑1 (MQ=255)                             gCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTggcggc                                                                  >  2:692745/1‑100 (MQ=255)                             gCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTggcggc                                                                  <  2:1250438/100‑1 (MQ=255)                              ccAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTggcggta                                                                 <  1:1072563/100‑3 (MQ=255)                                   aGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTg                                                            <  2:1101764/100‑1 (MQ=255)                                   aGTTACTGTCTCGTTATCCGGATCTTGATGCCAAACGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGTTATTCCTTGTCGGGATTGGCGGCAAGCTg                                                            <  2:3204/100‑1 (MQ=255)                                       aCTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCg                                                        >  2:616500/1‑100 (MQ=255)                                          tccgatctATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATg                                                     <  1:962010/93‑1 (MQ=255)                                          gTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATg                                                     >  1:868192/1‑100 (MQ=255)                                           tctcGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATgg                                                    <  1:940657/100‑1 (MQ=255)                                            ctcGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATgg                                                    <  2:400170/99‑1 (MQ=255)                                            ctcGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGt                                                   >  1:1453979/1‑100 (MQ=255)                                              cGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGGTTGGCGGCAAGCTGAGCGATGGTCa                                                 >  1:1023081/1‑100 (MQ=255)                                              cGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCa                                                 >  2:273497/1‑100 (MQ=255)                                              cGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCa                                                 <  1:143408/100‑1 (MQ=255)                                                 tATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCg                                              <  2:371543/100‑1 (MQ=255)                                                  aTCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATgagatcgga                                             >  2:962005/1‑92 (MQ=255)                                                  aTCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGc                                             <  1:1093717/100‑1 (MQ=255)                                                    ccGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCa                                           >  2:50913/1‑100 (MQ=255)                                                    ccGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCa                                           >  2:253159/1‑100 (MQ=255)                                                     cGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCac                                          >  1:1367314/1‑100 (MQ=255)                                                     cGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCac                                          <  1:110613/100‑1 (MQ=255)                                                     cGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCac                                          <  1:870926/100‑1 (MQ=255)                                                       gATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCacga                                        >  1:730920/1‑100 (MQ=255)                                                         tCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCacgacg                                      <  2:1242765/100‑1 (MQ=255)                                                         tCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCacgacg                                      <  2:105163/100‑1 (MQ=255)                                                           ttGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTg                                    <  1:1033335/100‑1 (MQ=255)                                                           ttGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTg                                    >  1:361993/1‑100 (MQ=255)                                                             gATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTgcg                                  <  1:271103/100‑1 (MQ=255)                                                              aTGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTgcgc                                 <  1:1053440/100‑1 (MQ=255)                                                              aTGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATGCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTgcgc                                 <  2:1154302/100‑1 (MQ=255)                                                                 ccAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTgcgcgc                               >  2:597758/1‑99 (MQ=255)                                                                 ccAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCa                              <  1:692747/100‑1 (MQ=255)                                                                  caagtggCGTGAGCGGGCGATAGCGAAAGAGCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCAc                             <  1:129392/95‑1 (MQ=255)                                                                  cAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGAGATGGTCATCGCCACGACGTGCGCGCAc                             >  1:941831/1‑100 (MQ=255)                                                                   aaaGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTGATCGCCACGACGTGCGCGCAcc                            <  2:315274/100‑1 (MQ=255)                                                                       ggCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGAt                        <  2:975486/100‑1 (MQ=255)                                                                       ggCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGAt                        <  1:905789/100‑1 (MQ=255)                                                                        gCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGAtt                       >  1:124108/1‑100 (MQ=255)                                                                         cGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTa                      >  2:643861/1‑100 (MQ=255)                                                                           tGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTatg                    <  1:686347/100‑1 (MQ=255)                                                                           tGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTatg                    <  2:310560/100‑1 (MQ=255)                                                                            gAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTatga                   <  1:732091/100‑1 (MQ=255)                                                                            gAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTatga                   >  2:750623/1‑100 (MQ=255)                                                                            gAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTatga                   >  2:1121715/1‑100 (MQ=255)                                                                             aGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTatgat                  <  2:1292719/100‑1 (MQ=255)                                                                               cGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTAtgatga                <  2:500697/100‑1 (MQ=255)                                                                                 ggCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACt              <  1:365684/100‑1 (MQ=255)                                                                                   cGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTgg            <  2:730918/100‑1 (MQ=255)                                                                                   cGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTgg            <  1:1539911/100‑1 (MQ=255)                                                                                    gATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGTCGTGCGCGCACCGGATTATGATGACTGGa           <  1:1142005/100‑1 (MQ=255)                                                                                     aTAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAg          <  2:582904/100‑1 (MQ=255)                                                                                      tAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGc         >  2:1276745/1‑100 (MQ=255)                                                                                       aGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCa        >  2:1258971/1‑100 (MQ=255)                                                                                             aGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCACCCCGt  <  1:307436/100‑1 (MQ=255)                                                                                                    |                                                                                              CTGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTGGCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCACCCCGT  >  NC_000913/3927635‑3927827

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