Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,672,639 A→C 8.4% Q198P (CAG→CCG)  hcaD → putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,672,6390AC 8.4% 106.7 / 2.2 44Q198P (CAG→CCG) hcaDputative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base C (2/2);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

TGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAG  >  NC_000913/2672499‑2672775
                                                                                                                                             |                                                                                                                                        
tGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAgg                                                                                                                                        <  1:907397/143‑1 (MQ=255)
    ggAACTATTCGTCTGGAACTCACTGCCAGCGCCACGCAGCGCAGATGTAACCTGACAGTGCATGACCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGACAGTGCAACGCGACCTTTTCCAGCGCCACCAGAAGGcccg                                                                                                                                    >  2:98293/1‑140 (MQ=255)
      aaCTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGgtg                                                                                                                                  >  1:1182865/1‑143 (MQ=255)
         ttttGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCCTGCTGgtggcg                                                                                                                               <  2:747620/141‑4 (MQ=255)
         tATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTAGCTTTTACAGCGCCACCACCAGGCTGGTGTgc                                                                                                                               >  2:1185357/1‑143 (MQ=255)
          aTTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCACACCAGCCTGCTGgtggcg                                                                                                                               <  1:968234/143‑4 (MQ=255)
            tGGTCTGGAACTGGCTGCCAGCGCCACGCGGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCa                                                                                                                            <  1:445626/143‑1 (MQ=255)
                 tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCCGGTGGGCGCATTCTg                                                                                                                       >  2:304458/1‑143 (MQ=255)
                      cgggcggcCAGCGCCACGCAGCGCGGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaa                                                                                                                  <  2:1220075/141‑1 (MQ=255)
                         gCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGCGTGCGCACTCCGCTCaataa                                                                                                               >  2:499592/1‑143 (MQ=255)
                                 cgcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAGTAATGCCATTg                                                                                                       >  1:138895/1‑143 (MQ=255)
                                    cACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACCGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAACGCCATTGAAc                                                                                                    >  1:948487/1‑143 (MQ=255)
                                          gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTgg                                                                                              <  1:1116116/143‑1 (MQ=255)
                                              agATGTAAGGTGACGGTGATTGGACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGa                                                                                          <  2:1152493/143‑1 (MQ=255)
                                              agATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACGACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAGTGCCATTGACCCTGTGGTCGa                                                                                          >  1:732398/1‑143 (MQ=255)
                                               gATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATCATGCCATTGAACCTGTGGTCGAt                                                                                         >  2:282811/1‑143 (MQ=255)
                                               gATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGAt                                                                                         <  2:927945/143‑1 (MQ=255)
                                                aTGTAAGGCGACAGTGACTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                        <  1:514843/143‑1 (MQ=255)
                                                        tGACAGTGATTGCACGGGCGGCAACCGTCACGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCCGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                <  2:813166/143‑1 (MQ=255)
                                                        tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                <  1:610185/143‑1 (MQ=255)
                                                        tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                <  1:398933/143‑1 (MQ=255)
                                                              tGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCTGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGCACACGGGGTCCATGGTGACCAAGTAGac                                                                          >  2:1034130/1‑142 (MQ=255)
                                                                  ggagcTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                      <  2:624944/139‑1 (MQ=255)
                                                                    aaCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGGGGTCGATGGTGAAAAAGTAGAACTGACg                                                                    >  1:329190/1‑143 (MQ=255)
                                                                      cTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGAc                                                                     >  1:209591/1‑140 (MQ=255)
                                                                      cTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGAc                                                                     <  2:209591/140‑1 (MQ=255)
                                                                      cTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCGAGCTTTGACAGCGCCACCAGCAGGCTGGCGTGGGCGTGCGGCGAAATAAGGCCGTCGAAGAGGTGGGTGATGGGGAAACAGTAGAAGTGACGCt                                                                  >  2:1067141/1‑143 (MQ=255)
                                                                         gcggcAACCGTCATGGGCCGTAATGCACCACCGACCGTGAAACGCTAGCTTTTACCGCGCCACCAGCCGGCTGGTGTGCGCATTCTGCTCGATAAAGCCGTTGAACAAGGGGCCGATGCTGAAACAGCAGAAAGGACGcagga                                                               >  2:913248/1‑139 (MQ=255)
                                                                               aCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                         <  1:756770/143‑1 (MQ=255)
                                                                                     aTGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCCCAACAATGCCATTGAACATGTGGTCGCTGGTGAACCAGTAGAACTGCCGCTGCACCGTCCGGAAAc                                                   >  1:295153/1‑143 (MQ=255)
                                                                                         gCCGTAATGCGCCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCAGTCTGCTCACTAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCtt                                               <  1:1045950/143‑1 (MQ=255)
                                                                                          ccGTAATGCACCGCCGCCCGTGCAACGCTATCTCTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACACTTc                                              <  2:1126734/143‑1 (MQ=255)
                                                                                           cGTAATGCAGCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACCTGTGGTCGATTGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                             <  1:541493/143‑1 (MQ=255)
                                                                                                               gCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                         <  1:1014186/143‑1 (MQ=255)
                                                                                                                 aaCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGt                       >  1:324281/1‑143 (MQ=255)
                                                                                                                  aCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTa                      <  1:724285/143‑1 (MQ=255)
                                                                                                                   cGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAt                     >  2:348110/1‑143 (MQ=255)
                                                                                                                     cTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGACCATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTg                   <  2:941432/143‑1 (MQ=255)
                                                                                                                        tCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCCTTCTGCTCAACAATGCCATTGAACAGGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTa                <  1:43474/143‑1 (MQ=255)
                                                                                                                          ttttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATc              <  1:953366/143‑1 (MQ=255)
                                                                                                                          ttttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATc              >  1:842654/1‑143 (MQ=255)
                                                                                                                              ccAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcg          <  1:965327/142‑1 (MQ=255)
                                                                                                                                    cacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAAc     >  1:216229/3‑143 (MQ=255)
                                                                                                                                      caccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGa   >  1:978637/1‑143 (MQ=255)
                                                                                                                                       tccggcagGGTGGGGTGCGGGTTCTGCGCAACAATGCCATTGGACATGTGGTCGATGGGGAAAAAGTAGAACTGGCGCTGCGAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAg  <  2:449572/139‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                        
TGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAG  >  NC_000913/2672499‑2672775

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