Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,519,159	A→T	T45S (ACA→TCA)	curA →	curcumin/dihydrocurcumin reductase, NADPH‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,519,159	0	A	T	98.7%	263.0 / ‑3.9	78	T45S (ACA→TCA)	curA	curcumin/dihydrocurcumin reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (33/44);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTCGGCG  >  NC_000913/1519063‑1519255                                                                                                 |                                                                                                 cTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAGGGACAGGTGTTACTGCGCTCAg                                                                                               <  1:1384485/100‑1 (MQ=255) cTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAg                                                                                               >  2:230812/1‑100 (MQ=255)  tGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGt                                                                                              <  2:1126316/100‑1 (MQ=255)  tGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGt                                                                                              >  2:1219661/1‑100 (MQ=255)  tGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGt                                                                                              <  2:1331146/100‑1 (MQ=255)  tGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGt                                                                                              <  2:1331173/100‑1 (MQ=255)   ggCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGtt                                                                                             <  2:1214579/100‑1 (MQ=255)    gCCTCGCGTCGACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGTGTGAAGGACAGGTGTTACTGCGCTCAGttt                                                                                            <  2:1540434/100‑1 (MQ=255)    gCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGttt                                                                                            >  1:1031969/1‑100 (MQ=255)     ccTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTa                                                                                           <  2:251813/100‑1 (MQ=255)     ccTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTa                                                                                           >  2:616953/1‑100 (MQ=255)      cTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTAGAAGAAGATGATGTGGCCACACCGGGGGAAGGACAGGTGATAATGCGCTCAGGTTAt                                                                                          >  2:1421676/1‑100 (MQ=255)       tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAtt                                                                                         <  1:1368861/100‑1 (MQ=255)        cgcgTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAttt                                                                                        >  1:1823709/1‑100 (MQ=255)         gcgTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTg                                                                                       <  1:409592/100‑1 (MQ=255)             ccACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTccc                                                                                   <  1:1707963/100‑1 (MQ=255)               acaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTg                                                                                 >  2:1943033/1‑100 (MQ=255)                      gcACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGt                                                                          <  1:1331967/100‑1 (MQ=255)                      gcACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGt                                                                          >  1:2049977/1‑100 (MQ=255)                      gcACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGt                                                                          >  1:321501/1‑100 (MQ=255)                        aCCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGtat                                                                        <  1:635057/100‑1 (MQ=255)                             ttCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCg                                                                   <  1:1967339/100‑1 (MQ=255)                              tCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGt                                                                  <  1:1968524/100‑1 (MQ=255)                               ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                 <  2:1599316/100‑1 (MQ=255)                               ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                 >  2:952616/1‑100 (MQ=255)                                cGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTgg                                                                <  2:129844/100‑1 (MQ=255)                                 ggAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGa                                                               >  1:812901/1‑100 (MQ=255)                                   agaATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACg                                                             <  1:379607/100‑1 (MQ=255)                                   agaATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACg                                                             >  1:166264/1‑100 (MQ=255)                                     aaTTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTa                                                           <  2:1466456/100‑1 (MQ=255)                                       tttCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATg                                                         >  1:2317870/1‑100 (MQ=255)                                          ccGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGc                                                      >  2:355451/1‑100 (MQ=255)                                          ccGTCTGGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGc                                                      >  2:312843/1‑100 (MQ=255)                                           cGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCg                                                     <  2:2150737/100‑1 (MQ=255)                                            gTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGa                                                    >  2:51613/1‑100 (MQ=255)                                             tCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGAt                                                   <  2:1913798/100‑1 (MQ=255)                                             tCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGAt                                                   <  2:1031967/100‑1 (MQ=255)                                              cTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATg                                                  >  2:1310382/1‑100 (MQ=255)                                              cTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATg                                                  <  2:2170308/100‑1 (MQ=255)                                                 gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGc                                               >  1:1763444/1‑100 (MQ=255)                                                 gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGc                                               >  2:570499/1‑100 (MQ=255)                                                 gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTATTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGc                                               >  2:1443549/1‑100 (MQ=255)                                                  aagaagATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGcc                                              <  1:800334/100‑1 (MQ=255)                                                  aagaagATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGcc                                              >  2:1701621/1‑100 (MQ=255)                                                   agaagaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCa                                             >  1:1707609/1‑100 (MQ=255)                                                     aagaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATc                                           <  2:1552615/100‑1 (MQ=255)                                                      agaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCt                                          >  1:1891518/1‑100 (MQ=255)                                                      agaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCt                                          >  1:1901193/1‑100 (MQ=255)                                                       gatgatGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCCTCtt                                         <  2:1515131/100‑1 (MQ=255)                                                         tgatgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAt                                       <  2:1697490/100‑1 (MQ=255)                                                          gatgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAtt                                      >  2:1308753/1‑100 (MQ=255)                                                           atgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTc                                     <  2:876941/100‑1 (MQ=255)                                                            tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                    <  2:2221140/100‑1 (MQ=255)                                                            tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                    <  2:395652/100‑1 (MQ=255)                                                            tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                    >  1:133035/1‑100 (MQ=255)                                                                gCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCg                                 >  1:1001207/1‑99 (MQ=255)                                                                 ccACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGcc                               <  1:1366564/100‑1 (MQ=255)                                                                 ccACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGcc                               <  1:1443553/100‑1 (MQ=255)                                                                   acacCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTg                             >  2:1851625/1‑100 (MQ=255)                                                                   acacCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTg                             <  1:1842461/100‑1 (MQ=255)                                                                      ccGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTg                          <  2:1707604/100‑1 (MQ=255)                                                                       cGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGa                         <  2:1891511/100‑1 (MQ=255)                                                                        gggTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGat                        >  1:982701/1‑100 (MQ=255)                                                                        gggTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGat                        >  2:2104740/1‑100 (MQ=255)                                                                         ggTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGata                       <  1:1483725/100‑1 (MQ=255)                                                                             aaGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTgg                   >  1:62603/1‑100 (MQ=255)                                                                             aaGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTgg                   <  1:1219664/100‑1 (MQ=255)                                                                               ggACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcg                 >  1:1036697/1‑100 (MQ=255)                                                                               ggACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcg                 <  2:820123/100‑1 (MQ=255)                                                                               ggACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcg                 >  2:462766/1‑100 (MQ=255)                                                                                gACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcgg                <  1:51615/100‑1 (MQ=255)                                                                                   tggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGt             <  2:1483645/99‑1 (MQ=255)                                                                                    ggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGta            >  1:1483649/1‑99 (MQ=255)                                                                                     gtgtTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCATGa           <  1:1308756/100‑1 (MQ=255)                                                                                      tgtTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGAt          <  2:1870820/100‑1 (MQ=255)                                                                                       gtTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATg         <  2:1001205/100‑1 (MQ=255)                                                                                        ttACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATgg        <  2:583657/100‑1 (MQ=255)                                                                                              cgcTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcggcg  <  2:1763439/100‑1 (MQ=255)                                                                                                 |                                                                                                 CTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTCGGCG  >  NC_000913/1519063‑1519255

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