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	.	97.3%	329.0 / 2.7	73	coding (580/993 nt)	asnA	asparagine synthetase A
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (44/27);  total (46/27)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CTGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTGGCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCACCCCGTC  >  NC_000913/3927635‑3927828                                                                                                    |                                                                                               cTGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTg                                                                                                <  1:835938/100‑1 (MQ=255) cTGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTg                                                                                                >  1:201167/1‑100 (MQ=255)  tGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTgc                                                                                               >  1:579253/1‑99 (MQ=255)  tGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTgc                                                                                               >  1:902225/1‑99 (MQ=255)  tGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTg                                                                                                >  1:97099/1‑99 (MQ=255)      ggATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCgg                                                                                          <  2:510862/100‑1 (MQ=255)          cAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTAtt                                                                                      <  1:1469269/100‑1 (MQ=255)          cAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGCGCGATAGCGAAAGATCTT‑GCGCGGTAtt                                                                                      >  1:1206743/1‑100 (MQ=255)           aGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTc                                                                                     >  2:1191186/1‑100 (MQ=255)           aGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTc                                                                                     >  1:1846613/1‑100 (MQ=255)            gATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTcc                                                                                    <  1:844533/100‑1 (MQ=255)             aTCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCt                                                                                   <  1:858977/100‑1 (MQ=255)              tCCACTTCGTCCACAGCCAGGAGTTACTGTCGCGTTATCGGGATCATCATGCCAAAGGGGATGGGGGGGCGATCGCCAAGGATCATgc                                                                                               >  1:1679524/1‑87 (MQ=255)                 aCTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTc                                                                               >  2:792959/1‑100 (MQ=255)                 aCTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTc                                                                               >  1:324350/1‑100 (MQ=255)                   ttCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCgg                                                                             >  1:616464/1‑100 (MQ=255)                      gTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGAt                                                                          <  2:666418/100‑1 (MQ=255)                      gTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGAt                                                                          >  2:1581112/1‑100 (MQ=255)                       tACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGAtt                                                                         >  2:1670813/1‑100 (MQ=255)                       tACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGAtt                                                                         >  2:1701001/1‑100 (MQ=255)                         cacaGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTgg                                                                       >  2:1868838/1‑100 (MQ=255)                         cacaGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTgg                                                                       >  2:144312/1‑100 (MQ=255)                         cacaGACAGGAGTTACTGTCTCGTTATCCGGATCTTTATGCCAAAGGTCGCTTGCGTGCGATAGCGATAGATCTT‑GCGTGTTATTCGTTGTCGCGTTTgg                                                                       >  1:1663381/1‑100 (MQ=255)                           caGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTggcg                                                                     <  1:703782/100‑1 (MQ=255)                            aGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTggcgg                                                                    <  2:1613001/100‑1 (MQ=255)                             gCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTggcggc                                                                   >  1:619082/1‑100 (MQ=255)                                aGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAg                                                                <  2:1553638/100‑1 (MQ=255)                                 ggAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGc                                                               <  2:1728157/100‑1 (MQ=255)                                    gTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGa                                                            >  1:1858096/1‑100 (MQ=255)                                      tACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGc                                                          >  2:159421/1‑100 (MQ=255)                                      tACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGc                                                          <  2:380984/100‑1 (MQ=255)                                       aCTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCg                                                         >  1:975268/1‑100 (MQ=255)                                       aCTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCg                                                         <  2:324352/100‑1 (MQ=255)                                        cTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGa                                                        >  2:1566076/1‑100 (MQ=255)                                          gTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATg                                                      >  2:1217437/1‑100 (MQ=255)                                          gTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATg                                                      >  2:42570/1‑100 (MQ=255)                                          gTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATg                                                      >  2:944710/1‑100 (MQ=255)                                          gTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATg                                                      >  1:621526/1‑100 (MQ=255)                                             tcGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTc                                                   >  2:747916/1‑100 (MQ=255)                                              cGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCa                                                  >  1:125949/1‑100 (MQ=255)                                              cGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGCCAAGCTGAGCGATGGTCa                                                  >  1:1766391/1‑100 (MQ=255)                                               gTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCAt                                                 >  2:1714544/1‑100 (MQ=255)                                                 tATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCg                                               <  2:139171/100‑1 (MQ=255)                                                  aTCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCTTCAAGCTGAGCGATGGTCATCg                                               <  1:405004/99‑1 (MQ=255)                                                   tCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGcc                                             >  1:1070191/1‑100 (MQ=255)                                                   tCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGcc                                             <  1:1179193/100‑1 (MQ=255)                                                   tCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGcc                                             <  1:168267/100‑1 (MQ=255)                                                    ccGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTGGCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGcc                                             >  1:211685/1‑100 (MQ=255)                                                    ccGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCa                                            >  1:515921/1‑100 (MQ=255)                                                    ccGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGACGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCa                                            >  1:649280/1‑100 (MQ=255)                                                     cGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCac                                           >  2:559464/1‑100 (MQ=255)                                                     cGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCac                                           <  1:636627/100‑1 (MQ=255)                                                     cGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCac                                           >  1:400319/1‑100 (MQ=255)                                                     cGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCac                                           >  2:239602/1‑100 (MQ=255)                                                      ggATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCacg                                          >  1:1599289/1‑100 (MQ=255)                                                        aTCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCacgac                                        >  1:565556/1‑100 (MQ=255)                                                              ntGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTgcgc                                  >  1:1151396/2‑100 (MQ=255)                                                              aTGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTgcgc                                  >  2:1176822/1‑100 (MQ=255)                                                                gCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTgcgcgc                                <  2:201168/100‑1 (MQ=255)                                                                gCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTgc                                    >  1:109877/1‑96 (MQ=255)                                                                   aaaGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGTCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCAcc                             >  1:637440/1‑100 (MQ=255)                                                                    aaGGGCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCg                            <  1:1873194/100‑1 (MQ=255)                                                                       ggCGTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGAt                         <  2:1844775/100‑1 (MQ=255)                                                                          gTGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTat                      <  1:42570/100‑1 (MQ=255)                                                                           tGAGCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTatg                     >  1:590861/1‑100 (MQ=255)                                                                              gCGGGCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTatgatg                  <  2:1846625/100‑1 (MQ=255)                                                                                 ggCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACt               <  2:1295545/100‑1 (MQ=255)                                                                                  gCGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACtn              <  1:338946/100‑2 (MQ=255)                                                                                   cGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTgg             <  2:234459/100‑1 (MQ=255)                                                                                   cGATAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTgg             <  2:619083/100‑1 (MQ=255)                                                                                     aTAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAg           >  1:117344/1‑100 (MQ=255)                                                                                     aTAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAg           <  2:1232212/100‑1 (MQ=255)                                                                                      tAGCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGc          <  2:1223743/100‑1 (MQ=255)                                                                                        gCGAAAGATCTTGGCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCa         >  1:535703/1‑100 (MQ=255)                                                                                        gCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCAc        >  2:351960/1‑100 (MQ=255)                                                                                        gCGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCAc        >  2:183432/1‑100 (MQ=255)                                                                                         cGAAAGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCAcc       >  1:347247/1‑100 (MQ=255)                                                                                           aaaGATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCAcccc     >  1:1166681/1‑100 (MQ=255)                                                                                              gATCTT‑GCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCACCCCGTc  >  2:27006/1‑100 (MQ=255)                                                                                                    |                                                                                               CTGCCGGATCAGATCCACTTCGTACACAGCCAGGAGTTACTGTCTCGTTATCCGGATCTTGATGCCAAAGGGCGTGAGCGGGCGATAGCGAAAGATCTTGGCGCGGTATTCCTTGTCGGGATTGGCGGCAAGCTGAGCGATGGTCATCGCCACGACGTGCGCGCACCGGATTATGATGACTGGAGCACCCCGTC  >  NC_000913/3927635‑3927828

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