Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,773,941:1 +G 13.2% coding (945/1053 nt) calA → coniferyl alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,773,9411.G13.2% 122.8 / 17.3 38coding (945/1053 nt)calAconiferyl alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (18/15);  new base G (2/3);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

TCGGCCTGATCGAACCGATCGACCCGGCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCTGATTGATCGGTTGCCTGTACGGGCCTCTTCGCGG  >  NC_002947/2773799‑2774081
                                                                                                                                               |                                                                                                                                            
tCGGCCTGATCGAACCGATCGACCCGGCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATc                                                                                                                                        <  1:68847/149‑1 (MQ=255)
 cgtccTGATCGAACCGATCGACCCGGCTTTGCATGCGGGCAATCTGGGGATGCAGGGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTGAGGATGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCg                                                                                                                                       <  8:331569/146‑1 (MQ=255)
      tGATCGAACCGATCGACCCGGCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGGATATCcagca                                                                                                                                   >  5:252570/1‑147 (MQ=255)
           gAACCGATCGACCCGGCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTgg                                                                                                                              <  4:175443/148‑1 (MQ=255)
             aCCGATCGACCCGGCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTCGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGCACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGATGGata                                                                                                                           >  2:326248/1‑149 (MQ=255)
                gATCGACCCGGCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATcc                                                                                                                        >  7:159571/1‑149 (MQ=255)
                     aCCCGGCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGGATATCCAGCATCTCGATATCGCAAcg                                                                                                                  >  8:295364/1‑148 (MQ=255)
                         ggCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGGATATCCAGCATCTCGATATCGCAAcggatg                                                                                                              <  7:295364/149‑6 (MQ=37)
                          gCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGGATATCCAGCATCTCGATATCGCAAcggatg                                                                                                              <  8:337990/149‑6 (MQ=38)
                           cATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAAc                                                                                                             <  7:350275/149‑1 (MQ=255)
                            aTTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGGATATCCAGCATCTCGATATCGCAAcggatgtc                                                                                                            <  6:252570/149‑8 (MQ=37)
                            aTTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAAc                                                                                                             >  6:329858/1‑148 (MQ=255)
                                cATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAgg                                                                                                         <  6:188028/148‑1 (MQ=255)
                                cATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGc                                                                                                        <  3:84900/149‑1 (MQ=255)
                                   gCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTa                                                                                                     >  6:187848/1‑149 (MQ=255)
                                    cGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTAt                                                                                                    >  5:14392/1‑149 (MQ=255)
                                        cAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGa                                                                                                                                            <  7:99755/105‑1 (MQ=255)
                                        cAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGa                                                                                                                                            >  8:99755/1‑105 (MQ=255)
                                               gTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGgatga                                                                                          <  8:159571/148‑1 (MQ=255)
                                                    gaAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGCCTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCg                                                                                    >  2:328026/1‑149 (MQ=255)
                                                    gaAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCg                                                                                    >  3:327655/1‑149 (MQ=255)
                                                     aaGCGCCGGGTGCTGGCGGGGTCGTTGATAGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGAGCCGAGCACGACATCCGTTGCGAAATCGAGATGCTGGATATCCGCAACATCAACAAGGCCTATGAGAGGATGATTGcctg                                                                                   >  2:305897/1‑147 (MQ=255)
                                                         gacGGGTGCTTGCGGGGTCTTTGATGGGGGGGATTGCCGAGACTCATGAAGTGCTGGACTTCTGCGCCGAGCATGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGAt                                                                               <  1:305897/147‑1 (MQ=255)
                                                                    gCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTAc                                                                      <  3:156086/147‑1 (MQ=255)
                                                                                 tCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGa                                                        >  6:147252/1‑148 (MQ=255)
                                                                                 tCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGa                                                        >  6:145171/1‑148 (MQ=255)
                                                                                          ttGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACg                                              <  6:14392/149‑1 (MQ=255)
                                                                                                    tCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCt                                    >  4:142393/1‑149 (MQ=255)
                                                                                                       ggAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCTGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCtgat                                 <  1:326248/149‑1 (MQ=255)
                                                                                                       ggAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCtgat                                 >  1:280462/1‑149 (MQ=255)
                                                                                                             gCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCTGATTGATCg                           >  2:90038/1‑149 (MQ=255)
                                                                                                               tGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCTGATTGATCGGt                         >  2:95905/1‑149 (MQ=255)
                                                                                                                    ttCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCTGATTGATCGGTTGCCt                    >  4:190337/1‑149 (MQ=255)
                                                                                                                          gcCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCTGATTGATCGGTTGCCTGTACgg              <  5:329858/149‑1 (MQ=255)
                                                                                                                                 aCGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCTGATTGATCGGTTGCCTGTACGGGCCTCtt       >  8:212966/1‑149 (MQ=255)
                                                                                                                                      aTCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCTGATTGATCGGTTGCCTGTACGGGCCTCTTCGCgg  <  5:187848/149‑1 (MQ=255)
                                                                                                                                        ccGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATgg                                                   >  5:72930/1‑98 (MQ=255)
                                                                                                                                        ccGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATgg                                                   <  6:72930/98‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
TCGGCCTGATCGAACCGATCGACCCGGCATTGCATGCGGCCAATCTGGTGATGAAGCGCCGGGTGCTGGCGGGGTCGTTGATCGGGGGGATTGCCGAGACTCAGGAAGTGCTGGACTTCTGCGCCGAGCACGACATCCGTTGCGATATCGAGATGCTGGATATCCGCAACATCAACCAGGCCTATGAGCGGATGATTGCCGGTGATGTGAAGTACCGGTTCGTGATCGACATGGCCACGCTACAGGGCTGATTGATCGGTTGCCTGTACGGGCCTCTTCGCGG  >  NC_002947/2773799‑2774081

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