Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,330 C→T 53.0% intergenic (+79/‑117) atoC → / → atoD DNA‑binding transcriptional activator/ornithine decarboxylase inhibitor AtoC/acetyl‑CoA:acetoacetyl‑CoA transferase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,323,3300CT53.0% 4.5 / 41.0 44intergenic (+79/‑117)atoC/atoDDNA‑binding transcriptional activator/ornithine decarboxylase inhibitor AtoC/acetyl‑CoA:acetoacetyl‑CoA transferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (11/10);  new base T (9/14);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

CATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGAAAACAAAATTGATGACA  >  NC_000913/2323188‑2323467
                                                                                                                                              |                                                                                                                                         
cATCAGGCGCCGTGCATTGGTGTATACACTCCAGGACTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                                                                                                                           <  2:791827/143‑1 (MQ=255)
      tCGCCGGGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                                                                                                                     <  1:438657/143‑1 (MQ=255)
            tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCACTTTATCACt                                                                                                                               >  2:527364/1‑143 (MQ=255)
                  gATGTATAAACTACAGGAATACGGTATCGATCAGGCGGATGTATAACACTAAAACTTGATATGCAGAAATTTGCACAGTGCGCAATTTTCTGAATAGTCGATAATTTTTATTATAAATCACAATTAAACTTATACCGTAatat                                                                                                                         >  2:670401/1‑143 (MQ=255)
                    tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTGCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                                                                                                                       <  1:532360/143‑1 (MQ=255)
                    tGTATAAACTCCAGCAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                                                                                                                       <  1:542328/143‑1 (MQ=255)
                       ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAAt                                                                                                                    >  1:5222/1‑143 (MQ=255)
                        taAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATCTTCAAtt                                                                                                                   >  1:527119/1‑143 (MQ=255)
                               cAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGCTaaa                                                                                                            >  1:594336/1‑143 (MQ=255)
                                 ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                                                                                                          >  1:492999/1‑143 (MQ=255)
                                  gAATACGGTATCGATCCGGCGGATGTATAACACCGAAACTTGGTGTGCAGAAATGTGCACCGTGCGTACTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATaa                                                                                                         <  1:642756/143‑1 (MQ=255)
                                      aCGGTAGCGATGCGGCGGATGTATAACACCAAAAGTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                                                                                                     <  1:607518/143‑1 (MQ=255)
                                       cGGTATCGATCCGGCGGATGTCTAACTCCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaaa                                                                                                    <  1:733954/143‑1 (MQ=255)
                                          tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc                                                                                                 >  1:315603/1‑143 (MQ=255)
                                             cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat                                                                                              <  2:531008/143‑1 (MQ=255)
                                                tCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaataat                                                                                           <  2:336220/143‑1 (MQ=255)
                                                   ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTCGTTAAATAACTAAATCCAATAAtctc                                                                                        >  2:213036/1‑143 (MQ=255)
                                                      ggATGTATAGCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAACCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAtt                                                                                     <  1:780206/143‑1 (MQ=255)
                                                         tGTATCACATCAAATCTTGCTATGCAAAGATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAGTAGTTAAATAACTAAATCCAATAATCTCATTCTg                                                                                  <  2:437499/143‑1 (MQ=255)
                                                                  ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTcccc                                                                         >  2:526270/1‑143 (MQ=255)
                                                                    accccTTCCTATGCAGAACTCTGCACAGTGCGCCATTTTCTGCATAGCCGCTCATCCTCCTTACAAATCCCCATCCAATTTATCCCTTCACATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCtt                                                                       <  1:96074/139‑1 (MQ=255)
                                                                       tcTTGCTATGCAGAAATTTGCACAGTGCGCAATGTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCt                                                                    <  2:565098/142‑1 (MQ=255)
                                                                          ggCTATGCAGAAAATTGCGCACAGTGCAAATTTCTGCATAGCCGCTCATTGTCCTTATAAATCACCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATACAATAATCTCATTCTGGCACTCCGATTGCGAtt                                                                 >  2:573075/2‑143 (MQ=25)
                                                                           gCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTGAAATAACTAAATCCAATAATCTCATTATGGCACTCCCCTTGCTATTg                                                                >  2:565660/1‑143 (MQ=38)
                                                                                 cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGAc                                                          >  1:302418/1‑143 (MQ=255)
                                                                                  aGAAATCTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACt                                                         <  1:313615/143‑1 (MQ=255)
                                                                                        ttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTAccc                                                   >  1:472675/1‑143 (MQ=255)
                                                                                                  cgcAACTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATGAGTGAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGgtgt                                         <  2:167099/143‑1 (MQ=255)
                                                                                                   gcAATTTTCTGCATAGCCGCTCATTCTCCGTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTa                                        <  2:540953/143‑1 (MQ=255)
                                                                                                    cAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATATCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTAt                                       <  2:329469/143‑1 (MQ=255)
                                                                                                    cAATTGTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTAt                                       <  2:218459/143‑1 (MQ=255)
                                                                                                     aaTTGTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATg                                      <  1:592866/143‑1 (MQ=255)
                                                                                                          tCTGCATAGCCGCTCATTTTCCTTATAAATCCCCATCCAATTTATCCCTTCTTATTCAATTACTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAaga                                 <  1:414735/143‑1 (MQ=255)
                                                                                                                     gCTCATTCTCCTTATATCTCCGCATTCAACTTATCCCTTCATATGCAATTAGTTAAATAAGTAAATCCAATCATCTCACGCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATaaaaaa                      <  2:494641/143‑1 (MQ=255)
                                                                                                                         aTTCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGAGATAAAAAATgcc                  >  2:64268/1‑141 (MQ=255)
                                                                                                                         aTTCTCATTATCAATCTCCGTCCAATTTATCCCTTCATATTCAATTAGTGATATAAATAAAACTAATAGTCTCATTCTGGCACTCCCCTTTCGATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGaa                  <  2:202892/143‑1 (MQ=255)
                                                                                                                          ttCTCCTTATAAATCCCCATTCAATTTATCCCTTCATATCCAACTAGTTAAATAACTAAATCCAATAATCGCATTCTGGAACTCCCCGGGCTATTGCCTGACCGGACCGACATCGGGGGATGCAAGAGGGATAAAAAATGaaa                 <  2:810077/143‑1 (MQ=255)
                                                                                                                             tcCTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTAGAACTCGCCTTGCTATTGCATGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGAAAACa              >  2:715286/1‑143 (MQ=255)
                                                                                                                             tcCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGAAAACa              >  1:98687/1‑143 (MQ=255)
                                                                                                                              ccTTATAAATCCCCATTCAATTTATCCCTTCATATTCAATTAGGTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCCAGAGGGATAAAACATGCAAACaa             >  1:768097/1‑143 (MQ=255)
                                                                                                                                 tataAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAACAAAAGAAAACAAAAt          >  2:425821/1‑143 (MQ=255)
                                                                                                                                  ataAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGAAAACAAAAtt         >  1:778763/1‑143 (MQ=255)
                                                                                                                                    aaaTCACCATTCAATGTATCCCTTCATATTCAATTCGTTAAATAACTAAATCCAAGAATCTCATTCTGGCACTGCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGAAAACAAAATtga       <  2:592527/143‑1 (MQ=255)
                                                                                                                                     aaTCCCCATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGAAAACAAAATtgat      <  1:452007/143‑1 (MQ=255)
                                                                                                                                         cccATTCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGAAAACAAAATTGATGACa  >  1:825673/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                         
CATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAAGAGGGATAAAAAATGAAAACAAAATTGATGACA  >  NC_000913/2323188‑2323467

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