Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,822,334 G→A 9.5% intergenic (+19/‑94) PP_2475 → / → curA TetR family transcriptional regulator/NADPH‑dependent curcumin/dihydrocurcumin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,822,3340GA 9.5% 108.6 / 3.1 42intergenic (+19/‑94)PP_2475/curATetR family transcriptional regulator/NADPH‑dependent curcumin/dihydrocurcumin reductase
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base A (2/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGCGTCGTTCTGGCCTCACGTC  >  NC_002947/2822186‑2822476
                                                                                                                                                    |                                                                                                                                              
ccACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTCACAACATCTGAAGGCGAGGTTGCCCACCTg                                                                                                                                                >  1:126696/1‑149 (MQ=255)
            ggCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCAtt                                                                                                                                    >  4:15404/1‑149 (MQ=255)
               ccTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCAttttt                                                                                                                                 <  6:122552/149‑1 (MQ=255)
               ccTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCAACAGGTGGGCaa                                                                                                                                     >  5:347212/1‑144 (MQ=255)
               ccTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCAACAGGTGGGCaa                                                                                                                                     <  6:347212/145‑2 (MQ=255)
                    aCCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACAACTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCAc                                                                                                                            >  1:312812/1‑149 (MQ=255)
                       aGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGa                                                                                                                         <  2:312812/149‑1 (MQ=255)
                        gTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCTTTCTGACCTCTGACATCATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGAc                                                                                                                        >  3:303031/1‑149 (MQ=255)
                              gCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAAc                                                                                                                    >  2:145700/1‑147 (MQ=255)
                              gCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTATGACCTCTGACAACATCTGATGGCGAGGTTGCCCACCTGTTGGTTGGCATTTTTTTCACttccaaac                                                                                                                    >  2:174950/1‑140 (MQ=255)
                                tCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCAACAGGTGGGCaac                                                                                                                                    <  4:12735/129‑3 (MQ=255)
                                tCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCAACAGGTGGGCaac                                                                                                                                    >  3:12735/1‑127 (MQ=255)
                                        aGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGt                                                                                                        <  1:182620/149‑1 (MQ=255)
                                         gTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTa                                                                                                                  >  8:7235/1‑138 (MQ=255)
                                         gTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTa                                                                                                                  <  7:7235/138‑1 (MQ=255)
                                              ctgAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTaa                                                                                                    <  2:105217/147‑1 (MQ=255)
                                              ctgAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTa                                                                                                     <  7:87508/146‑1 (MQ=255)
                                                gAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCg                                                                                                >  1:206520/1‑149 (MQ=255)
                                                gAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCg                                                                                                >  1:212263/1‑149 (MQ=255)
                                                 aGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCgg                                                                                               <  2:212263/149‑1 (MQ=255)
                                                 aGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCgg                                                                                               <  2:206520/149‑1 (MQ=255)
                                                      gCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGAt                                                                                          >  3:14264/1‑149 (MQ=255)
                                                          cacCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATc                                                                                         >  8:260220/1‑146 (MQ=255)
                                                               tAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGACGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAg                                                                                 >  2:128032/1‑149 (MQ=255)
                                                                        gCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTc                                                                        <  3:15404/149‑1 (MQ=255)
                                                                         ccGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCa                                                                       <  4:27469/149‑1 (MQ=255)
                                                                                 aaaCGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAAccc                                                                >  1:115303/1‑148 (MQ=255)
                                                                                           gCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTg                                                      <  1:145700/148‑1 (MQ=255)
                                                                                                 aCTCAGACCATTCTGACCTCTGACAACATCTAATGGCGCGGTTGCCAACCTGTTGGTGGGCAATTTTTTCACTGACCAACTATTCGAATGGTCTAATTCGGGTGATCGAGAAAAGCCCAACTTCAGCAACCCCACGAGTTTGCACATGc                                               >  5:250542/1‑149 (MQ=255)
                                                                                                     aGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGcccc                                            >  3:308174/1‑148 (MQ=255)
                                                                                                     aGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATAGAGACAAGCCAAACTTCAGCAACCCCAGCAGTTTGCACATGcccc                                            >  4:190682/1‑148 (MQ=255)
                                                                                                      gACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCaa                                          >  5:203780/1‑149 (MQ=255)
                                                                                                                    ctGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCaaa                                     >  7:68095/1‑140 (MQ=255)
                                                                                                                    ctGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCaaa                                     <  8:68095/140‑1 (MQ=255)
                                                                                                                                gATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGcgtcgt                >  8:77740/1‑149 (MQ=255)
                                                                                                                                    gcgcGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGGCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGCGTCGTTCTg            <  1:174950/149‑1 (MQ=255)
                                                                                                                                      gcgGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGCGTCGTTCTGGc          <  7:260220/149‑1 (MQ=255)
                                                                                                                                      gcgGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGCGTCGTTCTGGc          <  1:128032/149‑1 (MQ=255)
                                                                                                                                      aaaaatGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGCGTCGTTCTGGc          >  3:177412/6‑149 (MQ=255)
                                                                                                                                          ttGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCACGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGCGTCGTTCTGTCCTCa      <  6:250542/149‑1 (MQ=255)
                                                                                                                                              ccACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGgcgc                     <  5:93185/130‑1 (MQ=255)
                                                                                                                                              ccACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGgcgc                     >  6:93185/1‑130 (MQ=255)
                                                                                                                                              ccACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGCGTCGTTCTGGCCTCACGTc  >  4:142239/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCACACTGCGACGGCCCTTTACCAGTTGTGGCTCGGTGCGAGTCTGCTGAGCAAGCTACACCGTAATAGTCGGCCGCTGGAAAACGCCATGGCCACAACTCAGACCATTCTGACCTCTGACAACATCTGATGGCGCGGTTGCCCACCTGTTGGTGGGCATTTTTTTCACTGACCAACTAGTCGACTGGTCTAATTCGGGTGATCGAGACAAGCCCAACTTCAGCAACCCCAGGAGTTTGCACATGCCCCAATCAAAGCTAATCAATCGGCGCGTCGTTCTGGCCTCACGTC  >  NC_002947/2822186‑2822476

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