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	96.0%	252.8 / 0.7	76	T45S (ACA→TCA)	curA	curcumin/dihydrocurcumin reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (39/34);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTCGGCG  >  NC_000913/1519061‑1519255                                                                                                   |                                                                                                 ttCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGcac                                                                                                 <  2:911557/100‑1 (MQ=255)   cTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAg                                                                                               <  1:167042/100‑1 (MQ=255)     ggCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGtt                                                                                             >  2:725305/1‑100 (MQ=255)      gCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGttt                                                                                            >  2:2004831/1‑100 (MQ=255)      gCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGttt                                                                                            >  1:71498/1‑100 (MQ=255)      gCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGttt                                                                                            >  1:1640301/1‑100 (MQ=255)      gCCTCGCGTCCACATGGCGCACCTGTTCCGGAAAATATCCGACTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAAGTGTTACTGCGCTCAGttt                                                                                            >  1:1843382/1‑100 (MQ=255)         tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGTGCTCAGTTTAtt                                                                                         <  1:1340021/100‑1 (MQ=255)         tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCTGTTTAtt                                                                                         <  2:1078403/100‑1 (MQ=255)         tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAtt                                                                                         <  2:1383263/100‑1 (MQ=255)         tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAtt                                                                                         <  1:249561/100‑1 (MQ=255)         tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCACCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAtt                                                                                         <  2:1646916/100‑1 (MQ=255)           gcgTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTg                                                                                       <  1:1349741/100‑1 (MQ=255)           gcgTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTg                                                                                       <  2:1285101/100‑1 (MQ=255)           gcgTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTg                                                                                       >  2:176542/1‑100 (MQ=255)             gTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTc                                                                                     >  1:97961/1‑100 (MQ=255)              tCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTcc                                                                                    <  2:570306/100‑1 (MQ=255)               ccACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTccc                                                                                   >  1:398545/1‑100 (MQ=255)               ccACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTccc                                                                                   <  1:176543/100‑1 (MQ=255)               ccACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTccc                                                                                   <  1:1522296/100‑1 (MQ=255)               ccACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTccc                                                                                   <  2:351590/100‑1 (MQ=255)                cacaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCt                                                                                  <  2:1318979/100‑1 (MQ=255)                cacaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCt                                                                                  >  1:1767996/1‑100 (MQ=255)                  caTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTgg                                                                                <  2:1640288/100‑1 (MQ=255)                   aTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGa                                                                               >  2:1597750/1‑100 (MQ=255)                       cgcACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGTCCCg                                                                           <  2:1398230/100‑1 (MQ=255)                       cgcACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCg                                                                           <  1:814517/100‑1 (MQ=255)                           ccTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGtata                                                                       >  1:1828299/1‑100 (MQ=255)                            cTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTatat                                                                      >  1:61704/1‑100 (MQ=255)                            cTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTatat                                                                      >  2:1455251/1‑100 (MQ=255)                              gTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGc                                                                    <  1:1052380/100‑1 (MQ=255)                                tCCGGAGAATTTCCGTCTTGAAGAAGATTATGTCTCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGAccc                                                                            <  2:1833934/90‑1 (MQ=255)                                tCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGt                                                                  >  2:1016102/1‑100 (MQ=255)                                 ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                 >  2:1107115/1‑100 (MQ=255)                                 ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                 >  2:767400/1‑100 (MQ=255)                                 ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                 <  1:739515/100‑1 (MQ=255)                                  cGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTgg                                                                >  1:1104981/1‑100 (MQ=255)                                   ggggaatcgccGTCTGGAAGAAGATTATGTCACCACTTTGGGTGAAGGATAGGTGTTACGGCGCACATTTTATTTGTCCCTGGACCCGTATATg                                                                     <  2:1893715/85‑1 (MQ=255)                                       aaTGTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTa                                                           <  2:1377674/100‑1 (MQ=255)                                         tttCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATg                                                         <  2:1162272/100‑1 (MQ=255)                                           tCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAg                                                       >  1:866244/1‑100 (MQ=255)                                             cGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCg                                                     >  1:130689/1‑100 (MQ=255)                                              gTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGa                                                    >  2:1671347/1‑100 (MQ=255)                                                cTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATg                                                  >  1:144058/1‑100 (MQ=255)                                                cTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATg                                                  <  2:398544/100‑1 (MQ=255)                                                  tGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAg                                                <  2:196610/100‑1 (MQ=255)                                                   gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAgn                                               <  2:965563/100‑2 (MQ=255)                                                     agaagaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCa                                             >  1:1640045/1‑100 (MQ=255)                                                     agaagaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCa                                             >  1:746061/1‑100 (MQ=255)                                                           tgatgTCGCCACTCCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAt                                       <  1:677519/100‑1 (MQ=255)                                                           tgatgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAt                                       <  2:400054/100‑1 (MQ=255)                                                             atgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTc                                     >  1:1491577/1‑100 (MQ=255)                                                              tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                    >  1:1446585/1‑100 (MQ=255)                                                              tgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCa                                    <  1:1016107/100‑1 (MQ=255)                                                               gTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCAc                                   >  2:50494/1‑100 (MQ=255)                                                               gTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCAc                                   >  2:1493669/1‑100 (MQ=255)                                                                 cGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCg                                 >  2:1883699/1‑100 (MQ=255)                                                                 cGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCg                                 <  1:725307/100‑1 (MQ=255)                                                                  gCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGATTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGc                                >  2:664236/1‑100 (MQ=255)                                                                   ccACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCg                                 >  2:524410/1‑98 (MQ=255)                                                                   ccACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGcc                               >  2:1423385/1‑100 (MQ=255)                                                                   ccACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGcc                               <  2:1491568/100‑1 (MQ=255)                                                                    cacaCCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCt                              <  2:1640032/100‑1 (MQ=255)                                                                    cacaCCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCt                              >  2:420489/1‑100 (MQ=255)                                                                        ccGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTg                          <  2:144057/100‑1 (MQ=255)                                                                         cGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGa                         >  2:2038804/1‑100 (MQ=255)                                                                           ggTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGata                       >  1:78941/1‑100 (MQ=255)                                                                            gTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGatat                      <  1:647879/100‑1 (MQ=255)                                                                               aaGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTgg                   <  2:514928/100‑1 (MQ=255)                                                                                  gACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcgg                >  2:416017/1‑100 (MQ=255)                                                                                      ggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTg            >  1:1516619/1‑100 (MQ=255)                                                                                      ggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTg            >  1:252361/1‑100 (MQ=255)                                                                                      ggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTg            >  1:525742/1‑100 (MQ=255)                                                                                            aCTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATAGTGGCGGCGTGATGGTc      >  2:1197249/1‑100 (MQ=255)                                                                                              tGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcgg    <  1:645600/100‑1 (MQ=255)                                                                                               ncgcTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcggc   >  1:1687733/2‑100 (MQ=255)                                                                                                cgcTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcggcg  <  1:1493678/100‑1 (MQ=255)                                                                                                   |                                                                                                 TTCTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTCGGCG  >  NC_000913/1519061‑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: ‑