Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,800,669 A→T 6.2% W66R (TGG→AGG)  PP_3358 ← hydroxycinnamoyl‑CoA hydratase‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,800,6690AT 6.2% 545.3 / 17.7 162W66R (TGG→AGG) PP_3358hydroxycinnamoyl‑CoA hydratase‑lyase
Reads supporting (aligned to +/- strand):  ref base A (72/80);  new base T (4/6);  total (76/86)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 163 total aligned reads displayed.

CCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTT  >  NC_002947/3800521‑3800818
                                                                                                                                                     |                                                                                                                                                     
ccGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAg                                                                                                                                                      <  1:4880730/150‑1 (MQ=255)
ccGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAg                                                                                                                                                      <  1:377230/150‑1 (MQ=255)
ccGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAg                                                                                                                                                      <  2:1260993/150‑1 (MQ=255)
  gTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGa                                                                                                                                                    >  1:3372826/1‑150 (MQ=255)
      aCCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGAccg                                                                                                                                                 >  1:1109634/1‑148 (MQ=255)
      aCCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCg                                                                                                                                                >  1:1618437/1‑150 (MQ=255)
       ccATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGc                                                                                                                                               <  2:5392932/150‑1 (MQ=255)
       ccATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGc                                                                                                                                               <  2:4915341/150‑1 (MQ=255)
       ccATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGc                                                                                                                                               <  1:4374494/150‑1 (MQ=255)
        cATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGcc                                                                                                                                              >  1:5188229/1‑150 (MQ=255)
           ggCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGACCGCCGGc                                                                                                                                           <  2:4074674/150‑1 (MQ=255)
               aTGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCg                                                                                                                                       <  1:3942269/150‑1 (MQ=255)
                tGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCGTGCCGGCGAGTCCTGGACCGCCGGCatgga                                                                                                                                      <  2:1489593/150‑5 (MQ=255)
                 ggTCGGTTTGTCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGt                                                                                                                                     <  1:1994882/150‑1 (MQ=255)
                  gTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTc                                                                                                                                    <  1:4296772/150‑1 (MQ=255)
                       tttGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGACCGCCGGCa                                                                                                                                          <  1:1957378/139‑1 (MQ=255)
                       tttGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGAGTCCTGGACCGCCGGCa                                                                                                                                          >  2:1957378/1‑139 (MQ=255)
                         tGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCa                                                                                                                             >  2:3369488/1‑150 (MQ=255)
                         tGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCa                                                                                                                             <  1:2856553/150‑1 (MQ=255)
                           gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGc                                                                                                                           >  2:730468/1‑150 (MQ=255)
                           gCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGc                                                                                                                           >  2:730469/1‑150 (MQ=255)
                              tACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACg                                                                                                                        >  2:4187833/1‑150 (MQ=255)
                                cAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGcc                                                                                                                      >  2:1424439/1‑150 (MQ=255)
                                 aGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCca                                                                                                                     >  2:2751042/1‑150 (MQ=255)
                                 aGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCca                                                                                                                     <  2:4233733/150‑1 (MQ=255)
                                   aCGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAATTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagc                                                                                                                   <  1:3152963/150‑1 (MQ=255)
                                   aCGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagc                                                                                                                   >  1:4767245/1‑150 (MQ=255)
                                   aCGCAGCAACTTCCATTACCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagc                                                                                                                   <  1:5222606/150‑1 (MQ=255)
                                     gcagcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgt                                                                                                                 >  1:3394713/1‑150 (MQ=255)
                                       agcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGTAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtca                                                                                                               <  1:517616/150‑1 (MQ=255)
                                       agcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtca                                                                                                               <  1:2825154/150‑1 (MQ=255)
                                       agcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtca                                                                                                               <  1:98727/150‑1 (MQ=255)
                                       agcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtca                                                                                                               >  1:4586529/1‑150 (MQ=255)
                                       agcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtca                                                                                                               >  1:4586564/1‑150 (MQ=255)
                                       agcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtca                                                                                                               <  1:965542/150‑1 (MQ=255)
                                        gcaACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtcag                                                                                                              >  2:2992252/1‑150 (MQ=255)
                                           aCTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCcagcgtcagcgt                                                                                                           >  2:5199745/1‑150 (MQ=255)
                                            cTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCagcgtcagcgtc                                                                                                          <  1:2474645/150‑1 (MQ=255)
                                             ttCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTcc                                                                                                         >  1:3288441/1‑150 (MQ=255)
                                              tCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCt                                                                                                        <  2:585741/150‑1 (MQ=255)
                                                 attgccattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCt                                                                                                     <  2:4244385/150‑1 (MQ=255)
                                                     ccattgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGNTTTCCGGGCCGGCGTCCACCTCGCNGAAGTACTCCTTCAGGTCCATGCCNGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCNTCCTGCTCAAg                                                                                                 >  2:1241214/1‑150 (MQ=255)
                                                      cattgcGACGCTACGCGACGAATCTTTTCCTGGAGGATATCCGGGCCGGCCTCAACAACGCGGAAGAAATCCTTCAGGTCCATGCCGGCGGTCCAGGACTCTCAGGAACCCCACAATACCAGCACGCAAGCGTCAGCGTCCTGCTCAAgg                                                                                                >  2:822302/1‑150 (MQ=255)
                                                      cagtgcGAGGCGACGCGACGACTCTTGGCCTGGAGGATTTCCGGGCGGGCGTCCACCTCGCGGAAGTACTGCTTCAGTTCATAGCCGGCGGTCCAGGACTCGCCGGCACCGGTCATTACCATCACGCCAGCGTCAGCGTCCTGCTCAAgg                                                                                                <  1:3923350/147‑1 (MQ=255)
                                                         tgcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGtt                                                                                             <  2:1089702/150‑1 (MQ=255)
                                                          gcGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGttt                                                                                            >  2:3794602/1‑150 (MQ=255)
                                                           cGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTc                                                                                           <  2:307024/150‑1 (MQ=255)
                                                           cGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTc                                                                                           <  2:525170/150‑1 (MQ=255)
                                                            gAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTcc                                                                                          >  2:1537670/1‑150 (MQ=255)
                                                             aGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCa                                                                                         <  2:3413604/150‑1 (MQ=255)
                                                              ggCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAg                                                                                        <  2:3952685/150‑1 (MQ=255)
                                                              ggCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAg                                                                                        <  2:1725854/150‑1 (MQ=255)
                                                                cTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCa                                                                                      <  1:3262669/150‑1 (MQ=255)
                                                                 ttCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCAc                                                                                     >  2:2024305/1‑150 (MQ=255)
                                                                    gcgACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTc                                                                                  >  1:266961/1‑150 (MQ=255)
                                                                      gacgaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGa                                                                                <  1:2177994/150‑1 (MQ=255)
                                                                      gacgaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGa                                                                                <  2:4736817/150‑1 (MQ=255)
                                                                         gaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAg                                                                                                 <  2:4852015/130‑1 (MQ=255)
                                                                         gaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCa                                                                             >  1:21072/1‑150 (MQ=255)
                                                                         gaATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCa                                                                             >  1:3624248/1‑150 (MQ=255)
                                                                           aTCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCAtt                                                                           <  2:3769973/150‑1 (MQ=255)
                                                                            tCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCAttt                                                                          >  1:3427998/1‑150 (MQ=255)
                                                                            tCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCAttt                                                                          >  1:84374/1‑150 (MQ=255)
                                                                               tttCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTccc                                                                       <  1:2732636/150‑1 (MQ=255)
                                                                                ttCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCg                                                                      >  2:5351786/1‑150 (MQ=255)
                                                                                    tGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTc                                                                  >  1:442619/1‑150 (MQ=255)
                                                                                       aggaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGc                                                                                       <  2:316911/126‑1 (MQ=255)
                                                                                         gaTTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGt                                                             <  2:131927/150‑1 (MQ=255)
                                                                                            ttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTggg                                                          <  1:1776495/150‑1 (MQ=255)
                                                                                            ttCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTggg                                                          >  2:330246/1‑150 (MQ=255)
                                                                                                gggCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTc                                                      >  1:1041766/1‑150 (MQ=255)
                                                                                                   ccGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATg                                                   <  1:1027496/150‑1 (MQ=255)
                                                                                                   ccGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATg                                                   <  2:3604717/150‑1 (MQ=255)
                                                                                                   ccGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATg                                                   <  2:2719627/150‑1 (MQ=255)
                                                                                                   ccGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATg                                                   >  2:1737403/1‑150 (MQ=255)
                                                                                                      gCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCa                                                >  1:5436786/1‑150 (MQ=255)
                                                                                                          ccACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGc                                            <  2:3456602/150‑1 (MQ=255)
                                                                                                          ccACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGc                                            <  2:5436786/150‑1 (MQ=255)
                                                                                                             ccTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCACCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGtt                                         <  2:3950112/150‑1 (MQ=255)
                                                                                                              cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                        >  1:5330281/1‑150 (MQ=255)
                                                                                                              cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                        >  2:2149648/1‑150 (MQ=255)
                                                                                                              cTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttt                                        <  2:5240759/150‑1 (MQ=255)
                                                                                                               tCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCAttt                                                                          >  1:831425/1‑115 (MQ=255)
                                                                                                               tCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGtttt                                       >  1:3965847/1‑150 (MQ=255)
                                                                                                                cgcgGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGttttt                                      >  2:2116740/1‑150 (MQ=255)
                                                                                                                     aaGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCggg                                 <  1:1163264/150‑1 (MQ=255)
                                                                                                                       gTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCg                               >  2:310813/1‑150 (MQ=255)
                                                                                                                        tACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGa                              >  2:4304502/1‑150 (MQ=255)
                                                                                                                           tCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTg                           <  1:4819383/150‑1 (MQ=255)
                                                                                                                            ccTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGa                          <  1:2416122/150‑1 (MQ=255)
                                                                                                                            ccTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGa                          >  2:2192696/1‑150 (MQ=255)
                                                                                                                             cTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCAc                                                                                     <  1:145889/90‑1 (MQ=255)
                                                                                                                             cTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAg                         <  1:1121829/150‑1 (MQ=255)
                                                                                                                                   gTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCAc                   >  2:5169239/1‑150 (MQ=255)
                                                                                                                                       aTGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAg               <  1:1565109/150‑1 (MQ=255)
                                                                                                                                          ccGGCGAGTCCTGGACCGCCGGCATCGGTCAATACCAGCACGCCAGCGTCCGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCg            >  1:3817511/1‑150 (MQ=255)
                                                                                                                                           cggcgagTCCTGGACCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCACGGCAGTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGa           >  2:3880614/7‑150 (MQ=255)
                                                                                                                                           cggcggTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGa           >  2:4852671/1‑150 (MQ=255)
                                                                                                                                            ggcggTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGAt          <  1:2646295/150‑1 (MQ=255)
                                                                                                                                              cggTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGc        <  2:1618437/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTTGACCATGGCGATGGTCGGTTTGGCATACAGACGCAGCAACTTCCATTGCCATTGCGAGGCTTCGCGACGAATCTTTTCCTGGAGGATTTCCGGGCCGGCGTCCACCTCGCGGAAGTACTCCTTCAGGTCCATGCCGGCGGTCCAGGACTCGCCGGCACCGGTCAATACCAGCACGCCAGCGTCAGCGTCCTGCTCAAGGGTTTCCAGCACGTCGACCATTTCCCGGTTCAGGGTGGGGCTCATGGCATTGCGTTTTTCCGGGCGATTGAGGGTCACCCAGGCGATGCCCGCTT  >  NC_002947/3800521‑3800818

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