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_0009131,519,1590AT98.2% 190.8 / ‑3.6 57T45S (ACA→TCA) curAcurcumin/dihydrocurcumin reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (31/25);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

CTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTCGGC  >  NC_000913/1519063‑1519254
                                                                                                |                                                                                               
cTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAg                                                                                              >  2:577995/1‑100 (MQ=255)
   gCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGttt                                                                                           >  1:769036/1‑100 (MQ=255)
     cTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAt                                                                                         >  2:965109/1‑100 (MQ=255)
      tCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTAtt                                                                                        <  2:1281705/100‑1 (MQ=255)
          gTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTc                                                                                    >  2:901811/1‑100 (MQ=255)
           tCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTAATGCGCTCAGTTTATTTGt                                                                                     >  2:838748/1‑98 (MQ=255)
            ccACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTccc                                                                                  <  1:988459/100‑1 (MQ=255)
            ccACATGGCGCACCTGTTCCGGAGAATTTACGTCTTGAGGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTccc                                                                                  >  2:1279755/1‑100 (MQ=255)
             cacaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCt                                                                                 <  1:94693/100‑1 (MQ=255)
             cacaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCt                                                                                 >  1:10993/1‑100 (MQ=255)
             cacaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCt                                                                                 >  2:275624/1‑100 (MQ=255)
              acaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTg                                                                                >  2:639017/1‑100 (MQ=255)
              acaTGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTg                                                                                >  1:375333/1‑100 (MQ=255)
                    cgcACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCg                                                                          <  2:769031/100‑1 (MQ=255)
                    cgcACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCg                                                                          <  2:924981/100‑1 (MQ=255)
                      cACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGta                                                                        <  1:75506/100‑1 (MQ=255)
                        ccTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGtata                                                                      <  1:1044731/100‑1 (MQ=255)
                        ccTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGtata                                                                      <  1:1065847/100‑1 (MQ=255)
                         cTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTatat                                                                     <  2:1301518/100‑1 (MQ=255)
                              ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                >  1:1110863/1‑100 (MQ=255)
                              ccGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTg                                                                <  1:948678/100‑1 (MQ=255)
                               cGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTgg                                                               >  2:1166627/1‑100 (MQ=255)
                                 gagaATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGAc                                                             >  2:532887/1‑100 (MQ=255)
                                 gagaATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCGGGACCCGTATATGCGTGGAc                                                             <  2:798143/100‑1 (MQ=255)
                                      tttCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATg                                                        <  2:375330/100‑1 (MQ=255)
                                        tCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAg                                                      <  1:863569/100‑1 (MQ=255)
                                          cGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCg                                                    <  2:652564/100‑1 (MQ=255)
                                          cGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCg                                                    >  2:1100366/1‑100 (MQ=255)
                                                gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGc                                              <  1:893744/100‑1 (MQ=255)
                                                gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGc                                              >  1:285512/1‑100 (MQ=255)
                                                gaagaaGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGc                                              >  2:111626/1‑100 (MQ=255)
                                                   gaagaTGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCAt                                           >  2:379386/1‑100 (MQ=255)
                                                        tgatgTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAt                                      <  1:763792/100‑1 (MQ=255)
                                                        cttccgatctacacCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAt                                      <  1:93575/90‑1 (MQ=255)
                                                               gCCACACCGGGTAAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGc                               >  1:854848/1‑100 (MQ=255)
                                                                  acacCGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTAtagatcggaag                            >  2:93573/1‑90 (MQ=255)
                                                                      cGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGa                        >  2:1156636/1‑100 (MQ=255)
                                                                      cGGGTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGa                        >  2:274203/1‑100 (MQ=255)
                                                                       gggTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTg                         <  2:315959/98‑1 (MQ=255)
                                                                        ggTGAAGGACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGata                      <  1:1168347/100‑1 (MQ=255)
                                                                              ggACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcg                >  1:802680/1‑100 (MQ=255)
                                                                               gACAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTggcgg               >  1:690701/1‑100 (MQ=255)
                                                                                aCAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCCTCTTATTCACCGCCTGTTGATATTggcggc              >  2:324761/1‑100 (MQ=255)
                                                                                 cAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGgcggcg             >  2:871450/1‑100 (MQ=255)
                                                                                 cAGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGgcggcg             <  2:743429/100‑1 (MQ=255)
                                                                                  aGGTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGt            <  1:532891/100‑1 (MQ=255)
                                                                                   ggTGTTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTg           >  1:472424/1‑100 (MQ=255)
                                                                                      gtTACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATg        <  2:745635/100‑1 (MQ=255)
                                                                                      gtTACTGCGCTCAGTTTATTTGTCCCTGGACCCGGATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATg        <  2:538469/100‑1 (MQ=255)
                                                                                        tACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGt      <  1:1294617/100‑1 (MQ=255)
                                                                                        tACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGt      >  2:218368/1‑100 (MQ=255)
                                                                                        tACTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGt      >  2:1279874/1‑100 (MQ=255)
                                                                                         aCTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTc     <  1:577999/100‑1 (MQ=255)
                                                                                         aCTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTc     >  1:676118/1‑100 (MQ=255)
                                                                                          cTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcg    <  2:1077772/100‑1 (MQ=255)
                                                                                          cTGCGCTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcg    >  2:158855/1‑100 (MQ=255)
                                                                                            gcgcTCAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTcggc  >  2:1120463/1‑100 (MQ=255)
                                                                                                |                                                                                               
CTGGCCTCGCGTCCACATGGCGCACCTGTTCCGGAGAATTTCCGTCTTGAAGAAGATGATGTCGCCACACCGGGTGAAGGACAGGTGTTACTGCGCACAGTTTATTTGTCCCTGGACCCGTATATGCGTGGACGTATGAGCGATGAGCCATCTTATTCACCGCCTGTTGATATTGGCGGCGTGATGGTCGGC  >  NC_000913/1519063‑1519254

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