Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,574,062	A→T	C273S (TGC→AGC)	asd ←	aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,574,062	0	A	T	97.5%	274.7 / ‑1.8	79	C273S (TGC→AGC)	asd	aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (35/42);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCT  >  NC_000913/3573965‑3574156                                                                                                  |                                                                                               cttccgatctcccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTgc                                                                                              <  1:1915631/90‑1 (MQ=255)  cGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTgcg                                                                                             <  2:693718/100‑1 (MQ=255)   gACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTgcgc                                                                                            <  2:882065/100‑1 (MQ=255)    aCTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCa                                                                                           <  1:282653/100‑1 (MQ=255)     cTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTATTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCaa                                                                                          >  2:1328990/1‑100 (MQ=255)          gCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGccc                                                                                     >  1:1794012/1‑100 (MQ=255)          gCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGccc                                                                                     >  1:1509228/1‑100 (MQ=255)          gCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGccc                                                                                     <  1:1611342/100‑1 (MQ=255)           cccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTgcagatcggaag                                                                                    >  2:1915613/1‑90 (MQ=255)            ccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCg                                                                                   <  2:1042780/100‑1 (MQ=255)            ccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCg                                                                                   <  2:1427042/100‑1 (MQ=255)             cACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGTATGCCTGGCTGTGGCTGCGCAATGCCCCGa                                                                                  >  2:2020880/1‑100 (MQ=255)              aCGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGac                                                                                 <  2:1280728/100‑1 (MQ=255)                 gATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACg                                                                              >  1:278730/1‑100 (MQ=255)                   ttGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGca                                                                            <  2:295577/100‑1 (MQ=255)                   ttGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGca                                                                            <  1:1276287/100‑1 (MQ=255)                     gtgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcaca                                                                          <  1:1931004/100‑1 (MQ=255)                       gcgcAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGcacaca                                                                        >  1:293053/1‑100 (MQ=255)                        cgcAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACAt                                                                       >  1:688350/1‑100 (MQ=255)                         gcAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTTTTTTTATTATTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATa                                                                      >  2:1093466/1‑100 (MQ=255)                         gcAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATa                                                                      <  2:234006/100‑1 (MQ=255)                            gccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAAc                                                                   >  2:1675086/1‑100 (MQ=255)                                gcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATc                                                               >  2:623551/1‑100 (MQ=255)                                   gTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTAc                                                            >  1:1711773/1‑100 (MQ=255)                                    ttcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTAcc                                                           <  2:1742421/100‑1 (MQ=255)                                    ttcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTAcc                                                           <  1:672573/100‑1 (MQ=255)                                     tcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCg                                                          >  1:1044991/1‑100 (MQ=255)                                     tcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCg                                                          <  2:35374/100‑1 (MQ=255)                                      cttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCgg                                                         >  2:1536848/1‑100 (MQ=255)                                        tcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGaa                                                       >  1:815784/1‑100 (MQ=255)                                         ccACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGaa                                                       >  2:1204860/1‑99 (MQ=255)                                           aCGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTa                                                    <  1:1243671/100‑1 (MQ=255)                                             ggTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAAttcc                                                   >  2:1708173/1‑97 (MQ=255)                                                cGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGaa                                               <  2:497075/100‑1 (MQ=255)                                                cGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGaa                                               <  2:497390/100‑1 (MQ=255)                                                 ggAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAg                                              >  1:124721/1‑100 (MQ=255)                                                  gAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAg                                              >  2:977002/1‑99 (MQ=255)                                                   aaTAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAt                                            >  2:1148796/1‑100 (MQ=255)                                                      aGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtgtg                                         >  2:411032/1‑100 (MQ=255)                                                      aGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtgtg                                         >  2:462744/1‑100 (MQ=255)                                                       gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATgtgt                                        >  1:442079/1‑100 (MQ=255)                                                        tcttccgatctttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGtt                                       <  2:1771058/90‑1 (MQ=255)                                                           caTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAg                                    >  1:697206/1‑100 (MQ=255)                                                            aTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAgg                                   <  2:876377/100‑1 (MQ=255)                                                            aTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAgg                                   <  1:1186176/100‑1 (MQ=255)                                                             tCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggtttatgtg                                  <  1:1328999/100‑10 (MQ=25)                                                             tCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGa                                  <  2:1509219/100‑1 (MQ=255)                                                               ttttttttttttATTAGTGAATGCCTGGCTGTGGCTGCGCAATGCCTCGACACGCACACATAAACCATCTACCGGAATTACGGATGTTGTTTTTAGGAt                                 >  2:1855463/10‑99 (MQ=255)                                                               ttttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATc                                <  2:278729/100‑1 (MQ=255)                                                                tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                               <  2:169161/100‑1 (MQ=255)                                                                tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                               <  1:1597994/100‑1 (MQ=255)                                                                tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                               <  2:293052/100‑1 (MQ=255)                                                                tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                               <  1:1711956/100‑1 (MQ=255)                                                                tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                               <  2:483311/100‑1 (MQ=255)                                                                tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                               <  2:773612/100‑1 (MQ=255)                                                                tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGAATCt                               <  2:508558/100‑1 (MQ=255)                                                                    tCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                           <  2:1044986/100‑1 (MQ=255)                                                                    tCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                           <  1:1148802/100‑1 (MQ=255)                                                                     cAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgttg                          >  2:133797/1‑100 (MQ=255)                                                                     cAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgttg                          >  1:1139960/1‑100 (MQ=255)                                                                     cAATTTAAGAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgttg                          <  1:924307/100‑1 (MQ=255)                                                                        tttAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGtt                       <  1:1610235/100‑1 (MQ=255)                                                                        tttAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGtt                       >  1:1553245/1‑100 (MQ=255)                                                                          tAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTc                     >  2:1540487/1‑100 (MQ=255)                                                                           aaTAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTcc                    <  2:141928/100‑1 (MQ=255)                                                                             tAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCTGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGc                  >  2:565684/1‑100 (MQ=255)                                                                             tAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGc                  >  1:1258461/1‑100 (MQ=255)                                                                              aGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGcc                 >  2:1203613/1‑100 (MQ=255)                                                                              aGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGcc                 <  1:34770/100‑1 (MQ=255)                                                                                tGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTg               <  2:502890/100‑1 (MQ=255)                                                                                tGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTg               <  2:2027191/100‑1 (MQ=255)                                                                                 gAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGc              >  2:814439/1‑100 (MQ=255)                                                                                       cTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTc        <  1:623553/100‑1 (MQ=255)                                                                                        tggctgtggctgCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTcc       <  2:1793996/100‑1 (MQ=255)                                                                                         ggCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCATTTCCa      >  2:1124323/1‑100 (MQ=255)                                                                                          gCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCAc     >  1:1295496/1‑100 (MQ=255)                                                                                          gCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCAc     >  1:930400/1‑100 (MQ=255)                                                                                           cTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGCAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCAct    >  2:1136618/1‑100 (MQ=255)                                                                                             gtgGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCActct  <  2:815781/100‑1 (MQ=255)                                                                                             gtgGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCActct  <  2:917812/100‑1 (MQ=255)                                                                                                  |                                                                                               ACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCT  >  NC_000913/3573965‑3574156

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