Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,798,742:1 +A 13.9% coding (201/831 nt) PP_3358 ← hydroxycinnamoyl‑CoA hydratase‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,798,7421.A13.9% 114.9 / 15.1 36coding (201/831 nt)PP_3358hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base . (18/13);  new base A (2/3);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGA  >  NC_002947/3798599‑3798887
                                                                                                                                                |                                                                                                                                                 
ccGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAg                                                                                                                                             <  2:501306/150‑1 (MQ=255)
             cGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGACCGCCGGCat                                                                                                                                <  2:194815/150‑2 (MQ=255)
                tGGTCGGTTTGGCATACAAACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATACCGGCGGTCCAGGACTCGCCGGCACCgg                                                                                                                             >  1:627428/1‑150 (MQ=255)
                   tCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCa                                                                                                                          >  1:569636/1‑150 (MQ=255)
                       tttGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCCGGACTCGCCGGCACCGGTCAAGAc                                                                                                                      <  2:53354/150‑1 (MQ=255)
                         tGGCATACAGACGCATCAACTTCCCTTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGTAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCa                                                                                                                    <  2:159490/150‑1 (MQ=255)
                             aTACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCAc                                                                                                                >  2:801546/1‑150 (MQ=255)
                              tACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCcggcg                                                                                                                                                    <  1:335615/114‑1 (MQ=255)
                              tACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCcggcg                                                                                                                                                    >  2:335615/1‑114 (MQ=255)
                                            cTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATAc                                                                                                                      >  1:145568/1‑129 (MQ=255)
                                            cTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATAc                                                                                                                      <  2:145568/129‑1 (MQ=255)
                                            cTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCagcgtcagcgtc                                                                                                 >  2:72284/1‑150 (MQ=255)
                                                 attgccattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTGCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCAGCACCGGTCAATACCAGCACGCCGGCGTCAGCGTCCTGCt                                                                                            >  2:395211/1‑150 (MQ=255)
                                                     ccattgcGAGGCTTCGCGGCGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGACCGCCGGCat                                                                                                                                <  2:1211/110‑2 (MQ=255)
                                                     ccattgcGAGGCTTCGCGGCGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGACCGCCGGCat                                                                                                                                >  1:1211/1‑109 (MQ=255)
                                                       attgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAgg                                                                                       >  1:840412/1‑149 (MQ=255)
                                                       attgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAgg                                                                                       <  2:840412/149‑1 (MQ=255)
                                                             aGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCa                                                                                <  2:569636/150‑1 (MQ=255)
                                                                              ttttCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGACCGCCGGCat                                                                                                                                >  1:256710/1‑84 (MQ=37)
                                                                              ttttCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGACCGCCGGCat                                                                                                                                <  2:256710/85‑2 (MQ=37)
                                                                               tttCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTccc                                                              <  1:801546/150‑1 (MQ=255)
                                                                                             tCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTgggg                                                >  2:124785/1‑150 (MQ=255)
                                                                                                           cACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCg                                  >  2:777124/1‑150 (MQ=255)
                                                                                                              cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                               >  1:652394/1‑150 (MQ=255)
                                                                                                                    gAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCCGCGTCCGCGGCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                         <  2:627428/150‑1 (MQ=255)
                                                                                                                    gAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                         <  1:72284/150‑1 (MQ=255)
                                                                                                                         aCTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                      <  1:627547/100‑1 (MQ=255)
                                                                                                                         aCTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGAc                                                                      >  2:627547/1‑100 (MQ=255)
                                                                                                                             cTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                         <  1:96901/141‑1 (MQ=255)
                                                                                                                             cTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCgg                         >  2:96901/1‑141 (MQ=255)
                                                                                                                             cTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAg                >  2:417927/1‑150 (MQ=255)
                                                                                                                             cTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAg                >  1:5890/1‑150 (MQ=255)
                                                                                                                                   gTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTc                                                                >  2:534020/1‑96 (MQ=255)
                                                                                                                                   gTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTc                                                                <  1:534020/96‑1 (MQ=255)
                                                                                                                                    tCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCAcc         >  1:72629/1‑150 (MQ=255)
                                                                                                                                       aTGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAg                <  2:80121/140‑1 (MQ=255)
                                                                                                                                       aTGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAg                >  1:80121/1‑140 (MQ=255)
                                                                                                                                        tGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAgg     >  2:383782/1‑150 (MQ=255)
                                                                                                                                           cggcggTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGa  <  2:652394/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGA  >  NC_002947/3798599‑3798887

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