Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,456,682 C→A 9.3% R85S (CGC→AGC)  paaE → ring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,456,6820CA 9.3% 125.9 / 3.2 43R85S (CGC→AGC) paaEring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component
Reads supporting (aligned to +/- strand):  ref base C (18/21);  new base A (2/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GCCCCGGTCAACATTTGACCTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTATTACGCCAATGCTGGCGATTATCGCC  >  NC_000913/1456539‑1456825
                                                                                                                                               |                                                                                                                                               
gCCCCGGTCAACATTTTACCTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGaa                                                                                                                                            <  6:162356/149‑1 (MQ=255)
gCCCCGGTCAACATTTGACCTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGaa                                                                                                                                            >  1:32008/1‑149 (MQ=255)
gCCCCGGTCAACATTTGACCTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGaa                                                                                                                                            <  6:146491/149‑1 (MQ=255)
 ccccGGTCAACATTTGACCTTAAAAGCCAGCTTTGATGGCGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAac                                                                                                                                           <  7:121849/149‑1 (MQ=255)
  cccGGTCAACATTTGACCTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAaca                                                                                                                                          <  7:321553/149‑1 (MQ=255)
   ccGGTCAACATTTGACCTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAacac                                                                                                                                         <  8:295663/149‑1 (MQ=255)
                gACCTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTa                                                                                                                            >  7:35927/1‑149 (MQ=255)
                   cTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGATGTTTCTCCCGCTATGCCAGCGaa                                                                                                                                            <  1:334043/130‑1 (MQ=255)
                   cTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGATGTTTCTCCCGCTATGCCAGCGaa                                                                                                                                            >  6:221474/1‑130 (MQ=255)
                   cTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGATGTTTCTCCCGCTATGCCAGCGaa                                                                                                                                            >  2:334043/1‑130 (MQ=255)
                   cTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGATGTTTCTCCCGCTATGCCAGCGaa                                                                                                                                            <  5:221474/130‑1 (MQ=255)
                       aaaGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCt                                                                                                                     <  1:210553/149‑1 (MQ=255)
                        aaGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCGGCCAGGGTATGACGCTg                                                                                                                    <  5:287597/149‑1 (MQ=255)
                        aaGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTg                                                                                                                    >  5:163726/1‑149 (MQ=255)
                         aGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTgg                                                                                                                   <  1:149773/149‑1 (MQ=255)
                           ccAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAg                                                                                                                 >  7:80925/1‑149 (MQ=255)
                            cAGCTTTGATGGTGAGGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAgg                                                                                                                >  7:232957/1‑149 (MQ=255)
                                             aaTTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAgg                                                                                                 >  1:122470/1‑147 (MQ=255)
                                             aaTTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAgg                                                                                                 <  2:122470/147‑1 (MQ=255)
                                              aTTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGc                                                                                              <  5:32496/149‑1 (MQ=255)
                                                  cggcGTTGTTACGCCATTTGCCGCAGCTATCTGCCTGGCTAAAATCGTGTGGCGGTGAAAGGTATTGAAGGCGGTCGTTTCTGCCGCTATGCCCGCGAACAGATCCGCGAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCAttt                                                                                          <  1:278144/146‑1 (MQ=255)
                                                    cctttggtACTCCATTAGCCGCAGCTATCTGACTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAATGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGTTCATGGTGCCGCAGGGGCATTTCg                                                                                        <  6:32008/146‑1 (MQ=255)
                                                    ccGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCg                                                                                        >  1:68277/1‑149 (MQ=255)
                                                     cGTTGTTACTCCATTTGCCGCAGCTATCTGACTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCgg                                                                                       >  8:2485/1‑149 (MQ=255)
                                                            aCTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAg                                                                                >  3:68452/1‑149 (MQ=255)
                                                                cATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCg                                                                             >  1:64622/1‑148 (MQ=255)
                                                                  tttGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAg                                                                          <  2:32008/149‑1 (MQ=255)
                                                                     gCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGc                                                                        <  8:35927/148‑1 (MQ=255)
                                                                             tATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGc                                                                 >  2:249529/1‑147 (MQ=255)
                                                                              aTCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAg                                                              >  7:162982/1‑149 (MQ=255)
                                                                                   ccTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAgg                                                             <  4:113062/145‑1 (MQ=255)
                                                                                        cGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATc                                                    <  8:80925/149‑1 (MQ=255)
                                                                                         gAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCt                                                   <  8:162982/149‑1 (MQ=255)
                                                                                                 tgtgGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATg                                                                                                                          >  7:233055/1‑70 (MQ=255)
                                                                                                      cGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTgcag                                        >  7:90375/1‑147 (MQ=255)
                                                                                                                cATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTa                            >  3:238447/1‑149 (MQ=255)
                                                                                                                  ttGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTAt                           >  2:310685/1‑148 (MQ=255)
                                                                                                                    gAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTATTa                         <  2:68277/148‑1 (MQ=255)
                                                                                                                      aGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTATTACGc                      <  6:163726/149‑1 (MQ=255)
                                                                                                                               tttCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTATTACGCCAATGCTgg             >  3:245240/1‑149 (MQ=255)
                                                                                                                                     ccGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCAGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTATTACGCCAATGCTGGCGa          <  8:232957/146‑1 (MQ=255)
                                                                                                                                       gCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTATTACGCCAATGCTGGCGATTATc     >  3:85116/1‑149 (MQ=255)
                                                                                                                                           tGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTATTACGCCAATGCTGGCGATTATCGcc  <  8:90375/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GCCCCGGTCAACATTTGACCTTAAAAGCCAGCTTTGATGGTGAAGAATTACGCCGTTGTTACTCCATTTGCCGCAGCTATCTGCCTGGCGAAATTAGTGTGGCGGTGAAAGCCATTGAAGGCGGACGTTTCTCCCGCTATGCCCGCGAACACATCCGCCAGGGTATGACGCTGGAGGTCATGGTGCCGCAGGGGCATTTCGGCTATCAGCCGCAGGCCGAACGCCAGGGGCGCTATCTGGCAATTGCAGCAGGATCAGGTATTACGCCAATGCTGGCGATTATCGCC  >  NC_000913/1456539‑1456825

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